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BEHAVIOUR OF PLASMA IN A ROTATING MAGNETIC FIELD 


A. LeGatTow1cz 
INSTITUTE OF NUCLEAR RESEARCH 


Warsaw, PoLAND 


An investigation is made of the non-relativistic motion of a charged particle (in a plasma) in 
an external rotating electromagnetic field of the form: Hx =H, cos (my/c) cos wt, Hy 
cos (wx/c) sin wot, Hz = constant, Ey = Ey= 0, Ez= H, [sin (@x/c) cos wt + sin (wy/c) sin of), 
with wx/e<1 and wy/c<1. If the condition — 1 < (eHz/mew) <—1 + (eH,/mcw)? is satisfied, then 
the particle moves away from the Z-axis. The particle energy is ~ H,2e?r?/(mc?k) in °K where 7 is 
the average distance from the Z-axis. 

The motion of the plasma is then investigated taking into account its proper electromagnetic 
field. The following transport equation is used: nm +VV-v=nqg (E+ vx H/c) —Vy—nmVq+p. 
The assumptions are: a plasma consisting of equal populations of electrons and deuterons with 
density ~10'5/em*, < 10!/sec, H,~10° G, v<0.1¢c. At t=0 it is assumed that =10% °K 
and v=0 and that the derivatives of the plasma density with respect to the space variables are 
negligible compared to other terms of the transport equation. An expansion in powers of v/c is 
used. Zeroth, first and second order approximations are calculated using the Laplace transformation. 
Up to the second order of approximation, the field causes neither a durable change in plasma 
density nor a charge separation. 

Oscillations of four different frequencies appear in the plasma. At a definite frequency of the 
rotating field there appears a resonance phenomenon in which the amplitude of oscillation in- 
creases linearly with ¢. At resonance the mean energy per ion (in °K) transferred directly to the 
ionic part of the plasma increases with time as follows: (1/192 =?) (e?m/c4+M*k) (H,* H.?/n,?) @. 
This means, for example, that in an axial field of 104G and a rotating field of amplitude 10°G, 
the time necessary to provide energy corresponding to 10° °K (disregarding losses) is ~0.3 sec. 


1. Introduction 


To obtain controlled fusion it is necessary to have 
very hot plasma; besides, the plasma must be con- 
fined for sufficiently long time within a limited space. 
Practically, the only way to satisfy these requirements 
is to provide adequate magnetic fields. Therefore it 
is important to know the behaviour of plasma in 
particular configurations of magnetic fields. 

One possible configuration might be a rotating 
magnetic field in the presence of a constant homo- 
geneous magnetic field perpendicular to the former. 
A preliminary investigation of the behaviour of plasma 
in this configuration of fields, under certain simplifying 
assumptions, is undertaken in the present paper. 


2. External magnetic field 


We assume that the external magnetic field has 
the following form 


Hx = Hy cos (wy/c) cos wt 
Hy=H, cos (w2/c) sin wt 
H.=constant (1) 
E,=E,y=0 


H, [sin (w2/c) cos wt + sin(@y/c) sin wf). 


Near the Z-axis this field corresponds to the con- 
figuration of the rotating magnetic field: 


cosmt 
H,=H,sinot 
H.=constant (2) 
E,=E,=0 

E.=(H, o/c) (x cosmt+-ysin wt), 


provided the following conditions are satisfied: 
(3) 


3. The motion of a charged particle in the rotating field 


As an initial problem let us consider the motion 
of a charged particle in the field defined by Eq. (2). 
Starting from the non-relativistic equation of motion: 

mf =eE + (e/e)txH, (4) 


we obtain, after taking into consideration Eq. (2), 
the equations for the separate components: 


d@x dy dz 
“dr? =2; ($2, sin dr 
ay dx dz 
dr 2, dr (2g cos dr (5) 


=2, (x sin T—y cosT), 
where 


t=ot, 2,=e Hz/(mew). 
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Integrating the third equation of Eq. (5) we get 


=, (x sint — y cost) +C;. (6) 
Substituting Eq. (6) into the first two equations of 
Eq. (5), we obtain a system of equations that can be 
written thus 


dU 1 


i234, 20? U+ 2 U *e2it + 


(7) 
where +iy. 
Substituting then 
U=feit, (8) 
we obtain an 7 with constant coefficients 


(9) 


Solving Eq. (9) and substituting the result into 
Eq. (8), we obtain 


+ MCs les (10) 


where 


m=—(2+2))[ 5 (2 T+) 


2.C, 


sin (Q, t+ g,) + Cz sin (Q,7 + 


— 1 
1/2 [2 + Q,2 + + 22,4 VA]? 
A=(24 2,2 + 2,7 + 22,)?—4 (Q, + 1) (Q3 + 1— Q,?). 


C,, Cy, Cs, Py, Pg are constants dependent on the 
initial state of the particle. 

The condition for the particle to move within a 
limited distance from the Z-axis is that {2,,, must 
be real. This leads to the inequality 


+ 2,2 + 2.2, + 2)? >4 (Q, + 1) (234+ 1—2,2)>0. 
(11) 


The left side of the inequality is always satisfied. 
The right side leads to a condition which can be 
written as 


either 
e H.|/(mew)<—1 
or 
e (12) 


If there are two kinds of particles, negative electrons 
and positive ions, by a proper choice of H:, H, and 
® we can obtain the four possible cases, viz: 
(a) both kinds of particles move within a limited 
distance from the Z-axis; 


(b) both kinds of particles move away from the 
Z-axis; 
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(ec, d) the particles of one kind move away and 
those of the other kind remain confined. 


The energy achieved by a particle can be estimated 
on the basis of Eq. (6). Assuming that the average 
distance from the Z-axis is r, the average energy of 
the particle will be of the order 


Tk = mek (K), (13) 


where k denotes the Boltzmann constant. If 
H,, = 1000 G and r = 10 em, the energy of the deuterons 
(M =3,3 -10-*4g) will be of the order 


Tx ~5.107(°K). 


4. The motion of plasma taking into account its proper 
electromagnetic field 


4.1. INITIAL EQUATIONS AND ASSUMPTIONS 


Let us start from the transport equation for char- 
ged particles, page 97 of [1], 


nm (or -~+vV- -y) =ng(E+**<*) +p. 


(14) 
We assume that: 


1) The electron density n and the ion density N 
are of the order of 1015 

2) The rotating field frequency m< 10" sec~!. 

3) The amplitude of the external magnetic field 
is of the order of 10° G. 

4) Plasma temperature at t=0 is of the order of 
10° °K; this temperature has been achieved by 
other means. 

5) Particle velocities are at least one order of 
magnitude lower than the velocity of light. 

6) Deuterons serve as ions. 


Let us introduce new variables 
w=vi/c_ (electrons) 
W=Vi/e (ions) 
Op=V4 7 M|(m+ M) (15) 


XL, Y, > 


=(/c) xe, T= . 


N is the plasma density at t=0. 

Accepting the above assumptions, we note that 
the first terms of the left and right side of Eq. (14) 
are the most important because the other terms 
are smaller by at least three orders of magnitude. 
In the notation of Eq. (15) the transport equations 
for the electrons and ions will be, respectively: 


ow 


Cr Wp (E+wxH) 

(16) 


To the above equations the equation of continuity 
and the Maxwell equations should be added. In the 
new variables of Eq. (15) the equations are 


F 
\ 
| 
| 
| 
Le 
wa 
4 


ity 
he 


On o 
“Os + =p 
oN +2 V.(NW)=0 
ér Op 


= e(N— n), j=ec(NW— nw) 


228 (17) 
CT 
VxH =— E+ 
@ 
V-E= ese 0 


V-H=0. 


The system, Eqs. (16) and (17) is nonlinear. To solve 
we expand all variables in powers of a small para- 


meter 6 

w=)'w.d', W=)'Wd 

H= 

n == N=)'N, 

i= 


To determine 6 we assume that w,=W,=0 and 
w, is of the order of unity — then the parameter 6 
will be of the order of v/c for electrons. 


(18) 


4.2. THE INITIAL CONDITIONS 
We assume that at ¢=0 the plasma is at rest, i.e. 
w (0) =WO)=j@=09, 
and is neutral 
n(0)= NO), 0(0)=0. 


We assume also that at ¢=0 the plasma is practically 
homogeneous; this means that the derivatives of 
plasma density with respect to spatial variables, are 
negligible in comparison with the other terms of the 
equations. 


4.3. ZERO-ORDER APPROXIMATION 


Substituting the expansion, Eq. (18), into Eqs. (16), 
(17), and comparing the coefficients for the zero power 
of 6 we obtain a zero-order approximation: 


+ V (m9 Wo)=9, (Ny — no) 

V(NoWo) =9. jo=ec (NoWo— 
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and the Maxwell equation for the quantities with 
the index 1=0. From the assumption w,= W,= 0 
it follows that E,=0 and H, does not depend on 
time. We are going to assume that the homogeneous 
field along the X,-(or Z-) axis is of zero order and is 
an order of magnitude larger than the amplitude 
of the rotating field. Thus we obtain in the zero- 
order approximation the following solution: 


E,=0 w,=0 W,=0 j,=—0 o,=0 
20) 
H, (0.0.H;). H,=constant, n )—N,=constant, 
where H, is 2,- or z-component of the external 
magnetic field, see Eq. (2). 


4.4. FIRST-ORDER APPROXIMATION 


Proceeding in the same way as in Section 4.3, and 
taking into account Eq. (20) we obtain the equations 
of the first-order approximation 


=— (E, + Ww, Hy) 

ow, (E, + W, x H,) 
+ = V-(m)W,)=9 
+ =0 j,—ecn,(W,—w,), 


and the Maxwell equations for the quantities with 
index 

Consider the electromagnetic field E and H to be 
a superposition of the external field E,, H, expressed 
by Eq. (1) and of the proper plasma field ¢, h, 


S=e+E,, H=h+H,. (22) 


Assuming, further, that the proper plasma field 
at ¢=0 is equal to zero, we can obtain the solution 
of the system of equations, Eq. (21), without any 
difficulty after appropriate calculations using the 
Laplace transformation : 


Ey, = Ey.=0, Hy3=0, Wy = vy =0 
ju=)x=9, 
Hy (. 
E\3= (sin &, cos 2, 7+ o, &, sin 
H, (| 1 
cos &, cos 2, t + a, é,) 
Hy, (2 
(5- cos &, sin 2,1) 
eH, (23) 
dmeapQ, 
[sin &, sin 2,7 — sin &, (cos 2,7 — 1)| 
1 
eH, 
dmc opQ, 
[sin &, sin Q, 7 — sin (cos 2,7 — )| 
1 
Q="-, 
Op 
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4.5. SECOND-ORDER APPROXIMATION 


In the second-order approximation, having taken 
account of the solutions of the zero and first order, 
we obtain the equations 


+ - aX € 2 + b) 


Ot me MC Wp 
ow, e 
Or M cop 


QV: (ng W2)=0, (Nz — m2) 
ON 
and the Maxwell equations for quantities with the 
index |= 2. In this approximation 

Weg= &3=0, (25) 
which results from the fact that 

b,=B,=0. 

Applying the Laplace transformation to Eq. (24) and 
introducing complex variables: 


W= Wy, + 1 W= Wo, +i Woe 


-+2V-(N,W,)=0, cn, (W,—w.), 


= +1 Egg 


J=Jat+ije 
eH, _ eH, _ 
Meop 


~ Meap 


1 ~ ~ ~ 


4x. 
=P ij 


(+ B) 
(27) 


where p is a transformation parameter, the asterisk 
denotes a conjugate variable, # is the transform of 
w, ete., and the subscript symbol / denotes @/0é). 

From the Eq. (27) we can obtain the equation for 7 
which can be written in the form 


l ~ 1 ~ 1 ~ “~** > 


where 


p+ i Qie 
mM P (p+ iQ) (p—i2) 
eH, 
ic = (M — m) ape 


(28) 
2 (Pp + (p — +1 

(p — (p — iQe) 
G=iDH 


dc op, 


[sin &, sin 2, — sin (cos 2,7 — 1)| ‘ 
1 
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Retaining in the expressions for D and H only the 
largest terms (at least two orders of magnitude larger 
than the remaining ones), we obtain 


1 


G= iF [sin (€, + &) + sin (§, — &,)] p+? (29) 


Assuming 7) in the form 


= F, sin (€, + &) + F, sin (&, — &), 


and substituting it into Eq. (28), we obtain, after 
proper calculations, 


where 
Q=2 2 [ (p? + 93?) (p? + Q.?) + p? + Qi Qe] 
+ (p? + (p? + Qo?) + 2 (p? + Qj Qe) + 1] 
= 2. [p? + —Qie) 
+ (2,0, + Qie — Qie Qe) P— 26 Qe] (31) 


a= 


- [p? + (Qe— Qi + Qie) 
+ (QQe+ Qie Qi — Die Qe) QQ. 


(Qie + Qe — Qi) P— Qe QB Qe} 
{p* + (Qj Qe + Qj Qie — Qc Qic) p+ Qie} | 


Now we substitute the expressions for PF, and F, 
into Eq. (30) and obtain the Laplace transform of 
the electric field 7 and the transform of electron 
and ion velocities of the plasma: 


[ 8, sin (&, + &) + B,sin(é,—é,) | 


P 1 


1 


i (32) 


+8; 


After the application of the inverse transformation 
we obtain the solution of the second-order approxi- 
mation. However, as such expressions would be too 
complicated and awkward for discussion we shall 
discuss the transform instead, Eqs. (30)—(32), and 
apply the inverse transformation only in such cases 
where it is necessary to obtain numerical results. 


5. Discussion 
5.1. THE PROBLEM OF THE CHANGE OF PLASMA DENSITY 


From the expressions Eq. (31) and (32) it follows 
that in the velocities of both plasma components 
the periodical terms as well as the non-periodical 
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are present, the latter being independent of time. From 
the point of view of a durable change in plasma density 
the latter are the most important. Introducing p=0 
into Eq. (31) and (32) we find that the time-independent 
terms of the plasma velocities satisfy the equations 


V-w.=0, V-W,=0. (33) 


This means, that up to the second order of approx- 
imation, the field concerned causes neither a durable 
change in plasma density nor a charge separation. 


5.2. FREQUENCIES OF PLASMA OSCILLATION 


The frequencies plasma oscillations are given by 
the singular points of the expressions, Eq. (32). For 
electrons these will be the points 


p= = i 2. 
(34) 
p= +iQ,, 
and the singular points of F, and P,. 
For ions the points will be 
p= +i1Q 
(35) 
p= + i22,, 


and the singular points of F, and F,. 

The latter are found by the assumption Q=0 in 
Eq. (31). This expression is of the third degree in p?; 
therefore, zero points of this polynomial can be found 
without difficulty. We shall not perform this computa- 
tion because it is complicated and not interesting; 
we shall confine ourselves to the statement that, on 
the basis of the form of the expression for Q, it can 
be shown that all the three roots are imaginary and 
that each root is different. This means that in the 
plasma there appear oscillations of three additional 
different frequencies 


p=+iQ, 
p=+iQ, (36) 
p=+iQ,. 


5.3. RESONANCE PHENOMENA 


There arises an interesting question whether it is 
possible to find a frequency of the rotating field such 
as to cause a resonance of the ion part of the plasma. 

Putting p?= — 2;? in the expression for Q in 
Eq. (31), we obtain for the resonance frequency of 
the rotating field the formula 


Q2= Qj [1 + (Be — 


2 (Q¢e — Q) (37) 


With the accepted assumptions this frequency in 
ordinary units is approximately equal to 


Ore meer (38) 


It follows from the above formulae that the reson- 
ance frequency of the rotating field is not equal 


PLASMA IN A ROTATING MAGNETIC FIELD 


to the cyclotron-frequency of the ions but has another 
value, practically independent of H, (=Hz). This 
value is of the order of the proper frequency of the 
ion part of the plasma. 

Also from Eq. (37) it follows that, for plasma in a 
dispersed state when the condition 


is satisfied, the resonance frequency of the rotating 
field is approximately equal to the ion cyclotron 
frequency, 


Qe ~ 


Thus, in this case we have an ordinary cyclotron 
resonance phenomenon. 


5.4. THE PROBLEM OF ENERGY EXCHANGE BETWEEN 
THE ELECTROMAGNETIC FIELD AND THE PLASMA 


Power density transferred to plasma by the electro- 
magnetic field is expressed by E-j. Taking into 
consideration the power density transferred to the 
ion part of the plasma and taking into account only 
the resonance terms we get from Eqs. (30), (31), (32), 
an expression for average energy transferred per single 
ion, in kinetic temperature units: 

1 em H,'H,? 

T= +92 ng? (e K). (40) 

If the axial field Hz = H, = 10‘ G and the amplitude 
of the rotating magnetic field H,=10* G the time 
necessary to provide energy corresponding to 10° °K, 
disregarding losses, amounts to about 0,3 sec. 


6. Final remarks 


All the above results have of course an approximate 
character because the assumptions take into account 
only the volume effects. In practical equipment 
plasma is never homogeneous. In addition, the assumed 
initial conditions disregard the way of switching in 
the rotating field. 
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AXIAL CONDUCTION AND RADIATION LOSSES IN A STABILISED LINEAR PINCH 


A. H. pE Borpe anp F. A. Haas*, 


Evectric Co. Lrp., NELSON RESEARCH LABORATORIES, 


STAFFORD, ENGLAND 


This paper presents a theory of conduction and radiation losses in a linear pinched discharge 
under steady conditions. The model selected is one in which the ohmic heating in a thin skin 
of current is equated to the radiation and axial conduction losses, the discharge being considered 
at a uniform pressure determined by the Bennett relation modified to include the effect of a 
possible completely trapped axial magnetic field. Formulae for temperature and other physical 
quantities are presented and limiting forms considered in a variety of circumstances. The effect 
of thermoelectric phenomena is considered and conditions under which the treatment is likely 


to be applicable discussed. 


1. Introduction 


This paper investigates theoretically the significance 
of thermal conduction loss at the electrodes of a 
stabilised linear pinched discharge. BRaGINSkIT and 
SHAFRANOV [1], and Harnes [2] have also considered 
this problem but these authors have investigated 
models differing considerably from that discussed here. 

The present work has been developed from that of 
KavurMaNn and Furtu [3]. They consider a quiescent 
plasma in which axial conduction loss through the 
electrodes is maintained by joule heating. Adopting 
a one-dimensional formulation and assuming the 
electrodes to be held steady at some known tempera- 
ture, they derive a relation connecting voltage across 
the discharge and maximum temperature generated. 
They also obtain the temperature distribution as a 
function of distance from the centre of the discharge. 

Under suitable conditions, the axial current in a 
discharge is at first confined by inductive effects to 
a thin layer at the surface of the discharge. According 
to TAYLER [4], when the current has penetrated 
to an effective depth of about 15% of the radius 
of the discharge, instabilities cause break-up and 
the gas is no longer adequately confined. Within the 
main body of the discharge the pressure should be 
approximately uniform since in this region there is 
no appreciable current. These considerations suggest 
a simple model which is susceptible to mathematical 
treatment and based on the following assumptions: 

(1) The gas of the discharge is fully ionised 
throughout a cylinder, but the axial current 
only passes through a skin at the surface 
occupying a fraction A of the cross section of 
the cylinder. This assumption replaces the 
consideration of inductive effects. The general 
problem of how the axial current is distributed 
in a cylindrical conductor has been considered 
by [2]. 

(2) The current in the skin is determined by the 
usual resistivity formula. 


(3) The pressure is uniform throughout the 
discharge. 

(4) Any longitudinal magnetic field is considered 
to be completely trapped. This assumption 
is justified to some extent by the work of 
ALLIBONE ef al. [5]. 

(5) The discontinuity in the axial magnetic field, 

at the surface of the cylinder, gives rise to 

an azimuthal sheet current which we assume 
is uniform within a skin. This skin occupies 

a fraction yu of the cross section of the discharge. 

The only loss processes operative are axial 

thermal conduction and radiation losses through 

the emission of bremsstrahlung. 

Sheath effects close to the electrodes are 

neglected. 

(8) The temperature of the plasma is uniform 
across any cross section of the discharge per- 
pendicular to the axis but varies along the 
axis of the tube. 

(9) The ion and electron temperatures are equal. 
(10) The discharge has reached a steady state. 
(11) Temperature and current are related by the 

pinch relation of BENNETT [6], modified to 
include the effect of an applied longitudinal 
magnetic field. 

The assumptions (1) and (5) appear to contradict 
assumption (10) as well as the concept of appreciable 
heat loss, since the diffusion of a magnetic field into 
a conducting fluid proceeds at a rate proportional 
to that of the diffusion of heat. Thus it is necessary 
to show that there exists an interesting range of 
temperature 7’, pinch radius a and length L of the 
discharge in which 


(6 


(7 


— 


(1) 


where d¢ is the electrical skin-depth and 6, is the 
heat skin-depth. Generally 


de 1 


» 
«* (2) 


* “Now at A.W.R.E., Aldermaston, Berks., England 
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(2) 


where 7 is the electrical resistivity and x= K/(o Cg) 
is the thermal diffusivity, 9 and Cg being the density 
and specific heat per gram of plasma, respectively. 


Anticipating the formulae of Eq. (4) and the physical . 


constant of Eq. (32) one obtains 


=1.91n (3) 
where » is the particle density. As an example con- 
sider 7'=10%, n=10'% cm-? and Ind~l0 then 
de/dn ~0.002. Thus for a discharge of dimensions 
say, a=1 em and L=100 em the assumptions (1), 
(5) and (10) are compatible. Of course, in the vicinity 
of the electrodes these assumptions cannot be expected 
to hold, but then it is extremely unlikely that any 
of the assumptions (1) to (11) will be valid in these 
regions. 

These assumptions differ from those of KAUFMANN 
and Furrs [3] by the inclusion of radiation losses 
and the effect of a longitudinal magnetic field. In a 
later section the effect of including terms in the 
conduction equations relating to thermoelectric effects 
is investigated. 

While the assumptions are very restrictive and 
probably mean that the results cannot be applied 
directly to any known experiment, it is hoped that 
the results may be of value in estimating the impor- 
tance of axial conduction effects in a pinched discharge. 

It is assumed that the electrodes are maintained 
at T,—0° K, but it can be shown that electrode 
temperatures up to 1300° K do not appreciably 
affect the results obtained. Solution of the conduction 
equation enables an analysis of the stabilised linear 
discharge to be made in terms of seven parameters, 
namely, the electron density in the unpinched 
discharge, m9, the radius of the discharge tube, R, 
the length of the tube, 2Z, the initial axial magnetic 
field, By, two skin depth parameters /, uw, and the 
potential fall, 2q@, between the electrodes. 

In order to determine the time required for the 
establishment of a thermal steady state in a con- 
duction-cooled linear discharge, a rough numerical 
solution of the non-linear time dependent conduction 
equation was obtained.* A pulse of 10®° K was put 
on at the centre of the tube and maintained, whilst 
the electrodes were kept at 1300° K. For a tube of 
length 300 em containing a fully ionised gas of 
10! particles per unit volume the time required to 
achieve the thermal steady state is of the order one 
microsecond, that is some 4 or 5 times greater than 
the relaxation time to be expected from evaluation of 


L — 
V2 kT max/me 
To establish pressure equilibrium in deuterium gas 
requires a time of about one order of magnitude 
greater than that for the thermal steady state to be 
set up. Thus it appears that for a slow linear pinch, 


* This numerical solution was obtained by M. R. WetH- 
ERFIELD. 
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when loss processes other than bremsstrahlung radia- 
tion and axial conduction are insignificant, the 
steady-state problem is of genuine interest. 

In Section 2 the equations appropriate to the steady 
state are written down and the conduction equation 
solved. The results obtained are discussed in Section 3. 
The significance of thermoelectric effects is considered 
in Section 4. Finally, in Section 5 there is a general 
discussion of the results including a consideration of 
the regime in which they may be applicable. 

Electromagnetic units are used throughout except 
for the electronic charge e which is expressed in 
electrostatic units, and except where otherwise indi- 
cated. 


2. The steady-state problem 


The formulae taken for electrical resistivity 7, 
thermal conductivity A, and the amount of brems- 
strahlung radiation ¢, emitted per unit time per unit 
volume of plasma, are respectively, see SprTzer [7], 


° 
5 
K=KkK,T* (4) 
L 
T | 


with 
3 
Hjo= 6.53 x 10 In A emu (° K) * 


1.42 x 10-7 erg em* sec"! (°K) ?. 


The value for Ky, is given in Section 3. Since In A 
varies only slowly with the particle density (measured 
along the axis of the tube) n, and the absolute tem- 
perature 7’, this dependence is neglected. 

In the presence of a magnetic field both the thermal 
and electrical conductivities perpendicular to the field 
are reduced. It is possible therefore that the assumption 
of uniform temperature across the discharge might 
need refinement in order to deal with strong fields. 
However, for convenience we assume the uniform 
model to be appropriate both in the presence and 
absence of a magnetic field. 

The axial electric field EH, is related to the electro- 
static potential g by the familiar form, 


(5) 
x being the distance along the discharge measured 


from the anode. Other equations connecting the 
various axial electrical quantities are 


E.=jz (6) 
and 

djz 

dx 


Since the axial magnetic field B, is assumed 
completely trapped, it can be determined from 


By (x R?/ A), (8) 


where A is the cross-sectional area of the pinch. 


161 


ree 
kin 
red 
fa | 
cal 
ect 
ely 
che 
‘ed 
ion 
of 
Id, 
to 
me 
ies 
ge. 
ial 
gh 
ure 
rm 
er- 
he 
al. 
te. 
he 
to 
ial 
ict 
ito 
ial 
ry 
of | 
he 
(1) 
|| 
an 


A. H. DE BORDE, F. A. HAAS: 


The pinch relation of BENNETT [6] modified for the 
inclusion of an internal axial magnetic field is 
BA R? 

where N is the number of electrons (or ions) per unit 
length of discharge. 

The discontinuity in the B- field at the surface of 
the pinch gives rise to an azimuthal current of density 
je flowing in a thin skin of width 6 cm, such that 


B 
(10) 
If a is the radius of the pinch and noting that for 
d<a 


2nad 


(11) 
then the heat conduction equation to be solved is 


da =< joy. (12) 
Using the boundary conditions m—0, T= Tm when 
dT/dx—0, Eq. (12) may be integrated to give 


mKo m2 
where 
AI? B,? \2 t B,? R* 
(14) 


and J», a critical current analogous to that arising 
in the work of PEasE [8], is given by 


1 


( * 1.44 108 amp. (15) 


Inapractical case 7’, , the temperature of an electrode, 
is of order 1000° K while Eq. (13) would not be 
expected to hold for 7m<105° K, as the level of 
ionisation would be inadequate. Thus when the 
electrodes are maintained at+q,. a term of order 
(7'5/Tm)? can be neglected. It can be shown that 
allowing 7',--0 does not lead to any mathematical 
difficulty on account of the singularity in n. 

Using the boundary conditions d7'/dz =0 at T = Tm, 
T=0 at x=0, Eq. (12) can be made to give 


sin~ T/Ty, 
5 
= fsin® (16) 
0 
where 

Tm * jz——., (17) 

Fo 


thus exhibiting 7/7 as a universal function of «/L, 
that is 


F(;). (18) 


F (x/L) has been plotted by KaurMan and Furrtu [3]. 
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Since 7'=0 when x=2L it is simple to show that 


3 


=0.719= "Tm 


1 
_ 
2 (9/4) Pol L 


C.= 


Equation (19) shows the current density in terms of 
the average potential gradient and maximum temper- 
ature. C, appears as a correction factor to the current 
density to be expected from the supposition of a 
uniform temperature throughout the tube. 

Since the pressure throughout the length of the 
discharge is assumed to be constant, 


NT=NaT mn. (20) 


Nm, the number of particles per unit length at the 
centre of the tube, can be expressed in terms of the 
initial particle density, since assuming complete 
sweep-up of gas, the average number of particles 
per em. length of discharge before constriction must 
be equal to that during constriction. 

Using Eqs. (16), (18), (20), 


where 
1 
dy (9/4) 
C, = = ~1.22. (22 


C, is the factor by which the line density of the par- 
ticles is reduced from its original value at the point 
of maximum temperature. 

Using Eqs. (19) and (21), with the pinch relation, 
Eq. (9), for 7’=T7Tm, leads to the cubic equation 


4 
5 
Z (3°) Shae, (23) 
where 
1 
2 472 3 
Z=A Be (24) 
and the critical field B,., is given by 
3 1 
4L ng 2 on 


As an example for ~R?n,=10'* particles/em, 
4kV, L=225 ecm, R=20 em, A=0.3 then 
B.~4.4 gauss, so that in most cases of interest it is 
likely that By> Bo. 

After some elementary algebra one can derive the 
following result : — 


1 
ckT —8, + (8,2 + 48,) ? 
m = 1 ( 2) G (8;) (26) 
where 
2 
and 


1\ A 9 
= + ry 0.618 (I + ry (28) 
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Further algebra leads to the following results for 
and A, respectively, 


2 a 
_ ZR*C, By* Be? 
A= l6ngk Tm (30) 


If Q. is the heat dissipated through conduction 
via the electrodes and Qp that lost by bremsstrahlung, 
then it can be shown that 


I; 


This result indicates that for a sufficiently large 
axial magnetic field the radiation loss will be negligible 
compared with the thermal conduction loss except 
when unusually high currents flow. This effect is due 
to the “blowing out” of the discharge by the axial 
field trapped in the plasma. Numerical examples 
illustrating the above results are given in the following 
Section. 


(31) 


3. Examination of results 


According to SprrzEr [7] when no current flows the 
effective thermal conductivity is reduced due to 
thermoelectric effects. So far we have ignored these 
effects, but we shall see in Section 4, where a simplified 
model incorporating these effects is treated, that 
the reduced conductivity is appropriate in the present 
model. The appropriate value for K, is thus 


K, = 1.84 x 10-5 (In erg see! (°K)? . (32) 


In all calculations a mean value for In A equal to 
15 is adopted. 

We have seen in Section 2 that a typical value 
for B, is of the order 4.4 gauss so that usually for a 
stabilised discharge B,> B,. In these circumstances 
the Eq. (23) for Z can be solved approximately to 
give 


4 
1 ‘ 
(33) 
Eliminating Z between Eqs. (27) and (33) we obtain 
the following expression for s,: — 


4 3 
42 R* ny ck* i? 
B, 


Cy ng Ky 


This may be written 
8, ~ 0.783 ng) A? Vx (35) 
0 
where V is the voltage across the tube. 
The maximum temperature 7’, is given by Eq. (26). 
This may be simplified if s,<4/s,, that is when the 
discharge is primarily cooled by conduction, 


kT m= *) (vs.— +) (36) 
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The first term of the bracket is that which would 
arise if thermal conduction were the only heat loss; 
the second is a correction due to radiation losses. 
In the examples which follow we take A~y and 
since B.< By, by Eq. (33), we may take Z~1. As a 
particular example we may consider the following 
conditions : — 


Ng = 108 em-3, gy=1kV, A=—0.3, By= 1000 gauss 


(37) 
L=180cem R=45em. 
The dimensions L and R, for the sake of definiteness, 
are those occurring in “Columbus T2” (see BAKER 
et al [9]. By Eq. (25) this gives B.~6.3 gauss; also 
Eq. (35) gives s,~1.910-* The maximum tem- 
perature 7, can be determined from Eq. (36) and is 


Tm ~ 7.1x 10°°K. (38) 


The current J, is of the order 4.2 » 10° amp and 4, 
the cross-sectional area of the pinch, is of the order 
18 sq. ems. 

Applying Eq. (31) to the above example we find 
that Qp/Q.~1.2 « 10~* indicating that radiation losses 
are negligible. 

Consider, as second example, the following condi- 
tions : — 


Ng=10%em-3, gy= 10kV, A=0.3, By= 125 gauss (39 
39) 

L=180cem, R=45em. 

Under these conditions B.~11.1 gauss and s,~.003, 


thus 
Tm~71«10°°K, 10% amp, 


the last following from Eq. (29). The results arising 
for no axial field present are also obtainable from the 
previous expressions. The general expression for s, 
becomes 


3 
4x ck? 
Cy € Ko 


= 0.783 ny) V- 10-2. 
(40) 


&, = 


From Eq. (36) we see that the presence of the mag- 
netic field diminishes the effect of radiation loss. 
(8, <+/8,). Considering the conditions quoted in the 
second example above but with no magnetic field, it 
is found that 


Tm ~ 4.7x107°K. (41) 


We see that the effect of the magnetic field is to 
produce a negligible rise in temperature. 

Using the same conditions, but with no magnetic 
field applied, J, is found to be of the order 1.2 » 10% amp. 
Thus we see that the presence of the magnetic field leads 
to an increase in J, by a factor of order 6. 

The second example considered was selected in 
order to demonstrate the conditions under which 
radiation losses begin to become significant and, in 
fact, the value x R? n,~10'* is somewhat high 
for a practical discharge. A value 10'*® cm! might be 
more realistic, but in these circumstances radiation 


163 


of 
ent 
fa 
. || 
20) 
the 
the 
ete 
sles 
ust | 
21) 
22) 
ar- 
int 
on, 
ion 
24) 
27) 
| |_| 


A. H. DE BORDE, F. A. HAAS: 


losses are insignificant and the term s, in Eq. (36) 
may be neglected. 

Equation (26) may also be simplified if the dis- 
charge is primarily cooled by radiation, that is, if 
8,°>s,. In this case 


2 
G ~ (42) 
and 
kT m= 43 
m= (e/¢) — (43) 


The current J, is given, using Eq. (29), by the ex- 
pression 


As before, when B,> B.,thenZ~1 and when B,> By, 
then Z4( By/ B.) ~1. 


4. Thermoelectric effects 


So far the existence of thermoelectric effects has 
been ignored in consideration of the discharge. These 
effects produce a heat flow in the opposite direction 
to an applied electric field, and an additional current 
in the direction of a temperature gradient, the 
equations controlling the heat flow and current being 
given by Sprrzer [7] as 


P dT 

Q= —pE—K'-_. (46) 


In these equations, it is assumed that the heat flow Q, 
current density j, electric field HZ, and temperature 
gradient are all in the axial direction. « and # are 
given by Spitzer and Harm [10]. In the skin of the 
discharge Eqs. (45) and (46) hold without modification. 
In the main body of the discharge as before, we 
assume that no current flows due to inductive effects. 
Ignoring dependence on 7' through In A, « and f 
have the temperature dependence 


3 

a=a,7'* 
(47) 

B=Bo 
K'=K,'T?. (48) 


K’ occurring in Eq. (46) is, of course, the thermal 
conductivity unmodified for thermoelectric effects. 
Making use of Eqs. (45) to (48), we have in the skin 


5 
dQ 247 d 
de {Ko iz}— ‘dx (Bo (49) 
where 
Ky= Ky (1— (50) 
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Whilst for the main body, we have 


dx da 


To proceed further we interpret Q as an average 
heat flux by adding Eqs. (49) and (51) with weights 
A and 1—A respectively. We equate dQ/dx to the 
joule heating in the plasma, for simplicity neglecting 
any losses due to bremsstrahlung radiation. 

We thus obtain 


(Ko da ABo No dz +AEzjz4 0. 


(52) 


In obtaining Eq. (52) we have also used the con- 
tinuity equation djz/dz =—0. Substituting for using 
Eq. (45), Eq. (52) reduces to 


3 3 


da 
with 
= jz (Py — %)/Ko 
_ A 1 (54) 
The substitution in Eq. (53) 
dT RR 
dx (55) 
leads to 
d 
(PT) +ap+b=0, (56) 
which can be integrated in the form 
2 
In T+ 3 in(p +ap-+ b) D (57) 


where D is a constant. 

It can be proved using the results of Sprrzer and 
Haro [10] that 46—a?> 0 for all values of 2 between 
0 and 1.* 

Hence 


1 
q= + (4b—a?)? (58) 
is real; and Eq. (57) may be integrated to give 


; In (p? + ap + b) — (a/q) tan “P+ o+mnT=D. 
(59) 
The constant D may be fixed from the boundary 


conditions 7’=7,, when p=T?dT/dx = 0, giving 


D=In T—(a/q) tan— (a/q)+ nb. (60) 


* Comparison of the present work with similar but in- 
dependent work of M. G. Hatngs [11] has revealed that 
use of the transport coefficients derived by W. MARSHALL 
[12] leads to 4b—a* <0, and this implies an inward heat 
flux at the cathode. Thus there appears to be a discrep- 
ancy between the Marshall and the Spitzer and Harm 
transport coefficients. As yet, the source of this dis- 
crepancy is unexplained. 
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Tm can be determined by considering the behaviour 
of Eq. (59) in the vicinity of the electrodes where 
7-0, obtaining the total heat flux outwards from 
the tube and equating this to the electrical power 
dissipated. Since at the electrodes, the heat flow 


3 
Q=-—K,T*dT/dx remains finite, p=T*dT/dx 
must tend to + >.0 at one electrode and —oco at the 
other. At the first electrode we have 


a 


In pT = D+ (61) 
and at the second electrode 


In (— pT) = (62) 


The total heat dissipated through conduction is thus 


9K Decah(? 4 1 
2 Kye cosh 3) = + (63) 


which simplifies to 


— 
cos hr Vee (64) 
where r = a/q. 

Thus the results agree with those for the situation 
in which thermoelectric effects are neglected with the 
exception of the factor [exp (r tan~'r)]|/cosh (x r/2) 
which tends to 2/e for r -oco and to | as r +0. 
Substitution of a value 0.3 for A, and values of a». 
By from Spitzer and Harm [10] leads to a value of 
r~0.41 when [exp (r tan~!r)]/cosh (x 7/2) becomes 0.97 
which is negligibly different from its limiting value 1.0. 


5. Diseussion 


Before making any statement of the general con- 
clusions reached it is important to consider some of 
the underlying restrictions. 

First, the model only applies fully when equi- 
librium is attained between ions and electrons. Before 
this state is reached, however, the model will still 
apply approximately so long as we consider the 
temperatures involved as electron temperatures. 

Second, the only loss mechanisms we have con- 
sidered are bremsstrahlung radiation and axial con- 
duction. Many other processes are in fact possible 
and may be dominant under suitable conditions, 
see, e.g. Tuck [13]. 

Third, it must be remembered that all previous 
deductions have been based upon macroscopic 
equations which are dependent for their validity upon 
the condition 


dx 
|X dz < a (65) 


where X is any macroscopic variable and / the mean 
free path. For the macroscopic variable X = 7’ it is 
a simple matter to show that the above condition is 
equivalent to the fact that the rate of transport of 
heat energy must be very much less than the velocity 
of sound in the medium. 


LOSSES IN A STABILIZED PINCH 


Also this condition, Eq. (65), is a necessary condition 
for an approximately Maxwellian distribution of 
velocities to be maintained locally; it should be re- 
membered however that the well-known “Langmuir 
-aradox”’ deals precisely with the situation that a 
Maxwellian distribution of electron velocities is found 
experimentally to exist in situations where Eq. (65) 
fails to hold, see Gapor, Ash and Dracorr [14]. 
Thus the condition may be too stringent. 

In general, the condition of Eq. (65) cannot be 
satisfied throughout the entire length of the tube. 
Thus for temperatures sufficiently small, i.e. 7 <7”, 
the above condition will be violated. Our previous 
results have been derived assuming 7’, > 7'"~ 10° °K, 
where 7'” is that temperature at which a high degree 
of ionisation is possible. 

If 7’<T” the violation of Eq. (65) will only occur 
in regions near the electrodes where ionisation is 
incomplete and the uncertainty introduced is no 
greater than through the neglect of sheath effects. 
Even if this condition is not satisfied, providing that 
Tm> T’, then Eq. (65) will be satisfied over the 
major part of the length of the tube and again it 
might be anticipated that the formulae will hold 
unless conditions are dominated by the neglected 
effects in the neighbourhood of the electrodes. 

As regards the importance of thermoelectric pheno- 
mena, when the discharge is principally conduction 
cooled (as is likely in stable linear pinched discharges 
of current interest) the chief thermoelectric effect is 
to replace the unmodified thermal conductivity A’ 
by the modified A, and to slightly reduce the tempera- 
ture reached for a given voltage across the tube. In 
addition the original symmetrical distribution of 
temperature and density in the tube is replaced by 
an asymmetrical one. The precise form of the distri- 
bution and the modification to the formulae for 
current and discharge cross-section would require an 
integration of Eq. (59) but it seems unlikely that 
previous results will be changed drastically. 

Within these limitations, we may draw certain 
conclusions. First, in nearly all circumstances the 
main loss from a linear pinched discharge is likely 
to be from thermal conduction. In order that this 
should be so the following condition must be satisfied : 


3 


1 ( Be\? \2 
a ( "No K, ) (66) 


In these circumstances the highest temperature that 
can be reached, taking the thermoelectric effect to 
be negligible, is 


Tm =4.6V VA(1 + Aju) 108° K, (67) 
V being the voltage across the tube. 
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LONGITUDINAL OSCILLATIONS IN A NEUTRALIZED ELECTRON BEAM WITH A 


BOUNDARY 


M. YosHIKAWA 


DEPARTMENT OF Puysics, UNIversiry oF Tokyo, Tokyo, JAPAN 


Longitudinal oscillation of a neutralized electron beam with a boundary in a strong axial magnetic 
field is discussed. The beam is assumed to be cylindrical and to be confined in a conducting cylinder, 
or to be planar and to be held between two conducting plates. A sufficient and, in some cases, 
approximately necessary condition for stability is obtained. To be stable, the beam current should 
be below a certain value which depends on the electron velocity, the ratio of electron to ion mass 
and the geometrical dimensions of the beam and the conductor. The validity of the approximation 


is also studied. 


1. Introduction 


In this paper we shall consider longitudinal oscilla- 
tions in cylindrical and planar electron streams 
neutralized by positive ions in the absence of collisions. 
Many authors have discussed similar problems for the 
case of infinite streams, taking into account the 
effect of thermal motions in the Boltzmann equations 
{1.2.3}. Their results show that, unless the energy of 
the thermal motions is comparable to that of the 
ordered motions, the system is unstable with respect 
to oscillation. This self-excited oscillation may prevent 
production of electron beams of high intensity although 
the oscillation may also be advantageous from the 
thermonuclear point of view in that oscillations of 
ions are stimulated and their random energy raised. 

The case of bounded streams may be of particular 
interest as the situation is really true in plasma 
betatrons [4], high intensity accelerators of neutra- 
lized electron beams, in the some thermonuclear 
machines. 

For tractability we shall hereafter neglect the 
random motions of electrons and ions and _ shall 
assume that the beam holds a strong axial magnetic 
field within it and that the magnetic field produced 
by the beam itself is negligible. In what follows, a 
stability criterion will be derived and discussed. The 
criterion states that the longitudinal oscillation can 
not be excited spontaneously if the beam current is 
below some critical value which depends on experi- 
mental conditions. The validity of the above assump- 
tions will also be discussed. 


2. Dispersion relations for the cylindrical case 


We consider a uniform cylindrical beam of neutra- 
lized electrons. The beam has a radius d, is streaming 
along a strong axial magnetic field and is enclosed in 
a conducting cylinder of radius D. The intervening 
space (d<r<D) is assumed to be vacuum. Using 
cylindrical coordinates coaxial with the cylindrical 
beam, we can find a stationary state of the neutralized 
beam. This state is characterized by the velocity 
vo (r, 6, z components=0, U9), the velocity 


of ions Vio (0, vig, 0), the partial pressure and the 
density of electrons peg and ny, the partial pressure 
and the density of ions pjg and ny, the electric field 
E, (Ew, 0,0) and the magnetic field B, (0, Byo, 0). 

In order to investigate the behaviour of the beam 
for small deviations from the stationary state we 
expand these quantities to the first order. The first 
order approximation for the electron velocity, the 
ion velocity, the density of electrons, the density 
of ions, the electric field and the magnetic field are 
to be denoted by 


Ve(Ver,Veo.Vez),Vi Ne.m,E( -),B( B,, Bo, Bz) 


respectively. Some of the basic equations are closed 
in themselves, which are given (in MKS units) by 


m Up Ver =—eE 1 
eV (5, + | ez ez ( ) 
(2) 
One Ne O Vez 
+ U9- No == 3 
ony CVjiz 
; = 0 4 
( ) 
(5) 
ozor 02? 
l GE, E, 
r or 
, 1 OR, diz 6) 
—€Ng Vez — ENe Vg + CNg Yiz —tz (7) 


where we have assumed that the deviations are in- 

dependent of 6. As for the case when they vary as 

e'™® some discussions are given in Section 4. Here 

iz is the z-component of the current of the first order 
1 


and y means, as usual, (1—v,2/c?) *. wand e represent 
permeability and dielectric constant of vacuum. 

In the right hand sides of Eqs. (1) and (2) the 
terms —eve, Boo + even B, and + evjr Boo, respec- 
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tively, have been omitted in assuming that the 
magnetic self-field of the beam is negligible. This 
means that the pressure peo + pio is also negligibly 
small. Consistently, we shall assume that the energy 
of the random motions of electrons and ions is much 
smaller than the energy of the unidirectional motion 
of electrons. Also in the left hand sides of Eqs. (3) 
and (4) we have dropped the terms ny @ve,/@r and 
Ny Ovj,/Cr, respectively, keeping in mind that the 
axial magnetic field constrains the ions and the 
electrons to move along the lines of force. 

Now we study a normal mode that varies with ¢ 
and z as exp (iwt—ikz). Substituting this into Eqs. (1) 
to (7) and rearranging, we obtain an equation for 
the amplitude of the axial component of the electric 
field, 


1 ds dE 
(8) 
where 
(w — kv,)? I, (9) 
k’2— k®?—(/c)?, (10) 


In vacuum (d<r< D) it follows 
1 a dE; 


‘ 
r dr dr =0. (12) 

The solution of Eq. (8) is J, (k’sr) and the solutions 
of Eq. (12) are linear combinations of J, (k’r) and 
K, (k’r). Furthermore, boundary conditions are 


E.z(D)=0, (13) 

(d+0) — Ez (d—0)=0, (14) 
dEz(d4+0)  dEz(d@—0) 


The first and the second conditions imply that parallel 
components of the electric field are continuous at 
the boundary ; these can be generally derived from one 
of the Maxwell equations. The third condition is to 
be deduced from the consideration that in case of 
pure longitudinal oscillations there appears no surface 
current at the boundary (r=d) of the beam. 
Finally, substituting the solutions of Eqs. (8) and 
(12) into Eqs. (13), (14) and (15), we obtain a dispersion 
relation for the normal mode. Introducing new para- 
meters: 
J=¥D, (16) 


the dispersion relation becomes: 


J, (8) Ko (X) (x) + Lo (X) Ky _ D 

Ty(X) Ko (a) — Ky (X) (2) =¥( d 4D). 
(17) 

Roots of Eq. (17) are many-valued and we denote 

them as S (X/x, X;j), where the magnitude of the 

roots increases as j= 1, 2, ... and so on, corresponding 

to the number of loops of Z, in the radial direction. 
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Fig.1 S(D/d, kmin D; 1) versus kmin D for the 
cylindrical case. This characteristic parameter should 
be used to obtain stability conditions in Eq. (21). 


S (Did, kD;1) with k= kmin (see discussion of 
Eqs. (20) and (21)) is shown in Fig. 1 for the cylindrical 
case and in Fig. 5 for the planar case (which we 
shall treat in Section 3). The curves are shown as a 
function of kD, regarding D/d as a parameter. The 
values chosen for D/d are 1, 2, 10 and 100. kD can 
range from 0.1 to over 20. In the cylindrical case 
S (Did, kD;1) is a slowly varying function of kD. 
It can be seen from Fig. 1 that the largest value of 
d(log S)/d(log kD) is about 0.2. The dependence be- 
comes stronger in the case of the planar beam, see 
discussion in Section 3. Nevertheless the maximum 
value of d (log (log kD) for a wide range of 
parameters is 0.5 or less. Therefore it follows from 
Eq. (16) that, for a fixed value of D/d, the quantity s, 
which is S (D/d, k’ D;1)/kd, is a monotonously de- 
creasing function of kd because of the predominating 
dependence on the kd in the denominator. 

The dispersion relation composed of Eqs. (9) and 
(17) can be solved by graphical manipulation as 
follows. We write down the dispersion relation in 
the form of 

Ope? Opi? S? (Did, k’ D; j) 
— oe i+ k’? (38) 


(j=1, 2....) 


We shall plot the both sides of the equation as func- 
tions of for fixed values of D/k, d and D. We have 
shown in Figs. 2 and 3 two typical figures; the former 
has four real roots (stable case) and the latter has two 
real roots and two (mutually conjugate) complex 
roots (unstable case). The critical case, when we have 
two single real roots and one double (real) root, 
defines a boundary that separates the stable region 
from the unstable region in parameter space. This 
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0 ke 


oe, 


Fig. 2. Graphical method of solution is indicated 
schematically. The solutions are defined as the inter- 
section points of the curves corresponding to both sides 
of the dispersion relation, Eq. (18). Curve A is a plot 
of the left hand side of Eq. (18), curve B the right 
hand side and curve C the right hand side with k’ re- 
placed by k [see Eq. (18a)]. In this figure there are four 
real roots (stable case). 


can be obtained graphically without much trouble 
using Figs. 1 and 5 to obtain S (D/d, k’D;1). Here we 
have taken j= 1 because waves with higher values of 
j are more stable; this is so because in general 


S (Did, k’ D; j>1)>S (Did, k’ D;1) holds. 


0 = kc 


Fig. 3. Graphical method for an unstable case. The 
curves are labeled as described in the caption for Figure 2. 
There are two real roots as intersection points. In addition, 
there are two complex conjugate roots. 


However, it would be quite convenient if we could 
have a simpler solution of the dispersion relation. 
Now we can derive a sufficient condition for the 
stability. As shown in Figs. 2 and 3, when k’ in the 
right hand side of Eq. (18) is replaced by k, this 
side corresponds to a straight line parallel to the 
axis and is constantly below the original function 
1 +[S? (D/d, k’ D;1)/k’ *d?|. Therefore if the equation 


Ope ~3 Opi? 
(w— kv,)? 


S?(Djid, kD; 1) 


k? d? (18a) 


has four real roots, the original Eq. (18) also has four 
real roots. The sufficient condition can be derived if 
we note that the minimum of the left hand side as a 
function of @ is equal to 
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and 
2(Did, kD; 1 
S (Did, ) Ope * + Opi 
k2 d2 = k? (19) 
for all k. 
In this connection, when mey* < mi, 
a 
k.2e2 2 212 


is a small quantity and the former stability condition 
can be regarded as an approximate necessary con- 
dition, where «, and k, are the values associated with 
the equality sign in Eq. (19). 

D and d are parameters that characterize an 
experimental device. So is kmin the minimum value 
of k, if we apply a periodic boundary condition in z. 
It turns out by numerical calculation that for fixed 
values of D and d, y(D/d, kD) and consequently 
S (Did, kd;j) decrease monotonously as k decreases. 
Therefore the sufficient condition for the stability of 
the j-th mode becomes 


2 272 
S? (Did, kmin D; }) Lape wpi 


1 + . (20) 
Inequality (20) is satisfied when and only when 
inequality (20) for j7=1 is satisfied. Introducing a 
dimensionless parameter, v = d*)(ny where r, is 
the classical radius of electron e?/4 ~ eme, Eq. (17) for 
j=1 becomes 


l 9 
[S? (D/d, kmin D; 1) + kmin? d?] 


S (Did, kmin D:1) as a function of kpinD for typical 
values of D/d is shown in Fig. 1. When kypinD ap- 
proaches zero, S(D/d, kminD;1) tends to a finite 
value S (D/d, 0;1), as shown in Fig. 4. Furthermore 
S (D/d, 0;1) has the following limiting properties, 


S(D/d, 0; 1) +2.405[1—log (D/d)], {log(D/d) (22) 
1 


2 2 
S(Djd, 0; 1) > | {log (D/d) +00}. (23) 


Hence, the sufficient condition for stability for 
min D? <1 is 
1 


my 


log (D/d) <1, 


169 


| 1+ 
A 
| | 
| 
| 
| 
i 
| 
A | 
c 
y 8, 
de- 
A 
A 
as 
in 
7 
ne- 
ave 
ner 
lex 
ave 
ion 
his 
|_| 


on 


M. YOSHIKAWA: 


0.5 flog (1 + [mest * (25) 


log (D/d) >1. 


$(0/4,0;1) 


Fig. 4 S(D/d, 0;1) versus D/d for the cylindrical 
ease. This characteristic parameter should be used to 
obtain stability conditions in Eq. (21). 


3. Dispersion relation for the planar case 


For a planar neutralized electron beam of thick- 
ness 2d between two parallel conducting plates 
separated by 2D, where D>d, a dispersion relation 
similar to that for the cylindrical case can be obtained : 


Stan S = x cot (X— 2). (26) 


We can proceed as before and have only to present 
S (Did, kminD;1) and S (D/d, 0;1) in Figs. 5 and 6 


& 


$(D/4,K min O51) 


01 


1 0 100 
K min 
Fig. 5 S (D/d, kmin D; 1) versus kyjin D for the planar 


case. This characteristic parameter should be used to 
obtain stability conditions in Eq. (21). 
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Fig. 6 S(D/d, 0; 1) versus D/d for the planar case. 
This characteristic parameter should be used to obtain 
stability conditions in Eq. (21). 


respectively. Eq. (18) holds good in this case and 
instead of Eqs. (19) and (20) we have, for kj,inD? < 1: 
1 


0.62 {1—2[ (D/d)—1]}} {1 +( 


mj 


(Did) —1€1, 


3 + 
0.25 (D/d) (1—y-*) + mer >» (28) 
(D/d) >1. 


4. Discussion 


In Section 2 we assumed that the field E is uniform 
in the azimuthal direction. We can proceed as before 
when it varies as exp (im). Also we may assume 
that the field changes in the direction parallel to the 
conductor plane and normal to the direction of the 
static magnetic field. Then similar dispersion relations 
can be obtained, but they are, as numerical calcula- 
tions show, less restrictive conditions than Eq. (21). 

The stability criterion Eq. (21), is not a very severe 
restriction on a device of usual scale when the beam is 
relativistic. As [1 + is an increasing 
function of y with a value between [1 + (me/mj)'/*}-* 
and (m,/m;)-!, an electron current of hundreds of 
thousands amperes can flow without causing the 
instability. When the beam is non-relativistic, the 
effect of the conducting wall can no longer stabilize 
the beam and y decreases as S??/4. 

In Section 2 we have neglected —e (ver Boo —veoo Br) 
in the right hand side of Eq. (1), evj, Bay in the 
right hand side of Eq. (2), m9 @ver/@r in the left side 
of Eq. (3), and my @vj,/@r in the left side of Eq. (4). 
These simplifications, which are correct when the 
axial magnetic field is infinitely strong, have been 
convenient for us in extracting pure longitudinal 
oscillations and in studying the two-stream instability 
in a neutralized electron beam with a boundary. 
Without these assumptions the longitudinal oscilla- 
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tions would naturally couple with transverse oscilla- 
tions and the resulting dispersion relation would be 
much more complicated. In this connection, the 
conclusion reached in Section 2 and Section 3 should 
be altered for frequencies in the vicinity of the fre- 
quencies of the ion and electron cyclotron waves. 

It is difficult to express in a closed form how far 
the simplifications can be justified. Accordingly, we 
take two typical cases that can be treated analytically. 
The first case is when k-0 and the second is when 
k=ke, where k, is defined by Eq. (19). The former 
is significant for the reason that the two-stream 
amplification is more dangerous in the long-wave 
limit. Using new parameters, 


P=celMpe (29) 
q=(d/c) (30) 
Me| mj (31) 
b= Boo/ Bzo (32) 
Cee = (e/me) Bee (33) 


and further assuming for simplicity 


S2 


eae >! (34) 
d? wy 
St? <1, (36) 
we obtain 
7 
virlOr | | pty? 82) ‘ 


|W | <! (38) 


ls Bp? 
Br | » B 4 


ws (1— | 
Ver Boo | | bp @ve,/@r | 
E, | | |<! (40) 


vir Boo bpqu Ovj,/er 


Ez BS | <! (41) 
7 
W y? y2q ps yi q? (42) 


These inequalities can be satisfied by an infinitely 
large p. It follows further that gq, and therefore d, 
cannot be made arbitrarily small and that in this 
limit p? should be much larger than Iq 2u?). 

Ina contrasting case when k? = k,?, that is, when the 
system is near the critical thershold of stability, 
we obtain, 


| B%y4S2 


(43) 
| | pty? S? | 
‘00 B, | | 28 
Ver Boo | | OVer/Cr | 
Ez || BS (46) 
+ pty), (47) 


p* 


1 4 
where we have assumed y>1 and 
3 


By? Siq<1, in order to yield a brief result. These 
relations are again satisfied by an infinitely large p 
and do not hold when q is arbitrarily small. 
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RADIATION FROM A MODULATED BEAM OF CHARGED PARTICLES 
PENETRATING A PLASMA IN A UNIFORM MAGNETIC FIELD 


E. CanosBio* 


Max-Pianck-Instirut FUR PHystk UND ASTROPHYSIK 


MincHen, FEDERAL REPUBLIC OF GERMANY 


The radiation from a density-modulated beam of ions, which penetrates a plasma perpendicular 
to astrong magnetic field B,, is investigated in two simplified situations: (a) the beam is an 
infinite plane parallel to B,. and (b) the beam is an infinite cylindrical surface parallel to B,, 
the radius of the cylinder being the gyroradius of the beam particles. This latter beam can be 
ideally constructed by injecting into a plasma a linear beam, modulated at a frequency which 
is an integral multiple of the gyrofrequency of the beam particles and incident in a direction 
which forms a very small angle with a plane perpendicular to B,. 

In both situations some resonances of the Poynting vector are found. The resonance, which 
occurs when the modulation frequency is equal to the “ion-resonance”’ frequency, is specifically 
investigated, taking into account the finite electric conductivity of the plasma. It is shown that. 
under appropriate conditions, the beam-plasma interaction at this resonance becomes very strong. 


1. Introduction 


The radiation from an intensity-modulated beam 
of ions directed into a plasma in the presence of a 
strong magnetic field, constant in space and time, 
has recently been investigated by KippENHAHN and 
DE Vries [1]. The investigation was made in con- 
nection with the problem of the injection and heating 
of a plasma. Kippenhahn and de Vries considered 
modulation frequencies sufficiently low, plane beams, 
and infinite electric conductivity of the plasma. They 
found that the energy loss from the beam is always 
negligible, except when the velocity of the beam 
particles is equal to the Alfvén velocity. 

In the present paper we provide an extension which 
is useful in all ranges of frequencies and takes account 
of the curvature experienced by a two-dimensional 
beam in an external magnetic field. We show that 
the resonance of the kind found in [1] is practically 
destroyed by the beam curvature, but that other 
infinities occur in the Poynting vector at particular 


x 


Fig. 1 The plane beam. 


* On leave of absence from the Research Department, 


Fig. 2 The curved beam. 


values of the modulation frequency. By taking into 
account the finite electric conductivity of the plasma, 
the peak value of the energy loss from the beam at 
the “‘ion-resonance”” frequency is determined and 
briefly discussed. Such a peak value has, in many 
practical situations, a great order of magnitude and 
is a maximum when @,? is of the order of @j @e. 
The beams here considered are: 
i) an infinite plane parallel to the strong magnetic 
field B,, constant in space and time (Fig. 1). 
ii) a cylindrical surface parallel to B,, whose 
radius is the gyroradius of the beam particles 
in B, (Fig. 2). 

We assume that all the ions in each beam travel 
perpendicularly to By with the same speed, that the 
velocity is not modified by the interaction with the 
plasma and that the plasma is infinitely extended. 


Euratom, Bruxelles, Belgium. 
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Two-dimensional beams are considered here, as in 
the paper cited above, because of the great simplicity 
attained when beam and plasma oscillations have 
the same simple symmetry*. 

The beam of the case ii) is ideally constructed by 
injecting into a plasma a linear beam, modulated at 
a frequency which is an integral multiple of the 
gyrofrequency of the beam particles, in a direction 
which forms a very small angle with a plane perpen- 
dicular to By. 

The plane beam radiates by Cerenkov effect, ie. 
only when the ion-velocity is greater than the phase 
velocity V=c/n of the waves in the plasma; on the 
contrary, the cylindrical beam can radiate also when 
this condition is not fulfilled, because of the curvature 
of the trajectories of the beam particles. 


2. Basie equations 


Consider now the following linear, nonhomogeneous 
system of equations [2], [3]: 


ccurl B = + 4x (j + jy). (1) 

ccurlE = —B, (2) 

x Bo, (3) 

+7) B+ By (m—m)¥, (4) 
div B= 0, (5) 

div E = 4x (e + é9) ; (6) 


where B and E are the magnetic and electric fields 
respectively, j and j, are the electric current density 
of the plasma and the beam; V, 09. @p are the velocity, 
unperturbed mass density and frequency of the plasma; 
y is the mean collision frequency between ions and 
electrons in the plasma and e is the charge of the pro- 
ton; m,. and m; are the electron and ion masses res- 
pectively. For simplicity, Eq. (3) does not include 
the pressure term. Eq. (6), where ¢, is the electric 
charge density of the beam, is only a definition of the 
electric charge density ¢ of the plasma and will not be 
used in the following. The solutions of Eqs. (1), (2), 
(3), (4) and (5) which describe the fields radiated from 
the beam are obtained from the general integral by 
assuming “‘Ausstrahlungsbedingungen”, i.e. discard- 
ing fields which propagate from infinity towards the 
beam. In the plane case we chose Cartesian, ortho- 
gonal coordinates (a2, y, z) such that: 

B= 

(7) 

j=6 [J, + 


where 6 (x) is Dirac’s function, the constants J = New 
and J,—N,ew are the modulated and unmodulated 


* The only linear beam which has simple symmetry 
propertiesi.e., a beam parallel to By, will be studied in a 
future paper. 
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surface current densities of the beam, respectively. 
and w =«/l is the velocity of the particles in the beam; 
N,>WN and N are positive real numbers; @, / are 
positive and real. 

In the cylindrical case we introduce cylindrical 
coordinates (r, 9, z) such that: 


= 
and (8) 
jo (r — R) (1, + @,; 


where m is a positive integer and m=m wg; wg 
=e B,/ Mc is the gyrofrequency of the beam particles, 
whose mass is ./; the radius of the beam is R = w/mg 


and, as in the plane case, J= New and IJ, = N,ew. 


3. The plane beam 


In the reference system introduced above we look 
for solutions of Eqs. (1), (2), (3), (4) and (5) which are 
independent of z. 

If B, is a time-independent perturbation of the mag- 
netic field, we obtain 


B,= [—e (a) I,+ const. | (9) 


where ¢ (x) is equal to 1 if x>0, and —1 if x <0. 
We set the constant equal to zero and suppose that 
(2x/c) 1, < By, so that the steady magnetic fields 
on the left and on the right of the beam are ~B, 
and we can neglect B, in what follows. 
Now we look for solutions of Eqs. (1), (2), (3) and (4) 
of the form 
B= B(z)eily— (10) 


and analogous expressions for E, j and vy. In this way 
Eq. (5) is identically satisfied and the system Eqs. (1), 
(2), (3) and (4) separates into two independent parts. 
The first one contains only the components By, By, 
E;, jz and v; and may be reduced to 


Ez" (x) + (Ky? —P) Ez (x) =0, (11) 


where Ky? = (mng/c)? and ng? = 1 + is 
the square of the Eccles refractive index. The second part 
determines Bz, Ex, Ey, jx, jy, vx and vy and is reduced to 


Bz" (x) + (K?—BP) B. (x)= —4* (x) 
(12) 


where A?=(mn/c)? is the square refractive index, 
found by List [4] for plane waves ina plasma with an 
external magnetic field: 


1+ Mp? + -i yor) 
and (13) 


© Op? (14) 


Eq. (11) describes perturbations independent of the 
beam. They are either waves propagating from 
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+ coto — oo(or reciprocally), if K,? is real and greater 
than /?, or monotonic unbounded functions, if Ky? is 
real and positive but smaller than /?, or a product of 
such kinds of solutions, if Ky? is complex; so we choose 
the trivial solution 2,=0. It then follows: 


B, = B, = jz = 02 = = 


On the contrary, when we impose w> V (a restriction 
analogous to Gerenkov’s condition) and the “Ausstrah- 
lungsbedingungen” on Eq. (12), it then supplies the 
magnetic field radiated from the beam: 


B: (x) =0(—2) bye ~iVK*— Fx + Q(x) bei VK*— (15) 


where 4 (x) =1 if « >0, and 0 if 2 <0. The constants 

b, and b, satisfy identically the equation: 

+2iVK®—P (b, + by) (x) =— (x) 
(16) 

Remembering that 

f (x) 6 (x) = —6 (x) f' (x), 
provided that 
{(0)=0, 


we get 
4x — 4x 
—b,+b,= ——*1, iVK*—# (6,+ 6,) = — Ila, 


and (17) 
9 « 


Bz (x)=-— -{-6(- 


=q{-0(-2) 


(x) [VK?-2 + ila]ei VK-F x}. (18) 


Then from Eqs. (1), (2), (3), (4) and (18) we obtain 
Eq. (19), see bottom of page. 


We are interested in the time-mean value of the 
energy flux, given by the real part of the complex 
Poynting vector [5] 


S= Re (Ex B*), (20) 


where B* is the complex conjugate of B. 

In the remainder of this Section we shall consider, 
for simplicity, a non-dissipative plasma i.e. we will 
put y=0, and we shall study Eq. (20) only in the 
frequency regions in which n,?, defined as 

(21) 
is >1 in the other regions; there cannot be energy 
propagation: this fact is obvious for ny? <0. If, on the 
contrary, 0<n <1, we still assume that the energy 
flux is equal to zero due to the fact that the waves are 

monochromatic. In this way we find: 


( I K? — I? (1 — 

— rol K*—P(1— 2?) 
0. 


As a Cerenkov radiation, the vector S, as all other 
quantities, is independent of the mass of the beam 
particles; but we have to keep in mind that in a magnet- 
ic field a beam, carrying a given charge and current, 
may be considered as a plane beam only when the 
masses of its particles are very large. 

The Cerenkov angle @ between the vectors § and w 
is given by 


(a)V K2/P—1 =e (x)Vw?/V2—1 =tan8. (23) 


Hence the energy propagates with the velocity c/n, 
in the two directions which form angles # and —#@ 
with w. 


[l(1-a?) VK?-P-ia Px 4 (x) K? + ia K2]ei Px} 
— [0 (— 2) 4 0 (a) ef 


{0 (ila —VK?— P) [1(1 — n,?) +iaVK?— P| 


+0 (a) (ila +VK?— P) [1 (1— ng?) —iaV K?—P] ei Px} ; 


— {0 (ila —V K? —P) [la — i(1 —n,?)V iVK*— Px 


+ O(a) (ila +VK?—P) [la+i(1—n2)VK? —P]ei =} ; 


2xIw { 
(cK)*/K?— 
2xIw 
=~ 2 (1— a? 
‘ 2n/K? — 
Copw 
_ id, . 
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In order to discuss other properties of Sx, it is con- 
venient to use the form: 
x Ne\2 wV (w?— V?2(1— a? 
= Vw? — V2 
and to express 7,” and « by means of the dimensionless 
parameters 


(24) 


c 


(@e — a)? 25) 
— 


We We Me Me 


The results are 


zx/u—(l—x-+ o)? 
"0 +6) (1—x) 
and 


It appears from Eq. (24) that in addition to the reso- 
nance for w= V, which is of the kind found by Kippen- 
hahn and de Vries, S, has infinities in the resonance 
frequencies of n?, (and of « consequently), resonances 
which have already been studied by K. KOrPEr [6]. 

When o<1, n,? has poles at x~u and x~1/y; 
when o> 1/u the poles are at x~1 and x~o. In 
Kérper’s terminology the resonances at x~ and 
z~l are “ion-resonances’; those at x~Il/u and 
x~o are “electron-resonances”’. The last case, x~o, 
will not be considered because in the immediate 
neighbourhood of o there are also the two zeros of 

Supposing now w> V, in the first three cases we 
have the following simplified expressions for Eq. (24): 


on 


a) when o <1 and ; 


1 


Sx~ €(x) 5 (28) 


— 


b) when o > 1/u and wa <1, 


y 1 
c) when o< 1 and > 1, 
1 
Sy ~ (x) e( .) (30) 


Note that all these expressions are independent of 

the velocity of the beam particles. 

Outside the resonances of n,? we consider two ex- 

treme regions: 

a) Very high frequencies, i.e. x > Max (oc; 1/u): here 
Sx has to be assumed equal to zero, because of the 
inequality ?~1—(a/x) <1. 

b) Very low frequencies, i.e. x< Min [1; u(1 +<o)}: 


here we get 
(et) | 


Ne we 
Vw — [c?/(o + 1)] 
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In the case og>1, i.e. for high density plasma, this 
expression has already been found in [1]; it gives an 
extremely small energy loss from the beam far from 
the resonance w=V. For instance, when x <1 it is 
a? <I1/u; then S, becomes: 
= (Ne\2 uw? V 

the mean relative energy loss of an ion per em of path 
6, E/E, given by 


N, = us 


(33) 


is, in the case N — N,, smaller than (x) rj x 
N (V/w), where rj = e?/Mc? is the classical ion radius. 

For a deuteron beam with V ~10* 6,#/E has 
an extremely small value: V/w 

Unfortunately, measurable radiation cannot be 
obtained even in the resonance of Eq. (24) at w= V, 
as we shall see in the next Section, because this kind 
of resonance is destroyed by the curvature experienced 
by the beam in By. Appreciable energy might therefore 
be emitted only in the resonance of n,?. We shall 
discuss this point in Section 5, including, partially at 
least, the effect of the finite electric conductivity of 
the plasma. 


4. The curved beam 


As in the plane case, we suppose (27/c) 7, < By and 
look for solutions of Eqs. (1), (2), (3) and (4) of the 
form 


B = B(r) ei(mo— ot) (34) 


and analogous expressions for E, j and v. Then we 
satisfy Eq. (5) identically and get two systems of 
equations. The first one determines B,, By, Ez, jz 
and vz; because it turns out that they are independent 
of the beam quantities, as in Section 3, we choose the 
trivial solutions B,= By = Ez=jz=vz=—0. The second 
system of equations furnishes the components B,, 
E,, Eo, jz. jo, Vr and ve and reduces to the Bessel 
nonhomogeneous equation 


d 


dr dr 


- {5 + (35) 


Bz(r)) + B.(r)= 


where K? =(@n/c)*; n? and « are defined by Eq. (13) 
and Eq. (14). 

In order to simplify the calculations and discussions 
we assume in this Section, as we have done for the 
plane case, that n?=n?,>1, i.e. that the plasma is 
nondissipative and that the energy can propagate. 

By using the *‘Ausstrahlungsbedingungen’”’ and regu- 
larity conditions in r=0, we get the solution of 
Eq. (35): 


B:(r)=0(R—r) b, J m(Kr) +6 (r—R) b, H™ (Kr), (36) 
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where Jm(Kr) and H (Kr) are, respectively, the 
Bessel and Hankel functions of the first kind. The 
constants b, and b, satisfy the equation: 


I m(Kr)] (r—R) 
+ (Kr), n(Kr)) 
+ — R) 
—4n . —ma 
= 8(r—R) +2 (37) 
If we proceed as in Section 3, we obtain 


66 (R—r) hm (a, K R) Im (Kr) 
+0 (r—R) jm (a, KR) H®(Kr)} (38) 


where 
hn (%, KR)=m a H(K R)+ KR HY’ (KR) 
jm (a, KR)=madm(KR)+KRIn' (KR). 

In the above equations the functions J»’(K R) and 
H\)’ (K R) denote the derivatives with respect to 
Kr of Jm(Kr) and H') (Kr), respectively, evaluated 
at r= R. The function 6 (x), defined previously, is 
unity for x >0 and zero for x <0. 

The other components follow from Eqs. (1), (2), 
(3), (4) and (38) and are given at the bottom of the 
page, in Eq. (40). 

The components of the real part of the complex 
Poynting vector (20) are: 


o 
S,= 0(r— R) [jm (a, K R) |? 


(Zz) {0(R—r) \hm(x KR)? 
x |" Im? (Kr) KI m (Kr) (Kr)| 
+ O(r—R) |jim (a, KR)? 
Re (Kr))? + KH® (Kr) (Kr) |} 
S:=0. (41) 


As we are interested in the radiation close to the 


beam we consider S, on the external side of the cylin- 
der surface r= R: 
(> w 
J m 


2ix*] 
nes ong 
2x2] 


Now we consider (20) in the neighbourhood of the 
poles of n,?, where w> V and: 


w 2m+1 
Im(m 4 cos (m z). (43) 


By writing n,? and « in terms of x, o and y, as in 
Eqs. (26), (27), we get the following approximate 
expressions : 


a) “Ion-resonance” region: <1, 


1 1 
nf FH 2 2 w (on 2 2m-+1 
(44) 


b) “Lon-resonance” region: >1/u,au(1—z)<1 


1 
l-x 4 


(45) 
c) “Electron-resonance” region : < > 1 


1 
S?7 (R) ~ (Ne) 
1 
cos [m (aac) 4 (46) 


These expressions, as functions of no, have the form 


S? 
oplo(1-x)] 


C1 Ng (Cy + , 


where c¢,, c, and c, are constants. When n,->0o, they 
have extremely rapid oscillations. Hence it has perhaps 
more physical meaning to replace them by the mean 
value in a period 27: 


+ 27 
Ng COS” (Cy pq’ Cy) A(egny’) (47) 
Ng 


In this way, if c,ny> 27, Eqs. (44)—(46) become: 
1 


S?(R) ~ (Ne)2e (44a) 


S>(R) ~ (Ne)2e 


y 
S?(R) ~ —). (46a) 


{0( Rr) K R) Im Kr) +0 (Kr)| + O(r—R) jm(a, KR) [+ H (Kr) + KHO 


KR) (Kr) + Km (Kr) +0(r—R) jm (x, KR) Kr) 4+ KH’ 


-irl 
—iB, 
0, = = 


(40) 
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In each of these cases, the radiation is just the sum 
of the radiation emitted from the two sides of the plane 
beam in the corresponding situations, Eqs. (28)— 
(30). 

With the exception of these cases, when R<oo, 
S? (R) has no poles and satisfies the inequality: 


wi 
SP (R) < 


which follows from Eq. (42) by remembering that [7] 


V2 


Moreover, it is easy to see directly from Eqs. (42) 
and (43) that, when o>I/u.w> V and x<1, i.e. at 
low frequencies and high density plasma*, S? (R) 
has the same extremely small order of magnitude as 
indicated in Eq. (32). At very high frequencies, i.e. 
“> max (0, 1/u), S? (R) is to be considered equal to 
zero because V >c (see last Section). 

In concluding this Section we want to verify that 
Eq. (42) goes into the equivalent equation for the 
plane case when R->oco. More exactly, we want to 
show that the same operation which transforms, in 
the limit Roo, the curved beam into the plane one, 
transforms S, into the x-component of Eq. (20). This 
can be done, provided we use in Eq. (20) a solution of 
Eq. (12)—where 6 (x) is replaced by 6 (2— R)—which 
doesn’t propagate in the region x < R, as in the cylindri- 
cal case; for instance, this could happen because of 
a “mirror” situated in an appropriately chosen plane 
x=a<R. Such a solution, valid for x>a, is 


J m (x) (im =i, 2, 3,...). 


B.(x)=b,[—e' 9 + iV P(x 
+b,0(a— Ryei VK*-Px_ (48) 


When the constants b, and 6, are determined as in 
Section 3, we obtain: 


2x1 


(x) cV 


{(V.K2—P— ila) ei 


+ (V + ila) e~ iV PRG — Ry ei 


(49) 


For the 2-component of Eq. (20), calculated on 
the side of the beam looking towards the positive 
x-direction, we obtain the expression* *: 


* In these conditions a density-modulated beam of 
electrons cannot be realized and in order to have a 
deuteron beam we must assume m of the order of unity. 


** Note that in the resonance regions of 7,2, when 
w>V, Eq. (50) takes the same form as Eqs. (44), (45) 
and (46) and that, when we calculate the mean value of 
these expressions as we did in Section 4, we again find 
the radiation equal to the sum of the radiation emitted 
from the two sides of the beam of Section 3 in the corres- 
ponding situations. 
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c Vu? —V?2 


x {v w? — V2 cos [ V/wt—V2 (R—a } 


/ 2 
+aV sin| Vu v2(R—a)}} (50) 
Now in order to transform the curved beam into 
the plane one, we put: 


R=mll; Weylinder = Vplane = constant. 


Then, in the case w=V, we note that [7] the 
functions J, (m) and (m) are positive decreasing 
functions of m and are given (for a sufficiently large m) 
by the Cauchy formulae: 


1 2 
‘(1 3 26 
23.38 2* x 
We then get from Eq. (42) 
S? (R)~ 0,2 ( v3. (52) 


Hence, in the limit, S? (R) has the resonance 
at w= V as in Eq. (50). However, in order to see how 
effective such a resonance is, let us consider the mean 
relative energy loss of a deuteron per cm of path, 
given by Eq. (33). For instance, in the case o> 1 and 
x sufficiently smaller than 1, so that « <1/u, the mean 
relative energy loss is 

1 
122 Nrja?m* < 
E 

For a beam with N = 10* deuterons per cm?, the 
value of this expression is <10-* cem-'. 

Hence, the resonance found in the plane case for 
w= V is practically ineffective. Measurable radiation 
can then only be emitted in the resonances of n,?, 
as we shall see in the next Section, when the finite 
conductivity of the plasma is taken into account. 

In the case w< V, we substitute in Eq. (42) the 
inequality [7] 
mw 


)<m where <1. (54) 


which follows from Carlini’s formula. We then have, 
according to the Cerenkov condition, 


Im 


lim S?(R)=0. (55) 


co 


On the contrary, if w> V, we introduce in Eq. (42) 
the Langer uniform expansion [8]: 


1 
mw 1 3 


- 4 


(z) cos (z) sin 4 + (56) 


where 


1 
w= 2=m(W—tan—W). (57) 
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Kk. CANOBBIO: 
If z>1, J: (z) is given by Eq. (43) and Y, (z) ~ 
3 3 
1 


~(2/xz)* sin[z —(5x/12)]. Then, by equating the 
so-obtained expression from Eq. (42) with Eq. (50), 
we have an equation for the mirror position a. 


5. On the peak value of S 


When the finite electric conductivity of the plasma 
is accounted for, the expressions for 8, given in 
Sections 3 and 4, are no longer appropriate. The 
squared refractive index n? is a bounded complex 
function [6] whose real part in the resonance regions 
differs very much from n,?; moreover, the imaginary 
part of n indicates energy is absorbed by the plasma. 

It is not our intention to study the problem in its 
most general form. In what follows we shall restrict 
ourselves, for simplicity, to the evaluation of the 
peak value of the radiation from the plane and from 
the curved beam at the “ion-resonance” frequency. 
This is practically the only resonance frequency of 
interest for experiments. Finally, the resonance 
breadth and the attenuation distance of the waves 
will be briefly examined. 

Consider now the general forms of n? and « given 
by Eqs. (13) and (14), respectively, and express 
them in terms of the dimensionless parameters x, ¢ 
and yw introduced in Eq. (25) and of the new para- 
meter tT=y/V we. By supposing t<1 and denoting 
with x, the value of x which corresponds to the 
ion-resonance frequency in absence of collisions, it is 
easy to find that Re (n?) and Re («) reach their peak 
values at x =x,—O (rt), where O (7) is positive and of 
the order of Tt. 

In the following we shall choose, in order to simplify 
the calculations, just the value x =x,—O (r) + O (7°). 
At this value we have (see Fig. 3) 


Re (n?) = Im (n?) = n’, 


2n* 


-7* 


Fig. 3. Schematic behaviour of n? in the neighbourhood 
of x. The behaviour of « can be obtained by substituting 
« for nr? in the ordinate. 


so that we can write 


n2= n2V 2eit/4 
(59) 
a =a V2eir/4, 
#r—O (t), nand « are real functions of 0; approximate 
expressions are given in Table I. It follows from 
Table I that «> n. 

With these assumptions, the 2-component of the 
complex Poynting vector, Eq. (20), for the plane 
beam, 

where dy is the derivative with respect to x, acting 
only on the continous parts of B-, reduces to 


(58 c B.* B 
Re («) = Im Im{ =}. (61) 
TABLE I. 
1 o a 1 
| 1 1 1 
| 
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By using Eqs. (18), (59) and (61) and supposing 
moreover w> V, the x-component of the complex 
Poynting vector, calculated on the side of the beam 
looking towards the positive «x-direction, becomes: 


cos (x/8) sin (x/8) 


Sz = (N . (62) 


On the other side of the beam it is SX ——S?. 
In the same way, the radial component of 8, for 


the curved beam, takes the form (61) where dy is 
replaced by d,. Remembering now that [8] 


(KR)]*=H2(K*R) (63) 


and that for |K R| >m the Hankel functions be- 
come [8]: 
1 


while J» (K R) is given by Eq. (43), from Eqs. (38), 
(59) and (61) we get the following expressions for S, 
on the two sides of the beam, 


a? cos (7/8) + sin (x/8) 


S? —8&, 2 (Ne)? n3 21/4 


(65) 

For simplicity, we indicate with S the common 
value of S7 and S? (R). Note here that the mean 
energy flux, Eq. (62) or (65), does not depend on the 
velocity of the beam particles and on the beam 
curvature. This means that, when the interaction 
between beam and plasma becomes very strong for 
a particular value of the modulation frequency, the 
large modifications induced by the beam in the 
surrounding plasma react on the field radiated by the 
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Fig. 4 Schematic behaviour of «?/n* as a function of By. 


beam in such a way that the resulting field comes out 
to be essentially determined by the plasma properties 
and no longer by the beam dynamics*. The dependence 
of Eq. (65) on the plasma properties is better dis- 
cussed by expressing «*/n* in terms of the magnetic 
field By (measured in Gauss), the particle density of the 
plasma d and the plasma temperature 7' (measured 
in°K). By using for y the expression given by SprrzER** 
[3] we get Table II. 


* Something similar happens in the problem of the 
energy loss of a fast particle traversing matter, see 
e. g. E. Fermi [9]: here, the interaction between particle 
and the medium becomes very strong when the particle 
velocity v approaches c; but then the radiation doesn’t 
increase indefinitely ; on the contrary, for velocity larger 
than a certain v, it takes a value which is independent 
of v and is determined by the dielectric properties of the 
medium; such a value diminishes when the density of the 
medium increases. Note in this connection that our 
expression, Eq. (24), for instance for a non-dissipative 
plasma and at low frequencies becomes, when w-—c, 
oc B,/V d, where d is the particle density of the plasma. 


** Remembering that in our cases the current is 
perpendicular to B, we have to choose the y used for the 
transverse conductivity. 


TABLE I]. 


d<10°- B,? 


d > 26-B,? 


. 1 
; 2 
Oo, = 2 4,8- 108. (B,2 + 3,8- 10-2. 8.7 - 10°. Byd 
sec (By? + 10-5d) 2 
3 2 8 
a? 102. B,? 2-10°. Bo? 
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Ion-electron collisions tend to destroy the resonance 
condition; hence it is clear that the mean energy flux 
must increase with the temperature. Then, from Fig. 4 
where the behaviour of «?/n® taken as a function of By, 
is sketched, we see that the peak value S becomes 
particularly high when the plasma frequency @, is of 
the order of Vmjwe. From Fig. 5 on the contrary it 


B o = constant 
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83Bo d 
Fig. 5 Schematic behaviour of «?/n? as a function of d. 


follows that the plasma has also a kind of screening 
effect on the beam radiation which increases with 
the ratio wyp?/wjmex d/B,?. The greatest value of S 
occurs when d~8,3 B,?: 


3 
4 
S ~ 0,6 - 10? 3 (Ne)? = ; 
(InA) ® By? 


(66) 


In this case the mean relative energy loss of a deute- 
ron per cm of path, for the plane beam as well as for 
the curved one, has the value: 


(InA)2 B 

In many cases of interest for experiments, however, 

o>1/u; then from Table II we can write 


~1,5-10-"4 = (67) 
== w 
2 
0 


24 
(In A) 
and, for deuteron beams, 
7 3 
2 4 
~ 0,6 - 1074 = (69) 
When B,=104 Gauss, VN = 108 em-?, d= 10% 
and 7'= 10° °K, Eq. (69) becomes 


The particle velocity w is, in this case, ~0,8-107 A 
and the beam radius R~A/2 x, where / is the wave- 
length of the modulation. 

If A is of the order of 10 cm, so that R~1 or 2 em 
and w~8-107 em sec~!, the value of Eq. (70) is ~1. 

Of course, the practical interest of the previous 
equations depends essentially on the value of the 
“breadth” of the energy release, taken as a function 
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of the modulation frequency. This breadth can be 
evaluated by looking at the form of n? and «. For a 
displacement Ax ~—O (rt) away from the resonance, 
Re (n?) and Re («) become ~O (r)-times smaller than 
n* and «, respectively, while they go at once to zero 
for a displacment Ax~O (rt); moreover, except in 
the neighbourhood of x,, it is always Im (n?) < Re (n2) 
and Im («)<Re (x). Consequently, the value of § 
(and of d92/E£) becomes negligible compared to S 
(and to 6 £/E) for displacements Aw of the modulation 
frequency away from the resonance such that: 


(71) 

Approximated expressions of Eq. (71) are given 

in Table II. In cases where Eq. (69) holds, expression 
3 1 


(71) takes the form: ~1-10-3[dIn A/(T?* B,)|*; 
when Eq. (70) holds, it is |Aw|/w~3,7-10-%. 

In concluding we note that the energy released by 
the beam is absorbed by the plasma over a distance, 
whose order of magnitude L is given by [3] 

z = Im(*"). (72) 

For the explicit form of L, see Table Il. When 
Eq. (70) holds, ZL is about 10-?cm. L increases 
with By, while it diminishes when 7 increases, be- 
cause of the fact that the collisions tend to destroy 
the resonance and thus deepen the penetration. 
(The wave attenuation due to the collisions is negli- 
gible compared to that due to the absorption at the 
resonance. ) 
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THEORY OF CERENKOV AND CYCLOTRON RADIATIONS IN PLASMAS 


Taro KrHara, Osamu Aono, Ryo SUGIHARA 


DEPARTMENT OF Puysics, UNIVERSITY OF ToKYO, ToKYO, JAPAN 


Radiation from a charge q moving in a helix in magnetoplasmas is investigated theoretically. 
When its speed v is much larger than the thermal velocity, (m~! kT)!/?, of the plasma electrons 
and the gyration frequency is much smaller than the plasma frequency w,, the radiation power 
from the charge is (q* @,?/2v) In (v?/m~! kT). Cyclotron radiation from an electron with non- 
relativistic speed decreases to zero as plasma density increases. For a positron in a dilute plasma, 
however, the radiation is strengthened. This strengthened radiation from a positron again de- 


creases with increasing ,?/@y? and becomes zero for 


»2 (wy = gyration frequency of the 


plasma electron). Damping of the Cerenkov radiation due to collisions of plasma electrons is 


also discussed. 


1. Introduction 


In this paper we investigate theoretically the radia- 
tion from a fast moving electron in a uniform infinite 
plasma. The interactions of the electron with the 
plasma can be classified into two parts: the microscop- 
ic one (close and distant encounters) and the macro- 
scopic one (resonances with the plasma as a continuous 
medium). 

The distant encounter is characterized by the impact 
parameter which lies between e?/mv? and the Debye 
length lp. Here —e, m and v are the charge, mass 
and velocity of the moving electron. The impact 
parameter plays an essential role in such irreversible 
processes as temperature relaxation and _ electric 
conduction. A fast moving electron loses its kinetic 
energy mostly due to the distant encounter with 
plasma electrons. The time rate of energy loss of this 
type is given by 

ea," mv lp 
-In (1) 
where @, is the plasma frequency. 

The close encounter with ions causes the brems- 
strahlung. 

The energy loss due to the Cerenkov radiation from 
the moving electron and the shielding effect on the 
cyclotron radiation are two examples of resonant 
interaction characterized by the macroscopic electro- 
magnetic field of wavelength longer than the Debye 
length. 

A charged particle moving uniformly along the 
magnetic field in a plasma radiates electromagnetic 
waves [1, 2, 3], called the Cerenkov radiation [4, 5, 6}. 
A charged particle gyrating in a magnetic field in 
vacuum radiates electromagnetic waves, the fre- 
quencies being close to the gyration frequency and its 
higher harmonics [7, 8, 9]. This is the cyclotron 
radiation. For a charged particle moving helically in 
a magnetoplasma these two types of radiation can 
not in general be treated separately, although for 
higher harmonics of the cyclotron radiation the 


surrounding medium can be approximated by vacuum 
{10}. Theories of radiation in a magneto-plasma are 
mostly formal and give few concrete results [3, 10}. 

A purpose of the present article is to show that under 
certain conditions we can treat these two radiations 
separately and thus obtain useful theorems. We use 
the Fourier-expansion method developed by StrTENKO 
and [12]. 


2. Preliminaries 
Starting from the Maxwell equations 
1 @ 
rot E(r,t)= — H (r,t), 
(2a) 


rot H (r,t) = D (r,t) + J (r,t), 


we make the Fourier expansion for each quantity, 
in such a way as 


E(r,t) = E (k, exp (ik-r—iwt)dkdo, 


(2b) 


E( 


Then for the Fourier components the following equa- 
tions hold: 


k, — = E* (k, ), etc. 


ikxE='°H, ikxH=— "D+ (2c) 


On eliminating H, we have 


o 


D+ (kk-E—#E) =—i (3) 


Using the connection between the Fourier compo- 
nents of the electric displacement and the electric 
field intensity 


D (k, m) = €(m)- E(k, 


(4) 
(—o) = (om). 
we obtain 
T.E=—i**J, (5) 
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where 
Ti = n? %k — + (6) 
ke no 
k= (7) 


By use of the inverse tensor T-! the solution of 
Eq. (5) can be given in the form 


(8) 


o 


We have therefore 


E(r,t) =—idn{ Jexp (ik -r—iot)dko da 
and hence (9) 


[nx x (10) 


The integration variable can be transformed according 
to 


dkw "dw = wda (11) 


where dQ is the solid angle element of the direction 
of the wave vector k. 

We consider the electromagnetic field produced in 
the medium in which a point charge q is moving with 
a velocity v(t). Denoting by ry (t) the position of 
the charged particle at timet, the current is represented 
in the form qvé(r—r,.) with the Dirac 6-function. 
Then, its Fourier component is calculated to be 


3 (k, w) = [¥exp(—ik- te + (12) 


The energy radiated from the charge is given by a 
surface integral of the Poynting vector S=(47)-!c¢ 
(Ex H). The Poynting vector, however, is compli- 
cated in general, and the calculation is performed only 
in particular cases (see Section 5). 

On the other hand, the time rate of energy loss by 
the moving charge, —d W/dt, due to interactions with 
the surrounding continuous medium can be calculated 
by use of the relationship 


dw 
dé = E (fe, t) 


| .Jexp doo , 
(13) 
in which the field E is taken at the point where the 
charge is located. This energy loss is not always 
equal to the energy radiated from the charge even if 
the electric resistivity of the medium is negligible. 
The charged particle, in some cases, does work on the 
medium in the neighbourhood of its track without 
radiating energy. There are two criteria in this respect. 
The condition that absorption is absent is that the 
tensor ¢ is hermitian. 


The time rate of radiation from a moving charge is 
equal to the time rate of energy loss given by Eq. (13) 
when the dielectric tensor €(m) is hermitian in the 
domain of w which contributes to the integral of Eq. (13). 
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An example of energy loss without radiation is given 
in the next Section. 

Plane electromagnetic waves which are propagated 
in a source-free magneto-plasma are given by Eq. (5) 
with J=0: 

T-E=0, (5a) 


the index of refraction, n, being given by the relation- 
ship 
det T= 0. (15) 


There exist for n? two finite, roots n,? and n,?. (These 
roots are real when the tensor €(@) is hermitian.) 


When these roots are real and when the integral with 
respect to n® is determined by the integrand at a positive 
root or roots, then the energy loss is entirely due to the 
radiation with these indexes of refraction. 


When, in particular, the particle velocity v is a 


constant, the Fourier component of the electric cur- 
rent, Eq. (12), reduces to 


j- fexp (—ik - vf +- ict) dt 


= 
v l 


the particle position at ¢=0 being chosen as the 

coordinate origin. Then, the time rate of energy loss, 
(13), becomes 

i -vo(" , 


(17) 


where Eq. (11) has been used. This expression was 
first given by Strenko and KoLomenskti [12]. 


3. The polarization loss 


As a preparation we first consider a collisionless 
plasma with no magnetic field in which a charged 
particle is moving. We assume that the speed v of 
the charged particle is much larger than the thermal 
velocity of the plasma electron, 


In(v?/m—kT) > 1, (18) 


where k is the Boltzmann constant, 7' is the plasma 
temperature, and m is the electron mass. 
For a scalar dielectric constant ¢, we have 


the z-axis being taken along the velocity vy. The energy 
loss, Eq. (17), thus takes the simple form 


aw 


(20a) 
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(18) 


sma 


(19) 


a) 


in which P= v/e and @ is the angle between v and x. 
Hence 
= | dnwdm. (20b) 


n>p 


The plasma can be treated as a continuous medium 
with the dielectric constant 


e=1— (21a) 


for wavelengths longer than the Debye length /p, 
namely for 


Bn < Ip). (21b) 


The limit v/(\@! Ip) is a large number by virtue of our 
assumption, and the integral with respect to Bn can 
be cut off at this limit. The integral of Eq. (20b) is 
determined by the values of the integrand at its 


joles +m, and we have 
I 0 
dW 
dt 
1/8 
24) 2 
2v 


This is the polarization loss mentioned by FErmi [13] 
and Pines and Boum [14]. 

For an electron the energy loss, Eq. (22), due to 
resonant interaction with the plasma as a continuous 
medium is less than one half of the energy loss Eq. (1). 

fin IpkT 


In (la) 


y2 
mo 
which is due to encounters with plasma electrons as 
particles. 

The energy loss is not at all equal to the radiated 
energy from the charged particle, which will be found 
in the following to be zero. In fact the dielectric 
constant ¢ is not hermitian at = where collision 
effects can not be neglected. 

By use of Eqs. (9), (11), (16), and (19) we can cal- 
culate the electric field at a position r. On choosing 
the z-axis along the particle velocity v we let 


v = (0,0, fe), r= (0,0, 2), 


(23) 
= (sin6 cos sin cos§) . 


Then we have 


E:(r,t) =— iff = cos6 —1) 


e(n? — 


exp i (gsinBcosy + zcosf) 
qv. cos 
J sin 0) sin db n?drnoda 


qi 1— pre 
laa Jo( V pn? —1 
n 


(24a) 


CERENKOV AND CYCLOTRON RADIATIONS 


On taking V /?n?—1 as an integration variable, we 
obtain 

1— Be 


Ky (ho)e oda , (24b) 

1 
where h = (1 — Be)? w/v and K,, is the modified 
Bessel function of the second kind, for which the 
following formula holds: 


| = Ly? Jn (ax) dx = y" Ky (ay). 


The g-component of E is zero. 
As regards the magnetic field, its @-component is 
calculated to be 


H,(v.t) Ky (ho) da. 25) 


The z- and 9-components are zero. 

The power radiated from a cylindrical surface which 
is at the distance o from the particle track is given by 
the surface integral of the Poynting vector 


Cc 


| S, | —E-H,2nodz 


= ho Ky (ho) Ky den, (26) 


Rv 


on 


which is calculated to be 
(S,2ngdz = K,( Ka ). (27) 


This result is in agreement with Eq. (22) when we let 
o be lp. However, it decreases exponentially for large o, 


S,2nodz ~ (28) 


and hence radiation is not emitted beyond v/m,. The 
energy loss Eq. (22) is not in the form of radiation but 


in a kind of work QE-dD done in the vicinity of the 
particle track. 


4. Cerenkoy radiation in collisionless plasmas 


In this Section we consider the energy loss of a 
moving point charge due to the Cerenkov radiation. 
Assumption 1. The collision frequency is negligible 
compared with the gyration frequency my of the 
plasma electron; but the latter is not much larger 
than the plasma frequency @p. 


On (29) 


Assumption 2. The speed v of the moving point 


charge is much larger than thermal velocity of 


the plasma electron, 


In kT) > 1, (30) 
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where k is the Boltzmann constant, 7’ is the 
plasma temperature, and m is the electron mass. 

Assumption 3. The radius of curvature of the 

orbit is much longer than the length passed in a 
period of plasma oscillation. 

The last assumption enables us to consider the 
particle velocity as uniform in the time interval 1/09, 
a characteristic time for the Cerenkov radiation, and 
to use Kq. (17). 

The z-axis being chosen in the direction of H, the 
dielectric tensor takes the form 


& O 
€= 0 (31) 
0 O & 


where ¢,, € and ¢, are real for a collisionless plasma. 
The determinant of the tensor 7' is calculated to be 


det T =n‘ [e, sin? y + €, cos? y]—? [e, €, (1 + cos? y) 

+ (€;?—e,?) sin? y] + €3, (32) 
where y is the angle between H and x. Let the compo- 
nents of x be 

x= (sin y cos «, sin y sin %, cos 7). (33) 


Then we have 
v-T-!. v = (det T)-! v?{ cos? 6 
—n*[e, (1 + cos?@) + (es—e,) sin? (1 —sin? y sin? «) 
+ €, €,8in® O + (e,2—¢,") cos?O}, (34a) 


where @ is the angle between H and vy, and @ is the 
angle between v and x so that 


cos 6 = cos O cos y—sin O sin y cosa. (34b) 


The integral with respect to Pn can be cut off at 
v|(\| lp), and the total integral is mostly contri- 
buted by the region of large fn, where det T can be 
replaced by the first term in the right-hand side of 
Eq. (32) and cos? @ is very small because of the factor 
6 (nf cos0—1). Moreover, the nonvanishing integral 
with respect to » comes from the point where det T 
is zero. From the inequalities 


sin? y +e, cos*y|  (34e) 
it follows that 
€, (1 + cos? 0) + €,)sin®@ (1- sin? ysin?«), = 
+ cos*y + 2 cos cosy + €,c0s6| 
(34d) 


and the right-hand side of Eq. (34a) is approximated 
by its first term. Thus we obtain our basic relation 

iW ge 29 

i [dean | cos? y 

—oo 0 
x 6 cos —1) n? (35) 
Nm lp) . 


As regards the integration with respect to the direction 
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of x. we introduce azimuthal angle q of * around y 
so that 


d2=sin0 cos y = cos cos 6 + sin sin cos g, 
(36) 

@ being the angle between H and y. 

The three components of the dielectric tensor are 
given by 

for wavelength longer than the Debye length. The 
frequency of the dominant radiation is given by the 
inequalities (34c), namely 


| (@?— wy? — wy”) + COS? y | 


| <€1,cos*9<1, (38) 


| (w? — wy?) 
Hence the domain of w? is composed of two parts 


0 <w? < Min (wy?, 


In particular, 
+ forvx H=0, 
forv -H=x-H=0, (40) 


forv - H=0, 


(See Fig. 1). The set of intervals, Eq. (39), for the 
radiation frequency was pointed out by KoLoMENSKI 
[2] in the case of v parallel to H. 


x 


for 


Yous +w? 


Fig. 1 Frequency of the dominant Cerenkov radiation, 
see Eq. (39). 


The dielectric tensor €(@) is hermitian excepting 
two points, wm? =0 and wm? =m", where the integrand 
of Eq. (35) vanishes. Therefore, according to one of 
the criteria mentioned in Section 2, the energy loss 
Eq. (35) is entirely due to radiated power. 

On substituting Eq. (37) for ¢, and e, in (35) we 
consider the integral 


i — wy") dw? 


(41) 


Let 
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Let us write the denominator in the form (j@?—«’?) 
for which 


0 <w? < Min (my?, , 


Max <0”? + on? , (42) 
Then we have 
‘2 2 “2 "2 2 
I=. (on on) (43) 


the first and second terms on the right corresponding 
to the first and second parts of the domain, Eq. (39), 
respectively. This expression is equal to w” + @”"*—wy” 
or ,?, and we finally obtain 


"m 


nee | cos? 95 cos? —1) n? sin9 dn 
0 


0 
Ip 
v n 2» 
1/8 


We have tacitly assumed that the dominant fre- 
quency of the Cerenkov radiation in a collisionless 
plasma is of the same order of magnitude as the plasma 
frequency @ ,. This is in fact true, because the first 
term of Eq. (43) is negligible when wy?<w,?. Thus 
we reach the following conclusions: 


1. The power lost in the form of Cerenkov radiation 
is given by Eq. (44). It depends neither on the 
absolute value nor on the direction of the 
magnetic field H; it is just the same as the 
polarization loss, Eq. (22), in a plasma with no 
magnetic field. 

2. The Cerenkov radiation is emitted in the direc- 
tion almost perpendicular to the particle velo- 
city. The power is uniform with respect to the 
azimuthal angle around the particle velocity. 

3. The frequency @ is in the domain Eq. (39). The 
radiation in the first part is negligible if wy? < w,?. 


5. Collisional damping of the Gerenkoy radiation 


Having discussed the Cerenkov radiation in a 
collisionless plasma, in this Section we consider 
damping of the radiation due to finite electric con- 
ductivity of the plasma. 


Assumption 1. The collision frequency, eon, is 
still much smaller than the gyration and plasma 


frequencies : 
Meo” < Min 9”) (45) 

Assumption 2. 
<€v<c?. (46) 


Assumption 3. The velocity vy of the moving par- 
ticle is parallel to the magnetic field H. 

The first assumption will be taken into account when 

we make expansions in @con/m. The first half of the 
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second assumption corresponds to the relation that 
the wavelength is much larger than the Debye length; 
the second half enables us to treat large values of the 
refractive index only. For velocities which are not 
parallel to the magnetic field, calculation would be 
formidable. 

We start from the expression, Eqs. (9) and (10), 
for the electric and magnetic fields with the electric 
current given by Eq. (16). On choosing the z-axis 
along the particle track, we let 

v = (0,0, fe), r = (0,0, 2), 
(47) 
= (sin cos sin sin g, cos 9) . 


Because of the 6-function large values of n need to 
be considered, for which 


= 


xB cos 
€18in? 6 + e, cos*@ 


6(nBcos6—1) 


exp (i cos 9) db dg n?dnodo. 


(48) 
Its z-component is calculated to be, Eq. (24b), 
= K, (ho) mdm, (49a) 
where 
&3 w? 
h? = -* Reh>0O, (49b) 


and K,, is the modified Bessel function as in Section 3. 
The g-component of E is negligible. 

As regards the magnetic field H, its g-component 
is found to be 


q Be, sin 6 
H, (r,t) or | cos? 6 (n B cos6—1) 


x 9 sin 0) dO6ndnoda 


c 


= [a 


(50) 


The power radiated from a cylindrical surface which 
is at the distance 9 from the particle track is given by 
the surface integral of the Poynting vector, ef. Eq. (26), 


h* 0 K, (h* 0) Ky (ho)oda. (51) 
& 
For collisionless plasmas, the domain 
<Min 
Max + on’, (52) 


in which h is purely imaginary, contributes to the 
radiation. Regarding the collisional damping, we use 
the asymptotic formulas for large o. 
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1 
Ky (ho) ~ (z/2ho)* e-"e, 


= 
K, (h*9) ~ Me. (53a) 
Then 
nh* 


(53b) 


because of the inequality Re h< |Imh| . Hence 


[Se de = 2eRehanda, (54) 
1 

where + is to be taken so that Re(-+1/e,)>0. 

Considering the relation (4), &,(@)=«,* (—w), we 
finally obtain 

1 

|S, 2modz = — (55) 


w>0 


By use of an effective collision frequency @coy let 
us adopt the expressions 
[(@ + imeon)? — on?) 
(56) 
For a small @eou the real part of h is calculated, in the 
domain (52), to be 
1 
ok _ coll — on?) 


w? (a>? + oR? — ow?) 
| — we? (@? — wy’) (am? + — w?) 


The method of steepest descent (Sattelpunkt- 
methode) enables us to calculate Eq. (55). It is found 
that Reh takes its absolute minimum near the 
midpoint in the first part of the domain (52): 


0 < w* < Min (@y?, . (58) 


In fact, the values of wm? which minimize Re h are the 
following : 


1 
wo = 5 for on? 
2 2 2 1 2 
w = 0.479 for oy? = Mo 
1 
: 3 on = 3 for wy? = Og” (59) 


= 0.574 2 


= 0.560 for my? = 3 


for wy” = 2 wy 


1 
= for wn? > 


The expression of Eq. (55) is calculated to be 


1 
[Se2nede exp (— 4 @con 


for on < 9" (60a) 
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= exp (— 4 @cono/V3 v) 
Moll 
(60b) 
for on 
1 
(60) 


for wy? > 


6. Non-relativistic cyclotron radiation in collisionless 
plasmas 


A characteristic feature of the non-relativistic 
cyclotron radiation from a point charge q in a colli- 
sionless plasma is that the variation of k-r. in Eq. (12) 
is negligible (for k not too large) in comparison with 
the variation of wt: 

q 
On choosing the z-axis in the direction of the magnetic 
field, we let 
= 
\q| (62) 
V=(+2%,8in v, cos V2) , 


where m, is the particle mass and the upper and the 
lower signs of + are for positive and negative q. From 
Eq. (61) we then obtain 


+q eio,t —e—iot 


Jy = £2, oy) + 


The Eq. (13) with Eq. (11) is similarly reduced to 


(64) 
By use of the vector 


a= (ats, dz) = (+ (65) 
od in Eq. (64) can be written as 
since wd (m)=0; and we have 
imi ¢ i | - ae~i™ 
n?dndQ. (67) 
On taking the time average over a gyration period, 
and making use of the relation, see Eq. (6), 


a-T-1(—«,)-a* (69) 


d 
é Supt 
whe 
ig 
for 
| Th 
as 
wl 
0 
is 
: 


less 


stic 
vith 


61) 


etic 


62) 


the 
om 


63) 
to 


64) 


65) 


66) 


67) 


od, 


68) 


69) 


we finally obtain 


00 


(70) 


where y is the angle between the magnetic field and 
wave vector. 

When the wave vector is taken in the yz-plane, the 
tensor T is of the form 


ig, 0 
T=|-ie, -n®sin?y+-e, n®siny cosy 
0 n’siny cosy —n*sin? y+ (71) 
for which we get 
a*-T-!- adet T 
= [— n? (e, sin? y + cos? y) + &3] 
+ 2 dy (n® y — (72) 


+ ay?( — n? + €,) (—n? sin? y + €) . 


The discriminant of the right-hand side of this equation 
as a quadratic form of |a,| and ay is just the det T 
given by Eq. (32). Let 


det T = (e, sin? y + €, cos* y) (n? — (n®? — ng"), (73) 


where n,? and n,* are real, then for n? = n;? (j= 1,2 
the right hand side of Eq. (70) is 


[— nj? (e, sin? y + €, y) + + ajay)? 


74 
=[—nj%(e, sin®y + +e,e,) (1a) 


where 


— nj? sin? y + €5) 
— nj® (e, sin? y + €, cos? y 
| + &3COS* y) + &3 


(74b) 


On carrying out the integration with respect to » 
in Eq. (70), we have 


dw @v,?a,? df 1+ a;\2 
dt 3 2 
0 j=1,2 
| — nj? (e, sin? y 
| (e,sin?y 


2 

cos*y) + 
Re nj \sinydy. 
+-€3 y) 14”) 


(75) 


This integral comes entirely from real positive 
poles of the integrand of Eq. (64) or from real positive 
roots of det T=0. Hence, according to the second 
criterion in Section 2, the energy loss, Eq. (75), is 
due to the radiation with these indexes of refraction. 

The ¢,, €, and ¢, are given by Eq. (37); n,? and 
n,* are given by Appleton-Hartree’s formula (see, 
e.g., [15]): 
nj? = 1 — 2erg? (w? — /{2 w? (@? — Wy”) — Sin? 


+ (-1))-! y + 4 @? (@?— cos? y] *} 
(76) 
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where (—1)/-!=1 
(—1)-1= 


for the ordinary wave and 
—1 for the extraordinary wave; and 

— nj? (e, sin? y 
(e, sin? y+ 


+€3 cos? y) 


1 


KP’ 
1 
K;=-” oysin® y+ y+ 4 cos?) 2 
2 (w*— w,”) cos y 
(77) 
Let us first consider an electron (q=—e) which 
is gyrating in a dilute plasma under the conditions 
mIkT < we? < won?. (78) 
The gyration frequency @, is then given by 
= on? (1— v?/c?) , (79) 
and furthermore, for 


K,=I1/eosy, K,=—cosy, 


n,=1 ny” 

ssin?y + 1 

ssin® y — cos? y 


ssin*y +1 


8 = (c?/v?) (80) 


By use of these relationships, the radiated power is 
found to be 


iw 92 2 9,2 
— a= fe (8) (Sla) 
s+l 4 1 + cos? y 
fe (8)= y+ )? 4(ssin? 1)? sin y dy 
0 
(8+ 2)V¥28+1 
= 3(e+ 1? (81 b) 
The function f,(s) decreases from f/f, (0)=2/3, the 


well-known value for vacuum, to f,(cc)=0 as 
the plasma density increases. (Fig. 2). For «,? not 


{ 


0 5 10 


Fig. 2 Cyclotron radiation from an electron, /,, and 
from a positron, fp. see Eqs. (81), (82). 


much smaller than my”, no radiation is emitted from 
an electron since «, = %, = 1, as previously pointed out 
by GryzBure [6]. 
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1.54 


1.0 4 


0s 4 
2 
0 1 , 2 
Fig. 3 Angular distribution of the cyclotron radiation 


from a positron in the non-relativistic limit, see Eq. (83). 


It is interesting to compare a positron (q= +e) 
with an electron (q=—e). The radiated power from 
a gyrating positron under the condition (78) is caleu- 
lated to be 


fp (8), (82a) 


fp(s)=2 + 2tan-*?V2s+1 3 
(82b) 
The function fp(s) increases from /, (0)=2/3 to 
fp(co)=2+7/2 as the plasma density increases 
(Fig. 2). For m,? not much smaller than wy?, the 


Si pay 


2 


0 1 2 
Fig. 4. Integrated cyclotron radiation from a positron 
in the non-relativistic limit, see Eq. (83). 
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cyclotron radiation power is given by 


=a | (88a) 
0 j=1 
where 
F Re nj 
i(y)= (83b) 
o=OH 


The integrand )"F; (y) is shown in Fig. 3 as a function 
of @»?/wmy? and y. Fig. 4 shows the integrals | F (y) dy 
which take, in particular, the following values: 


(y) dy =2—nx/2, (y)dy=n 


for v?/c? <1; 


(y)dy=0, (y) dy =2 


for =1; 


fF: (y)dy=0, | F,(y)dy—90 


for > 2. (84) 
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WUAKIOTPOHHOE HOHOB B NJIASME 


B. U. TIMCTYHOBHY, B. Jl. LIA®PAHOB. 
OPJIEHA JIEHHHA UHCTUTYT ATOMHOUM UMM. HU. B. KYP4YATOBA 
AKAJIEMUM HAyYK CCCP, MockBA, CCCP 


B paGote onpeneseHa MHTCHCHBHOCTE H3J1Y4eHHA ObICTPbIX MOHOB B XOJIONHOH Ma3Me. 3aa4a NA 
B CBA3H C Ha ycTaHoBKe Orpa B PexKHMAaX Pe3IOHAHCHbIX MMKOB 
Ha MOHHOM YacTOTe ee OGepTOHax. XOTA yKa3aHHbie 
HE B BOJIHOBOH 3OHE, ECTECTBEHHO YTO OGHAPyxKeHHAaA B IKCNEPHMeHTe 3ABHCHMOCTE 
MMKOB OT MIOTHOCTH XONOAHOM M1a3Mbl CBOe OTPaxeHHe M B MHTCHCHB- 
HOCTH H3J1y4¥eHHA. PacyeTbl YTO CKOPOCTH HOHOB MaKCHMYM 
MHTCHCHBHOCTH H3JIYYCHHA, CMeLUaeTCA B CTOPOHY BbICOKHX 4YacTOT. 

STO CMeLUIeHHe CMeCLUCHHIO MaKCHMYMa MHTCHCHBHOCTH TPH CHHXPOTPOHHOM H3J1y- 
WIEKTPOHA, CKOpOCTb, K CKOpOCTH cBeTa. B MOHOB B M1a3Me CKO- 
POcTH cBeTa CKOPOCTb BOJIH B M1a3Me, KOTOPaA B 
pacyetos o6nactu K CKOpocTH c4 = By/(42 mj no)'!?. Tloxvromy seicoxne oGepTousi 


MPOABIAIOTCA MPH CpaBHUTeNbHO CKOPOCrH HOHOB c. 


1. 


B Hactosusei paSote 
MHTCHCHBHOCTH OBbICTpbIx 
HOHOB, B Haxous- 
B MarHHTOM Movie. 06 HOHOB 
BO3HMKJ1a B CBA3H C Ha ycTaHoBke Orpa 
pe30HaHCHbIX MMKOB HallpsxKeHHOCTH 
Ha YacToTe ee OGepToHax [1]. 
(PexuMbI, B KOTOPbIX HHTCHCHBHbIe 
HalipsxKeHHOCTH Ha OGep- 
TOHAX HOHHOH WHKJIOTPOHHOH YacTOTbI, 
Hatllycke B KaMepy Orpbi aproHa 
p= 10-*-10-° mm. pt.ct.). XoTa ykKa3aHHbie 
He B 30He, TaK 4TO 
mia OObACHeHHA HeOOXOAMMO paccuHTaTb 
KBaapaT 
HOH 4YaCTH NOMA, 
TeM He MeHee eCTeCTBeHHO B Ka4ecTBe NepBoro wara 
pewiMTb Gosee Npoctylo MHTeH- 
rHOCTHKH [Jla3Mbl B YCOBHAX NpH 

MOHOB B B Ma- 
THHTHOM ose, OTMMYAeTCA PADOM 
CBA3AHHbIX C BeIMYMHOH NoOKa3aTeA Mpe- 
JOMJICHHA TakOH Kak H3BeCTHO, B OOacTH 
HM3KHX NOpsAaKa MOHHOM 
YaCTOTbI NOKa3aTeJIb MpesOMJICHHA M1a3Mbl paBeH 
IIpHMepHO OTHOWICHHIO CKOPOCTH cBeTa K aJibdpBe- 
HOBCKOH ckopocTH N~c/c,4. B H3BeCT- 
HbIX Cly¥aeB JIA Nla3Mbl, B MarHHT- 
HOM <c, Tak N> 1. Ho yxe MHTeHCHB- 
HOCTb JHMOJbHOTO B Cpese C MOKa- 
3aTeseM N BN pa3 HHTCHCHB- 
HOCTH B BakyyMe ([2], cTp. 367). 


BOJIHbI TO MHTCHCHBHOCTb MOXKeT 
OKa3aTbca TOpa3i0 HHTCHCHBHOCTH JHMOJIb- 
pacueTOB HOHbI MOTyT H3Jly4aTb He 4eM 
QICKTPOHbI, C TOM CKOPOCTbWO B 
BakyyMe. BeMYHHa NOKa3aTeNA 
N CKa3bIBaeTCA Ha CneKTpalib- 
TeOPpHH CHHXPOTpOHHOTO H31y4eHHA [3] H3BeECTHO, 
4TO TIpH CKOPOCTH 3apsaa K CKOpOCTH 
cBeTa MaKCHMYM MHTCHCHMBHOCTH H3J1y4eHHA CMeLla- 
eTCA C UMKIOTPOHHOH YacTOTbE Mpe= Bi(mec), 
Ha ee OGepTOHb] NpakTH4eCKH 
ITOT 3aBHCHT OT OTHOLUCHHA CKOPOCTH 
K CKOPOCTH BOJIH H He 
CBA3aH C peJIATHBH3MOM. N> |, koraa 
ckopocTb V=—c/N<c, cmeuienne 
MaKCHMYMa HHTCHCHBHOCTH axe 
CKOpOCTM BeCbMa OT penATH- 
BUCTCKOH v<c. KoneyHo, MaKCH- 
MYMa MHTCHCHBHOCTH MarHHTHOrO HOHOB 
He MOxeT ObITb TaKHM CHJIBHbIM, Kak CHHXpo- 
TPOHHOM 3IeKTpOHa, TaK Kak 
yactoT, rue N> 1, B naa3me orpannyena. [lpencta- 
BIAKOT HHTepec NOITOMY YaCcTOTbI =m (m= I, 
2, 3...) CO CpaBHHTeNbHO HOMEpOM 

3apalla B MarHHTHOM Movie B Cc 
noka3aTesiemM mpenomienua 
Llbitopuyem [4], KoTopbiii, OfHakO, 
cpeny w30TponHoi. B mla3mMe 


* B w30TponHOK cpene c mOKa3aTenem nmpenomsenun N 
OKPyXHOCTH, paBHa 


2 e? op? v2? N 


3 [1 — (vw N/c)?}* 
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B. MUCTYHOBHY, B. Jl. LIA®PAHOB: 


nocBaueHbt [5] u TsBucca u 
Po6eprtca [6]. Onnako, ux 
TOJbKO K W31y4eHHA B 
MpHBOAMTCA KpaTKO MeTO 
TeH30pa 3aTeM (pop- 
MYJIbI TIPHMCHAFOTCA K 
ABIAIOTCA MOHbI > me). 
B nHactosuei paSote 4TO 9TH 
MOHbI JBHKYTCA B XOJOAHOM Ma3Me. ITO O3HAaYAaeT, 
4TO HX CKOPOCTb 3Ha4HTeJIbHO MpeBbiiiaeT 
CKOPOCTb MOHOB M1a3MbI > vp. yes0- 
4eHHA B CTOKHOBH- 
TOYHO Mavia. 


2. Meron pacueta 


PaccMOTpHM cpesy B KOTO- 
poi ABHXKYTCA 3apAAbI C ONHOPOAHOM 
B Cpesle TOK C 
j none E, paSote Hag 3apanamnu j-E 
(pa6ota cuJibI TpeHHs), B3ATOM C OOpaTHbIM 
3HakOM. Metro pacueTta mona 
noTepb B H310%KeH B 
pa6otax [un36ypra [7] u op. [8], [9]. Ana yao6ctBa 
ero KpaTKO B H3MeHCH- 


HOH dopme. 
BCce E, j ap. no naockuM 
BOJIHaM 
E(r,t) = dkde. (1) 


Aina amnoutya Ex, o, jx, 43 ypaBHenuit 
COOTHOWICHHA 


Peuienve ypaBHennii (2) B 


I=1 


rae a — coOcTBeHHble BeKTOPbI NOJIA, 
CHCTeMe ypaBHeHHii 


{NP (dx = 0. (4) 


T.K. Mbi mpeHeOperaem TO TeH30P 
IPMHTOB M3 (4) TakxKe ypaBHeHHe 


{ Nu? (bag — Nang) —Epo} Ona = 0. (4a) 
Ycnosue B HYJIb JeTepMHHaTa CHCTeMbI 
(4) (wim 4a), onpenenseT, KaK H3BECTHO, 3Ha4e- 
HHA KBalpata noKa3aTesA nmpenomsenua (/— 1, 2) 
KOTOPbIM COOTBETCTBYHOT COOCTBEHHBIX BeKTOpa 
a). YMuHoxas (4) Ha ays, a (4a) Ha aig 
; 3aTem (4a) u3 (4), 


Ng) = 0. (5) 
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[pu m, MHoxutTenb nepea N? — Nn? oOpautaetcs 
B HYJIb, /=m ero paBHbIM 
Takum o6pa30M, MMeeM yCIOBHe OpTo- 
HOPMHPOBKH BeKTOPOB a1: 


Ay, p Np) = Om (6) 
= NP Om. (6a) 


Bropoe u3 ecTb (6) (4). 
Tenepb (3) B (2) Hcnonb3yA ycuOBKe 
OPTOHOpMHpoBKH, Ey ,, B Bue 


2 


= Lp jp» Lap = o “NP (7) 


Bbipaxkenne NOTepb OTHECeH- 
HbIX K OObema O = —j (Fr, nony- 
yaeTCA B B KOTOpbiii BXOAMT KOppe- 
(pypbe-KOMMOHeHT MJIOTHOCTH 
TOKa 


G_g NAOTHOCTH paBHa 
[10] 

Gop = 0, Gig, 


co 


3necb V(t) — 3aKOH 3apata, yrsoBbie 
CKOOKH O3Ha4alOT yCpeaHeHHe MO Ha4aJIbHbIM CKO- 
poctam v (0). 

Bce OJMHAKOBYHO CKOPOCTh, 
TO yruloBbie CKOOKH HyXKHO OnycTHTb. PaccmMoTpum 
9TOT 

(8) B BbIpaxeHue O = —j(r,t) E(r,t), 
HaXOJMM NOTepb OAHHM 3apsOM: 


dk 


Q 
Lap Akd (10) 


3necb MHTerpHpoBaHHA GeckoneyHbie. Yun- 
TbIBaA, MO (8) Gig = a Takxe, 
4uTo corsmacHo (9) Gig (—k, = (kK, m) ana- 
€ap (—k, = MoxHO 3anucaTb (10) B 
Bue MHTerpasa MO 4YaCTOTaM 


co 
(11) 
0 
BeuiecTBeHHadA 4aCTb mpi 
BbIpaxeHHA, TOUHEe B (T.K. Gz, He HMeeT 
(kc/w)? = B THX TOUKAX COOTBETCTBYHOT 
6) npu Hezz=0 (Ocb Zz Cc HampaBuie- 
k), BbIveTbI B 39THX COOTBETCTBYIOT 
H3J1YYCHHHO BOJIH, YTO TO CaMoe, 


Bbl 
Or 
HH 
Md 
Me 
BI 
; 


eTCA 
Hue, 
DTO- 


(6) 
(6a) 


(7) 


CeH- 
ppe- 


(8) 


ITO- 
iBHa 


(9) 


BbIe 
CKO- 


CTb, 


r,t), 


(10) 


OKe, 
aHa- 
rBue 
)) B 


(11) 


npu 
AceT 
npu 
BJIe- 
Moe, 


3leCb He MHTepecyioT. 

3aMeTHM, 4YTO BMeCCTO B3ATHA BbIYCTOB, MOXKHO 
BbLICIUTb IPMHTOBY YacTb TeH30pa 
Torna 3anucaTb (11) B Bune 


oo 


Lip = (Lap + Lia) 
0 


(lla) 


pacCMOTpeHveM NoTepb Ha H3J1y4e- 
4TO Cpena GeckoHe4HO 
MaJIbIM norsouenHem (Im N? >0). Torna 3nHa- 
MeHaTesib B MMeeT BU 
rac P — 
Takum o6pa30m, 


4x? 

I=1 


k2e2 
* 
Aly 


onepatop 3HayeHHe™, d(x) 


va) (12) 


U3 Bbipaxennit (11a) (12) cnenyet, NockoNbKy 
(kc/w)?>0, w3ny4¥eHHe NpoucxonHT B Tex 
yacTOT yrsioB, rae N?>0. 

B BLIpPaxenue (11a) TeH30p Gs, 
BbIMMCJIMTb B CHCTeMe KOOPHHAaT, NOCTOAHHOe 
MarHMTHOe Nose By HanpaBseHO NO OCH a BEKTOP 
k B MOCKOCTH Xp» Zo. Onyckaa 
pe3yJIbTaT 


co 
Cha = — key os) (13) 
m = co 


xoxo 


= — IND, = 
TT = (A) 


118) = (A), 


Xoo 


m2(A). (14) 
HO NOCTOAHHOrO MarHHTHOTO MOJIA) KOMMOHECHTbI 
ckopocTH k,, 
KOMIMOHEHTLI BOJHOBOTO BeHTOpa k, 

ki v e By J/ 

OB Mee —P. ( ) 
B Toi CHcTeMe KOOpAHHAaT BeKTOp a HMeeT BH 
(cos =k,/k): 


a = dy{ia,,,1,ia,,} 
1 


1+ 


= = cosf—a,, sin . (16) 


LIUKJIOTPOHHOE 


KOMIOHECHTHI 


— = 


= & =@, (16a) 


BCJIM4HHbI %:, H N* 


N?—e N?gsin@cos9 (17) 


2 (e? — g?) 
2 72 
}sint | : ry sin’ @ +49? cos? 


(18) 


N? = 


B (lla) caeprky a, ay*, 
HaXOJHM 

mv 


(19) 


o6pa3HO CyYMMHpOBaHHH NO OTPHUaTebHbIM 
3aMeHHTb m Ha B (11a) HHTerpupo- 
BaHve no k c NOMOLUbIO BXOAALUIMX B 
Lug’ (12), BbIpaxkeHve AIA NOTepb Ha 
yeHHe B 


1 


2 co 
ew? 
0 


+ 
(20) 


mv, 


xd —m op, Nip) 
mv, 
= a= — v)) Ja 


Al 
x + mog— Nip). 


3aecb apryMeHToM dyHkunit Beccena 
etca A) =(w/mp) (v,/c) 1 — Aprymentsi 6- 
NoKa3bIBaroT, 4TO /, ONpeneNAeT YepeHKOB- 
ckoe 3apsta, — WHKIOTPOHHOe H3u1y- 
yeHHe C y4eTOM HOpMaJIbHOrO 
mp>0), C yYeTOM aHOMAaLIb- 
Horo [11]. Mtr 6ynem B 
HHTepecoBaTbCA TOJIbKO Ha pe3o- 
HaHCHbIX “| =m Mp, KOTOPOe ONpeneAeTCA 
cularaeMbIM Im. 


3. HOHOB B XOJ10QHO N.1a3Me Ha 
HbIX 4acTOTaXx 


o~m opi (m=1,2, ...) 3aBHCHT OT 3Ha4eHHA Napa- 
MeTpa A = @oj/Mpi, PaBHOTO OTHOLICHHHO 
POBCKOH HOHHOHM YaCTOTbI K HOHHOH 
YaCTOTE, HJIM, TO Ke CAMOE, OTHOLUICHHIO CKOPOCTH 


19] 


| 
| 


B. HW. B. Jl. LIA®PAHOB: 


_cBeTa K ckopoctTu. 
la3Mbl 3TOT MapaMeTp paBeH 


By 


npegnonarath, A> [pu B,—3-10* 
raycc (O6bI4HO Nose B Orpe) yc10- 
BHEe BbINOJHACTCA > 5-108 cm~*. 3ameTuM, 4TO 
Ha Orpe B Tex pexkHMax, re CyLueCTBeHHO 
FOTCA Pe3OHAHCHbIe WieKTpH4e- 
ckoro noma Ha oOGeproHax (m~3-4), NAOTHOCTE 
Ma3Mbl MpeBOCXOHT — 10° 
A> 1 3aBeOMO BbINOHAeTCA 
JOCTHXKUMBIX 3Ha4eHHAX By. 

ycnosun A> g, 7, ONpeneNsA- 
MOKA3aTeJIb MpesIOMJICHHA BEKTOP MOJIApH3a- 

A? 


_ __ mi A? 
y= 


&= 
Me x? 


(22) 
(mj, Me — MaCCbI MOHA H WICKTPOHA, COOTBETCTBEHHO). 

OTHOWIeHHIO K By) 21eKTpH4ecKoe Mose 
Ez,=D:z,/yn Mano (OHO 
QWICEKTPOHOB BAOJIb MarHHTHbIX CHJIOBbIX JIMHHi). 
Kak BugHo (17) (22) 1. Ony- 
CTHM NO3ITOMY CylaraemMoe az, B (20). 
W31y4eHHA Pm Ha Pe3OHAHCHbIX YacTOTAX (CJIaraeMoe 


e o Bi 


Pu = v A Pm (23) 
Im (A) 4 In’ (ay 
Pm = | du, (24) 
(25) 
2 
= = ) 1 (26) 


3necb xn V1 — BbIpaxke- 
Hue B (24) cnenyeT OpaTb 3HayeHHH x, 
BOPAIOLUICM (HyJIb apryMeHta 0-(pyHK- 
B Jn) 


m 


(28) 


B cbopmy.iax (24—28) BBeneHbl 


c CA 


Jina Nonyyaem u3 (18) c yseTom (22) 


9 
=- 


m 

(30) 
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3aBHCHMOCTb KBapaTa 
OT YacTOTbI H306paxeHa Ha pue.1. Kak 


! 


x, 1 2 3 4 Wp, 


Puc. 1 


n,?>0 Bo paccmMaTpuBaemoii 4YacToT, 
Torga Kak n*>O x< x,, rae 
it 
(1 — p*) 
m 


i 


(31) 


B BbipaxkeHHH cnharaemble C MOXKHO 
onycTutTb. Jia 9TO MOXHO IA 
yrsloB BHe y3Koro KOHyca (m-/mi)"/*. 
Hoe BbipaxkeHHe umeet Buz [12] 


2 


= = 
1+ 4 42? 2 


(32) 


PaccMOTpuM CHa4aJia BOJIH, XapaKTepH3y- 
una F B 9TOM paBHa 


V/1 + 2(2a?—1) 
pm Oblia YHCCHHbIM HHTe- 

MpH 3HAYCHHAX MapaMeTpoB 

Bz B Tabsuue I. 

MaKCHMYM HMHTCHCHBHOCTH CMeLlaeTCA B 

cropouy m*. 

OG6pauiaeT Ha BHHMaHHe TOT (akT, 4TO NpH 
MaJIbIX ¥,/CA WHTCHCMHBHOCTb Ha 
WHKJIOTPOHHOH YaCTOTe 3HaYHTeIbHO 4eM 
MHTCHCHBHOCTb Ha COCeMHUX OGepToHaX. 
BbIACHHTh 3ITOTO pa3zJIM4HA PaccMo- 
A< 1. U3 (28) x =m, Tak 4TO 


* 3ameTuM, 4TO Kak (23), HHTCHCHBHOCTE H3J1Y- 
4CHHA HOHA MOXKET CPaBHHBaTbCA C MHTCHCHBHOCTbIO 
WICKTPpOHA B BakyyMe NPH Apm~(mi/me)’. 
Kak BuaHo 43 6y~1-1,5 m> 20, pm mocTuraet 
> Tak npu A> 10° cymMapHas MHTCHCHBHOCTb 
HOHA MO*KET HHTCHCHBHOCTH 
CKOPOCTbIO B BaKyyMe. 


lo 
1 
1 
| 
| 
| 
4 
| 
H 
E 
3 
= 
| 
ae 


HHA 
HO, 


OT, 


LIUKJIOTPOHHOE HOHOB 


TABJIMLIA I. pm- 


7 


{ 


¥ : 1 sis | 4 6 | 8 10 | 12 wt 16 | 18 | 20 | 22 | 24 | 26 | 28 | 30 32 | 34 | nl 38| 40| 42 | 
| By | 2 | | | | | } | 
| 1,4-10-4 | 4.3. 10-3 | 3,1-10-4 | 
70.5 [0,17 0,14) 5,8-10-* | 1,7-10- ry 4 
0.5 | 71-109 /013 O18 0,15) 6-10-? | 1,7-10-8 
0.7 | 0,07/2,8-107|021 1055 108 1095 | =| | | | 
0,7 | 2,3-10-* | 0,29 0.65 |09 | 1,03 10,86 | 0,62) 0,41 | 0,26 | 
0,8 | 0,08| 4,6-10-3 0,28 0,84 15 | 2.5 3.1 (3.2/3.1 | 24 | 2.0 | 1,67] 1,33) 1,03] 0,78 | 0,6 
0.8 | 3.8-10-? 0.4 1,04 1,7 |28 (3,3 3.4/3.2 |29 | 17 | 1,36 
[009/710 [035 (12 (24 [87 [12 [15 [22.9] 23. 22.8) 22 [21 | 20 
10/01 116 135 196 [30 | 103 | 128 [156 |187 |220| 255/208) | | 
1,0 | 9.6.10 0,75 2.2 43 | 11 | 20 32/46 |64 | 131 | #58 
13 |66)20 65 | 88 108 | 125 | 143 | 168 | 208 | 267 | 345 | 436 | 527 | 608 | 66x | 707 755, 797 
| 1.5 10,87 131 [68 |16 jee j29 |75 | 105 | 122 [131 | 152 | 202 | 265 |309 | 322 | 333 
(14 B,un <1. Monaras, cornacno (28), x=1+f, un-, 
Pm | (34) nomyyaem (30) np %KCHHO 
1 
n?(1—m?)—1 n_? = — . (36) 
m ‘ 28, n_ 4 
mj 
n? = (35) > (elm) (2 = we? 
1+ V/1 +2 (2m? — 1) + pu (me/mi ( Bs Ho 
*KUTb = — (2 n_ T.e. n_ By = B,?/2 
4TO ycnoBHe f, <1 < 1). Ucnon- 
3yYA 9TO 3HAYCHHE N_, 
m= 1, Kak BUQHO, «=—1 


p,—0. Takum o6pa30m, B AHNOJbHOM 
UMKOTPOHHOH YacTOTe OTCyTCTBYyeT. B 
YeHHA WICKTPOHOB pe3yNbTaT Obi NoNy- 
yeH B paOotax [5], [6]. Ycnosue «=— (axg=— 1) 
O3HadaeT, YTO NpOeKUMA BeKTOpa Ha 
ILIOCKOCTb, K NOCTOAHHOMY Mar- 
HHTHOMY MOJO, BpallaeTcA NO Kpyry B 
HallpaBJICHHH, BpallleHHO 
ABIACTCA TIPHYHHOH OTCYTCTBHA 
3aMeTMM, 4TO 3HAaYeHHAX Mm B MHTe- 
rpasie (24) ocHOBHOH MaJibie 3Ha4eHHA 
CKOpOCTb 6M3Ka K 
CkopocTH V~c4. Dopmyna WIA NoKa- 
3aTeA (32) WIA Ya- 
CTOT, He HWKHIOW 4a- 


ctoty [13] @ mj/Me 


‘mj/me+ A? 
m He OJDKHbI MeCHbUIyIO 

O6OpatTumca Tenepb K BOJIHaM C Apyroi nonsapH3a- 
uneit. 

Bonubl mpenomsenua N_=An_ 
MOryT Npw x<x,< 1 (31). 
(aBwxKeHHe NO OKpyxHOCTH), cormacHo (28) 
x= 1(@= pi) OTCYT- 
crsyeT. IIpu (apuxeHHe NO BHHTOBOM 
pe30HaHCcHas YacTOTa CMellaeTca <0 Bcnen- 
B cTopoHy HH3KHX 


*, TaK 4YTO 3Ha4eHHA 


Vo 


MHOXUTeNIO mi/mMe MHTCHCHBHOCTb H3- 
OKa3biBaeTcA (axe /,) 
oueHb DTOT pe3ybTaT ABIAeCTCA, OMHAKO, 
BOJALIMX K BO3MOXHOCTH Kak yrOHO GOJIbLIMX 3Ha- 
yeHnit N_. HetpyaHo 4TO BCe H31y4eHHe NpH- 
XOAHMTCA Ha 4YaCTOTbI, BeCcbMa K 4YaCcTOTe 
© = Mpi X, KoTOpoH N_*= co. 3Hayenne B 
YACTOT OKa3bIBaeTCA NOITOMY 
BecbMa N_*~ A? (mj/me). 
ONMCaHHH peasibHOH Ma3Mbl Yy4HTbIBATb 
(bakTopbi, K 3Ha4eHHii 
N_*. Takumu MOryT ABHTbCA 
ABUXKCHHE 3apAOB M1a3Mbl, HX CTOJIKHOBe- 
HHA, MarHHTHOTO Nowa. Orcona 
aHHbIX MO H3MepeHHlO pe3OHaHCHbIX MMKOB Halips- 
X*KCHHOCTH C MHTCHCHBHOCTbIO 
PAaCCYHTAHHOH B paGoTe, 
myeT M3 PaCCMOTPeHHA 
yacToTy. 

y4HTbIBaTb TONbKO OGepTOHbI (@ = m 
m>2), TO cMelleHHe MaKCHMyMa HHTCHCHBHOCTH 
MPH yBesMYeHHH NapaMetpa /,=(v,/c) A 
HaXOJHTCA KaYeCTBCHHOM CO 
[1] pe30HaHCHbIX MHKOB HalpsxKeH- 
HOCTH NOJA MIOT- 


HOCTH 11a3Mbi (HaNOMHUM, 4TO A ~ V 
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31) 
<HO 
eH- 
| 3 mj ‘ 
| 3 (37) 
32) mo 
3y- 
HK- 
33) 
iTe- 
0B 
el. 
[ca 
A B 
Ipu 
Ha 
ax. 
MO- 
1TO 
le)”. 
\CTb 
CTH 
*Ke 


B. U. NMCTYHOBHY, B. Jl. LIA®PAHOB: 


4. OGcyxjenne pe3yLTaTOB 


HHTCHCHBHOCTb H3J1y4eHHA ObICTPbIX HOHOB B XOJOJ- 
HOH pacTeT C yBeJIM4eHHeM MJIOTHOCTH 
XapakTep HHTCHCHBHOCTH 
M3J1Y4CHHA MO TAPMOHHKaM Takxe CYLICCTBCHHO 3aBH- 
CHT OT MJIOTHOCTH yBeMYeHHH MIOT- 
HOCTH MIpOHCXOAMT CABAL MaKCHMyMa HHTCHCHBHOCTH 
B CTOPOHY BbICOKHX ObepToHOB. 93TH 
MO%KHO HCIOJIb30BaTb 
B MarHOCTH4eCKHX LeIAX 
HOCTH [J1a3MbI. 


Cnenyert, OAHaAKO, HMCTb BBH JY, 
(bopmy.1bI HHTCHCHBHOCTb H3JIY4YCHHA Ha 
PaccTOAHHAX Rno CpaBHeHHbO C 
2. B a6opaTopHbix yCSIOBHAX TpeGo- 
R >A BbINOJIHACTCA (ecm aHTCHHa 
B mJ1a3Me) CpaBHHTeJIbHO NAOTHOCTH 
MJla3Mbl, HMCHHO 

e? Nj 


> 
mj c* 
(Ni=ny:2 — ,,mOroHHOe HOHOB™ B cTonGe 


nla3Mbl a, m — HOMep 
TrapMOHHKH). 
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B 3akHo4eHHe aBTOPbI BbIpaxKatoT G6.1a- 
ronapHoctb A. E. BAKAHOBOM 3a 4H- 
CJICHHbIX Ppac4eTos. 
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SIEKTPOHHOKM TEMMEPATYPbI NIA3MbI 


B MOWHONW BOJHE 


®ununnoBa T.U., H.B., B.B., BUHOrPAOB B. 
OPJIEHA JIEHAHA UHCTUTYT ATOMHOM UM. HU. B. KYP4UATOBA 
AKAJIEMUUW HAYK CCCP, MOCKBA, CCCP 


MeToO0M MarHUTHOrO NOTOKa H3MepeHa IIEKTPOHHAA 
3a yoapHoi BonHOK. ckopoctsax oT 0,9-10" no 1,25-107 cm/cek Halinena 
TpOHHas Temnepatypa B 50 u 90 eV. IIpopeneHo NONyYeHHOH C pacc4n- 
TaHHOH MO CKOpOCTH BOHbI. H3MePeHHA CHKATHE CHJILHbIM MarHuT- 


Baenenne 


TEPMOAEPHbIX peakuHii, NOKA3aIM, YTO 
UMIYJIbCHbIM MOXKET HCTO“YHHKOM 
CHJIbHBIX YapHbIX BOJIH [1]. 

H3 BapHaHTOB CO32aHHA CHJIBHbIX 
BOJIH ABJIACTCA HMMYJIbCHbIM B 
yeckoH Kamepe. Ilo KamMepbl O6pa3yetca 
MeHHbI WHYp, B KOTOPpOM Mla3Ma Nos 
OObIIMM aBeHHeM H3-3a Ecuu B 
KaMepbl OTBepcTHe MO OCH H B 
oTBepcTHe NOCTaBHTb TpyOky, TO B 
BaTb KaK KOTOPpbIM oOOpa3yeTca 
ylapHad BOJIHa, no TpyOKe c 
CKOpocTbWO. NapaMeTpbi 11a3MbI 
3a BOJIHOM, MOXHO HeKOTOpbie 
BbIBOJIbI OTHOCHTEJIBHO COCTOAHHA B 
Kamepe. 

B [1, 2, 3, 4] u3y4anocb ra3a, HarpeToro 
C NOMOLIbIO pa3paaa. 
OHako, BO BCeX yKa3aHHbIx paGoTax 
HeOCTaTKOM ABJIACTCA MpHCyTCTBHe B 
xyuleHCA 3a yapHOH BOJIHOM, MarHHTHOTO Moss, 
CBA3AHHOTO C OCHOBHbIM 

B MarHuTHOe Nose OCHOBHOTO pa3- 
pala He B TpyOKy, no KOTOpoOii pacnpo- 
cTpaHsslacb Ha TpyOky HaMaTbIBaJINCh 
BUTKH, NpOK3BOAHYIO MarHuTHOTO 
noua. 
4TO MMMyJIbc MarHHTHOrO 
NOJIA OYCHb Masi GbICTpO 3aTyXaeT NO Mepe pacnipo- 
CTpaHeHHaA BouHbI (3a 0,5 usec). 
Hanmune takoro cnaG6oro 
KOMMOHCHTbI MOXKHO OObACHHTb Kak pe3yJIbTaT 
Hekoero HepaBHOBeCHOrO Mpouecca, 
3a BONHOK cpa3y nocne oO6pa30BaHHA 
HOH BOJIHBI. 

OtMetTuM, 4TO B paGoTte [5] MarHHTHOe Nose, cBA- 
34HHOe C OCHOBHbIM TOKOM HMeeT 
pacnpocTpaHseTca BMecTe C NOTOKOM 
njla3Mbl, 3a (PPOHTOM BOJIbHBI. 
B onucbipaemoi paboTe reoMeTpHaA 9KCHepHMeHTa 
MIPOHHKHOBCHHE OCHOBHOTO B TpyOKy. 


1. TeMnepatypbi 7, 
3a yapHoli MeTOOM BbLITeCHeHHOrO MarHHT- 
Horo HoTOKa 


1.1. 


Tlna3ma, 3a BONHOH B NocTO- 
AHHOM HarHOCTHYeCKOM MarHHTHOM 
nose, AHAJIOTH4HA MPOBOHHKY C o, 
Temmepatypoii 

3aa4y O MpOHHKHOBeCHHH MarHHTHOTO NOMA 
BHYTPb HalieM ee MpOBOAHMOCTb B 
4TO: 

1. o=const. no ceyeHHtO MMeeT 
(bopmMy HeCxKMMaMeMOrO Apa, 
TOporo MeHbuie, HO MOxeT ObLITb CpaBHHM Cc 
TOJILUMHOH 

2. ¢=const. no BpemMeHH. 

3. Bpema BuleTa B MarHHTHOe Mose MHOTO 

MeHbIIIe, BPeMA 10/18 BHYTPb 3MbI. 

B CHCTeMY KOOPJIMHAaT, JBHXKYLILYIOCA 
CO CKOPOCTbWO), paBHO CKOPOCTH Vz 3a 
HOH BONHOM, TO CBOAMTCA K B 
MOMECHT BPeMeHH fy MpOAOMbHOe Mose 
Ho, KOTOpoe NPOHHKaeT BHYTPb M11a3- 
MbI. ypaBHeHuit: 


4 1 cH 
rot E =—=5, (1) 
4x, ‘ 
rotH = (2) 
E = ae (3) 
culeayeT, 4TO 


Baenem F(r,t) Takyto, 4TO 

H (r,t) = F(r,t) + Hy (5) 
rae H, — none, B kOTOpoe 

(4) 

oF 4-0 OF 

Gr) = (6) 
ypaBHeHHe TenAONpOBOAHOCTH C IpaHH4HbIM 
: 
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4H- 
KO, 
all, 
335, 
il 
ce, 
pa- 
CP, 
31 
4 
AH =" (4) 
Zz 


F =0 (7) 
rae ry — m1a3Mbl C Ha4aJIbHbIM YCJIOBHeM: 
F(r,0) = —H, (8) 

MMeeT BU: 


n=1 
THe fn — KOPHH ypaBHeHHa Beccesa: 
(tun) = 0 (10) 


Orctoa MOTOK MarHuTHOrO Mods, 
Cro HallpAKCHHOCTh MO 


co 

exp 

n 
Ss n=1 (11) 
ITO *e PaBeHCTBO MOXHO 

o6pa30m: 

=Po—% (12) 
rae Po = (13) 
MarHHTHbIi MOTOK Yepe3 ceYeHHe, paBHOe 
JO ee BllycKa B Nose Hy-g, — pa3HHlla Mex Ly 
M Po, T.€. TOT MOTOK, KOTOPbIii BbITeCHACTCA 
BBeeHHH B Hero 3a 


Bpema f. Tellepb MOXHO ceMeiicTBO KpH- 
BbIX 


n=1 


Kak (pyHKUHMIO BpeMeHH, MeHAA or,?. 
KpuBbie ry=0,55 cM no- 
CTpoeHbI Ha 1. 3Had MOKHO 


06 


' 6= 3x10" 
6= 2x10" 


10" 


04 


= o=05x10" 
0.0 
0.2 04 06 
Puc. 1 3apucumoctb OT BPeMeHH MpOHHKHOBeHHA 
MarHUTHOrO MOUIA. 
B kayecTBe nmapaMeTpa B3sTO ory”, — 


Yo — U3mepenna B 
MOMEHTHI fy OT BXOXKICHHA BOJIHbI 
B OOsacTb MarHHTHbIM TIOJIEM. 
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T. OUSIANITIOBA, H. B. B. B. KYPUH, B. Tl. BAHOPPAJIOB: 


slerkKO MpOBOAHMOCTb O KOHKpeTHOro 
OtMeTuM, 4TO B IKCHEPpHMeHTe Ope tesia- 
JIOCb H3 JIOMOJIHHTCJIBHbIX u3MepeHHit 4eTbIpemM 
KOHTPOJIBHbIM KaTYLIKaM. 


1.2. DKCMEPMMEHT 


Ka- 
Mepa (puc. 2). Bataped B 85 


© 


Puc. 2 Cxema yCTaHOBKH. 

— KamMepa, auamMetTp 40cm paccTOsHue 
45 cm. 

2 — cTeknsHHasa Tpy6ka, 

3 — kaTylKa 18 CO31aHHA MarHUTHOrO Ho, 

4, 5 — u3MepuTesbHbIe BUTKH, 

6, 7 — 


20 30 kV, ckopocTb HapacTaHua ToKa 
4-10" A/cex, paGounii ra3 Hayanbubie 
0,2 u 0,05 mM pr. cr. 

B paOote [6] CHIbHO HOHH30BaHHOrO 
aproHa, CO3ZaHHOrO B TpyGe, 
C MNOMOLIbIO MarHHTHOH KaTyLIKH. 
KasmOposka curHasioB 
BaHHeM Yepe3 MarHHTHY!W KaTYLIKy 
OOpa3llOB H3BECTHOM MpOBOAMMOCTH. 
OTHOCHTCJIBHO HH3KHX NapaMeTpoB Kora 
BpemA MarHHTHOTO MOA Maso cpaBHe- 
HHO CO BpeMeHEM B MarHUTHOE Mole. 

B onvcpipaeMbix 
MarHUTHOe MOJe CO3qaeTCA KaTYLUKOM, HAaMOTAHHO 
MeHOH 2 x | MM B 9 cnoes. Ha Konuax ee 
(POpMbI MarHHTHOTO Nona, K 
HoH. C WeNbIO Ha KOHUAX 
TpyOouky, KOTOpad 9IKpaHOM, MpensATcTBy- 


Tec 
oT 
Ip 
no 
Me 
| 
( 
4 
4 M 
(7) 
/ M 
K 
M 
9 
— 
| | 
| 
oy 


TO 


Ka- 
uF 


M3MEPEHUME 3JIEKTPOHHOM TEMIEPATYPbI MJIA3MbI 


MarHHTHOTO MOTOKa, BbI- 
TecHaeMOroO MarHuTHad KaTyuika 
oT GaTapew akkKyMyJIATOpoB, 12, 24 u 48 V. 
9TOM BHYTPH PpaBHAJHCh COOTBET- 
cTBeHHO 1150, 2300 4600 oe. 
W3MepsJIcA BHTKaMH, 
yeckOe MarHHTHOe Nose. B 3aBHCHMOCTH OT 
MeHTa BAHTKOB M3MeHsIOCb. CurHasl C BUTKOB 


O5us 


| 
| 
| 


Puc.3 OcunnnorpaMMa CHrHasIOB C 
BHTKOB 4, 5. 


(puc. 3) perucTpuposasica Ha 
com ocuMiorpade OK-21 nocie 

Ilo BbITeCHeCHHOTO TOTOKa Ha- 
XOAMJIACb CKOPOCTb ylapHoi BOJIHbI NO BpeMeHH 
ee JBYMA BHT- 
kaMu. Oua B npenenax oT 0,9-10* cm/cek 
1,25- 10° cm/cek. 3ameyeHO, 4TO yAapHoi 
BOJIHbI Ha TpyOKu (/= 50cm) 
Maso. [losTOMy pacueToB ee 
MOX%KHO MCKAY ABYMA BHT- 
KaMM TIOCTOAHHOM paBHoi 107 cmM/cek. 

PacyeT MpOBOAHMOCTH He TOMbKO B 
ca60M NOCTOAHHOM, HO H CHJIBHOM JMarHOCTH4eCKOM 
MarHHTHOM H 3Hay¥eHHe o =(5+2)-10" 
en. CGSE. MpoBOAMMOCTH 
WICKTPOHHaA TeMMepatypa, [7], 
no (bopmy.te: 

(cs 

T. — Temnepatypa B eV, 
T.= 50 eV ana o 5-10" en. CGSE. Makcumasbuoe 
3Ha4eHHe eg. CGSE, co- 
ei TeMMepaTypa 
T.~ 90 eV. 


2. 


K yAapHoro PpoHTa 3aKOHbI 
COXpaHeHHA Maccbl, IHEPrHH MOXKHO 
OWCHHTb TeMMepaTypy Mla3Mbl 3a 
3Had TOJIbKO MapaMeTp, HallpHMep, CKOpOCcTb 
TeMMepaTypa M1a3Mbl, ycs0- 
BHH PaBHOBeCHA 3a 
BOJHOH, He CHJIbBHO OT 3JICKTPOHHOM 
MOJyYeHHOH MeTOOM BbITeCHeHHOrO 
MarHHTHOrO MOTOKa paBHsAeTCA, mpHMepHo, 40 eV. 

pA 9KCNepHMeHTOB NO 
MarHHTHbIM 4TO HadHHaeT 
CKUMATbCA HallpsxKeHHOCTH MarHHTHOrO 
H=~10koe, cooTBeTcTByeT MarHutT- 


Horo nossa =4 Cue MOBaTeJIbHO, ra30KHHe- 
THYeCKOe BHYTPH Ma3Mbl, 3a 
BOJIHOH, ~ 4 aTM. M3 3aKOHa COXPpaHeHHA Ha 
BOJIHe HaxXOMHM p, 3a yiapHoi 
MO ee ckOpocTH 4, 


Ps = (1— (16) 

— OTHOWeHHe NIOTHOCTeH ra3a 3a 
HHA p, = 0,2 mM pt. cr. = 10’ cm/cek 
3Ha4eHHe Pp, OKa3bIBaeTCA PaBHbIM 3,5 aTM. 

CTaBHJINCh B TOpell CO cTeHKON TpyOKu 
Ha pa3JIM4HbIX PpacCTOAHHAX) p= 4aTM. 
Bce py, NONYYCHHHbIe 
OaMH, He CHJIbHO, a JexaT B 
TOYHOCTH 9KCNepHMeHTa. 

@MopmMa H BpeMsA 4), a Takxe 
pacueTbl, 4TO 


P atm 
5 


| 
Puc.4 MarHHTHOrO NOTOKa B 
TlopepxHoctb Ha 1 CM MIOCKOCTH BUTKa 5. 


PaBHOBECHe 3a BOHO! 
HacTynaeT BeCbMa ObicTpo (3a BpeMsaA NopsAKa 
10-7 cek) OTKPbIBalOT NepcneKTHBbI 
yeckoi ra30BOH CKOpocTeii. 


ABTOPpbI BbIPaxkaloT OKT. 
@PAHOBY B. JI. 3a none3Hbie OOcyxKaeHHA 
pa6otbi 6naronapat KOMMCCAPOBA B. B. u KO- 
JIECHUKOBA FO. A. 3a NOMOUIb B NpoBeAeHHH 9KCTIe- 
PHMeHTOB. 
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LOSS OF PARTICLES IN A PINCHED DISCHARGE IN AN AXIAL 


MAGNETIC FIELD* 


R. K. Jacer 


INSTITUTE FOR FLurp DyNAMics AND APPLIED MATHEMATICS 


UNIverSITY OF MARYLAND, COLLEGE ParK, MARYLAND, U.S.A. 


A calculation of the loss of deuterons from a pinched current to the wall of the container has 
been given by G. P. Thomson [Phil. Mag. 32 (1958) 886]. This calculation is extended so that 
account is taken of an axial magnetic field. It is found that such a field can materially reduce the 


particle loss. 


1. Introduction 


If a strong current is passed through a low density 
gas in a cylindrical container, the self attraction of the 
parallel currents carried by the electrons causes an 
electron concentration near the axis of the cylindrical 
discharge and a radial electric field is set up between 
the axis of the discharge and its containing surface. 
Such a configuration is, however, found to be unstable. 
The forces arising from the twisting of an axial current 
produce further bending of the discharge column and 
the constricted matter comes in contact with the wall 
of the vessel. An axial magnetic field is usually applied 
to stabilize such currents as occur in Zeta, Scepter, 
Stellerator and other fusion machines. However, it is 
found that the axial magnetic field stabilizes the pinch 
only for certain small wavelength disturbances [1, 2, 3]. 

THOMSON [4] studied the problem of the loss of heat 
from the pinch due to the collision of the particles with 
the wall. He found that the loss of heat is reduced by 
the binding of the paths of both nuclei and electrons 
by the toroidal magnetic field of the current, as well 
as by the electrostatic potential created by the con- 
centration of electrons near the axis of the cylindrical 
discharge. It is to be noted that the electric field, the 
particle density and the magnetic field (outside the 
pinched region) decrease towards the wall. 

The object of the present paper is to estimate how 
the loss of particles to the wall is affected by the 
presence of an imposed magnetic field parallel to the 
axis of the discharge. 


2. Equilibrium configuration 


Assuming that the velocity of the electrons is 
wholly axial in the steady state configuration and that 
the velocity of the deuterons is negligible compared to 
that of the electrons, we have the following equation 
for the radial equilibrium of the deuterons: 


ne By =k Ta , (1) 


and for the equilibrium of the electrons 


—ne E,— nev. H=kT. (2) 
Here n is the particle density of the electrons as well 
as deuterons, ve is the mean axial velocity of the elec- 
trons; H is the toroidal magnetic field in gauss; £, 
is the radial electric field in e.m.u. and is negative, 7'q 
and 7’, are the temperatures for deuterons and elec- 
trons, respectively. The presence of the axial magnetic 
field does not change the steady state described by 
Thomson. 

If Ta=T.=T, Eq. (1) can be written: 


(Inn), (la) 
which integrates to 


n =nyexp(—eV/kT), 


V= —[Bdr. (3) 
0 
Eq. (2) yields 
kT d 
ve H =E, (Inn) =—2E,. (4) 


Thus the result that the pinch force on the moving 
electrons is twice the radial electric field still holds 
when the axial magnetic field is present. 

The equilibrium radius of the pinch will be written 
ry. Following Thomson we shall assume that in the 
region outside the pinch V is of the form 2 a log (r/r,). 
i.e. the electrostatic potential due to a uniformly 
charged cylinder of radius ry, when the region outside 
the pinch is electrically neutral. Our calculations 
follow the same stages as those of Thomson. 


3. Equation of motion for a deuteron and its solution 


We calculate the chance that a deuteron of mass m 
and charge +e which has made a collision at r, 6, z 
(in cylindrical coordinates) will hit the wall. The 
equation of motion of the deuteron, in e.m.u., is 


d?r 
m az ~e(E+vxH). (5) 


* This research was supported by the United States Air Foree through the Air Foree Office of Scientific 


Research of the Air Research and Development Command under Contract AF 18 (600) 1315. 
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In our case this yields three equations 


m (¥ —r6?) = eR, + eH,rO + (6) 
= (r?9) =—eH, (7) 
mz=- —eHr. (8) 


Here H =(2 I/r) with J the total current strength in 
the pinch. H, is the uniform axial magnetic field 
applied externally. 
Eq. (8) integrates to 
where 2 is the initial axial velocity of the deuteron, 
r its initial radius and r’ its final radius. 
Eq. (7) integrates to 
: C eH 
where C is the constant of integration. 
At a point where r’ has its maximum we have 7 =0. 
The condition for the collision with the wall is 


r'max > R, (9) 


R being the radius of the wall. At such a point the 
equation of energy gives 


2 


2m 


= mo*®—2oeln + r) (10) 


where we have written + mq? for the initial kinetic 
energy of the particle in so far as it corresponds to 
its degrees of freedom in the axial plane. The equation 
can be written 


2Ie, 9 r’ 
m r m r 
2.2 
+= Ay 
4m? 


1 1 
Following Thomson we neglect C? (1/r?—1/r’), the 
error in neglecting this term being smaller than the 
error involved in taking m as corresponding to two or 
three degrees of freedom. With these assumptions 
Eq. (11) yields: 


1 


r’ r’ 2 
where 

a=2le/jm, b=4oelm, p=eH,/2m. 


The condition for the particle to hit the wall is that 

Eq. (11a) shall give a real value of z, for r’> R and 

r<R. We make the further assumption that b < p? R?; 

this condition is satisfied in most cases of interest. 
When r’= R, Eq. (lla) reduces to: 


1 
Z=aln + [w?— p?( R2—r?)]? , (12) 


PARTICLE LOSS IN PINCHED DISCHARGE 


where Z has been written for z). This will be real if 
1 


2 


All the values of Z given by Eq. (12) are not allowed, 
since we must have (see Thomson [4}): 


Zi<o. (14) 


4. Total number of deuterons reaching the wall 


When + log (r’/r) can be neglected in comparison 
with p? (r?—r*), Eq. (lla) can be written 


, 1 
Z=aln". + [w? — p? (r’? — r?)]*. (15) 


It will be seen that for r’> R the two values of Z lie 
between the values Z,, Z, of Z given by Eq. (12). If 
all the allowed values of Z are equally probable, the 
chance that a particle hits the wall is 


2,—2, 
(16) 


provided again Z,<m and |Z, <m. This can be 
shown to be true for most cases of interest, in particular 
if H, is large enough. 

Therefore the chance that the deuteron hits the 
wall is 


1 
2[o*—p*( 


1 
[1— 2, |? 


1 
= (17) 


If is the density of deuterons at any point and o,w 
the average value of the product of scattering cross- 
section and relative velocity of deuterons, the total 
number of collisions at the point is 


1 
> Wow. (18) 


Thus the total number of deuterons reaching the wall 
per second will be 


R 
N= [nto,w r2 2ardr. (19) 


We assume that is given by the Boltzmann law 


= noexp| In 4 = no| (20) 
where 
a=2oe/kT. 


We can also write Eq. (20) as 
» 
n=nr( (20a) 


Putting x=r/R and x =r/R, Eq. (19) becomes 


1 
- 
N =2no,w fx 1 x2 dx. (21) 
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R. K. JAGGI: 


If w/pR=2 mo/eH, is small compared with unity, 
r will be close to R and x will be close to unity. In 
this case the integral occurring in Eq. (21) will be of 
the order (w/pR)*, with a coefficient =1/3, if 
a=0 (ie. in the absence of a radial electrostatic 
field); the coefficient will be somewhat larger than 
1/3 if «>0. (Values of « can be estimated to be of the 
order of unity.) 

Thus one obtains 


- 2 
N = 2nBo,w nr? (22) 


5. Comparison of the heat loss to that calculated by 
Thomson 


We would like to compare our results to those 
obtained by G. P.THomson [4]. For rather weak 
axial currents the containment depends upon the 
electrostatic field and Thomson finds ([4], see cases 1, 2 
on page 891): 

(Nr), = R?/a = R?(k T/2e). (23) 


For strong axial currents, the electrostatic field is not 
important and could even be zero; for this case Thom- 
son finds (case 3 and Section 9) 


=270, w nk wla = 270, wn? R?(mo/2 Ie). 
(24) 


When a strong axial magnetic field is present our 
result, Eq. (22), gives 


Ny = 2x fo, wnr? R?(w/pR)? 
= npr? R?(2meo/eH, R)?. (22) 
We find 
Ni 
wor ~ (pr) er 
For Zeta, Thomson gives the data ([4] p. 894) 


I~2 x 10' e.m.u (2 x 105 amps) , 
T ~5-10® degrees K. 


Thomson’s figure o ~6 x 10° does not quite correspond 
to this temperature: from + mw?=kT (as given 
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by Thomson, p. 889) we find m~2 107 cm/sec. 
Thomson further mentions o=10!—10" (we shal] 
take o=3x10"); from which 
ro ~5.6 em; In hence we assume = 20 em. 
We take H, = 10‘ gauss. With m =3.3 « 10-*4 gm for 
a deuteron and e=1.6 x 10-28 e.m.u. we find «~1.4, 
p=24x10', (w/pR)~4.2 and N5j/(Ny),~ 
~2.4 10-3. 
Similarly 


Ny o \2 a Imo 
With the data given above: 
a~19 x 105, 


N; 
1.7 x 10-?. 25 
(25) 
This shows that an axial magnetic field of 104 gauss 
considerably reduces the loss of particles to the wall. 
We note that, with the data used here, 


b~58 x 104, 
p? R?~ 23 10", (26) 


so that the assumption b< p? R?, made in order to 
simplify Eq. (lla), is justified. 
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LETTER TO EpITOR 


THERMONUCLEAR REACTION RATES* 


Physicists working on controlled thermonuclear re- 
actions in the laboratory have had little occasion to 
consult tables of (ov) so far. Happily, the indications 
are that this situation will soon be changed. To pro- 
vide some measure of consistency, newly revised 
curves of (ov) for DD, DT and DHe® are offered, 
see Figure 1 on next page. The results given here 
supersede the writer’s ealier compilation [1], itself a 
revision of a still earlier one. 


The <ov) values are based on cross-sections which 
comprise in part the values given in Table I. The 
values of o below 120 keV are those of ARNOLD, 
PHILLIPS, SAWYER, STOVALL and Tuck [2] except 
that the neutron branch DD, has been modified to 
take into account the angular distribution of Pre- 
ston, SHAW and Youne [3] in accordance with the 
method of McNett [4]. The latter moves o (DD,)/ 
a(DD,)) upward at low energies so that the branching 
ratio now declines to unity at the lowest energies, 
instead of to 0.92. The values of o at energies above 
120 keV come from many sources tabulated else- 
where [5]. 


The Maxwell averages given here are in close agree- 
ment with values given at 5, 10, 25, 50, 100 and 
150 keV temperatures by HESSELBERG JENSEN, 
KoroED-HANSEN, SILLESEN and WANDEL [6]. 
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TABLE I Basic cross-section values 


Deuteron energy 
in laboratory 


in millibarns 


system 
(keV) DD, DD, DT DHe® 
15 0.065 0.065 15 
30 | 1.15 1.10 278 
60 6.8 6.5 2180 2.3 
120 23 20 4700 31 
250 46 38 1720 290 
500 74 59 660 630 
1000 96 78 280 230 
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THERMONUCLEAR REACTION RATES 


LO 50 100 400 
MAXWELL AVERAGES FOR 
THERMONUCLEAR 


3 
LOS ALAMOS SCIENTIFIC LAB. DHe 


JUNE 1961 J.L.T. / / 
DT 


REACTION RATES 


-3 
cm sec 


/ / 
= 
V J 
 taBLE FOR ESTIMATING 
NEUTRON YIELD FROM 
TOTAL DD YIELD 


| 
| V>DOn T keV 


| 0.508 2 
| bs 0.516 10 
0.535 20 
a 0.543 50 
10 0.550 100 
= 0.555 200 
R= nn, <ov> 
| <ov> 


Ke) 5 10 
T keV 


Fig. 1. Maxwell averages for DDyo7ay. DT and DHe? reaction 
rates as a function of incident deuteron energy in keV. 
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ABSTRACTS IN ENGLISH 


Behaviour of plasma in a rotating magnetie field, A. Lecatow1cz (Institute of Nuclear Research, Warsaw, 


Poland) Nuclear Fusion 1 (1961) 155—159 


An investigation is made of the non-relativistic motion of a charged particle (in a plasma) in an external 
rotating electromagnetic field of the form: Hx, =H, cos (wy/c) cos wt, Hy=H, cos (wx/c) sin wt, Hz 
=constant, H,=Ey=0, Ez=H, [sin (wx/c) cos wt +sin (wy/c) sin wt), with wx/e< 1, wy/e <1. If the con- 
dition —1< (eHz/mem) <{—1+(eH,/mcw)*| is satisfied, then the particle moves away from the z-axis. 
The particle energy is ~H,? e? 7?/(mc? k) in °K where 7 is the average distance from the z-axis. 


The motion of the plasma is then investigated taking into account its proper electromagnetic field. The 
following transport equation is used: nm ¥ + vy -v=ng (E+v x H/c) — yy — nm yp +p. The assumptions 
are: a plasma consisting of equal populations of electrons and deuterons with density ~ 10!°/em®, w < 10!°/sec, 
H,~ 10° G, v< 0.1 c. At t=0 it is assumed that 7’ = 10° °K and v=0 and that the derivatives of the plasma 
density with respect to space variables are negligible compared to other terms of the transport equation. An 
expansion in powers of v/c is used. Zeroth, first and second order approximations are calculated using the 
Laplace transformation. Up to the second order of approximation, the field causes neither a durable change 
in plasma density nor a charge separation. 


Oscillations of four different frequencies appear in the plasma. At a definite frequency of the rotating field there 
appears a resonance phenomenon in which the amplitude of oscillation increases linearly with ¢. At resonance 
the mean energy per ion (in °K) transferred directly to the ionic part of the plasma increases with time as 
follows: (1/192 =*) (e#m/ct M4 k) (H,* Hz?/n,?) t®?. This means, for example, that in an axial field of 104 G 
and a rotating field of amplitude 10° G, the time necessary to provide energy corresponding to 10° °K (dis- 
regarding losses) is ~0.3 sec. 


Axial conduction and radiation losses in a stabilised linear pinch, A. H. pe Borpg, F. A. Haas (English 
Electric Co. Ltd., Nelson Research Laboratories, Stafford, England) Nuclear Fusion 1 (1961) 160—166 


This paper presents a theory of conduction and radiation losses in a linear pinched discharge under steady 
conditions. The model selected is one in which the ohmic heating in a thin skin of current is equated to the 
radiation and axial conduction losses, the discharge being considered at a uniform pressure determined by 
the Bennett relation modified to include the effect of a possible completely trapped axial magnetic field. Formulae 
for temperature and other physical quantities are presented and limiting forms considered in a variety of 
circumstances. The effect of thermoelectric phenomena is considered and conditions under which the treatment 
is likely to be applicable discussed. 


Longitudinal oscillations in a neutralized electron beam with a boundary, M. Yosurkawa (Department of 
Physics, University of Tokyo, Tokyo, Japan) Nuclear Fusion 1 (1961) 167—171 


Longitudinal oscillation of a neutralized electron beam in a strong axial magnetic field is discussed. The beam 
is assumed to be cylindrical and to be confined in a conducting cylinder, or to be planar and to be held between 
two conducting plates. A sufficient and, in some cases, approximately necessary condition for stability is 
obtained. To be stable, the beam current should be below a certain value which depends on the electron velocity, 
the ratio of electron to ion mass and the geometrical dimensions of the beam and the conductor. The validity 
of the approximation is also studied. 


Radiation from a modulated beam of charged particles penetrating a plasma in a uniform magnetic field, 
E. Canossio ( Max-Planck- Institut fiir Physik und Astrophysik, Munich, Federal Republic of Germany) 
Nuclear Fusion 1 (1961) 172—180 


The radiation from a density-modulated beam of ions, which penetrates a plasma perpendicular to a strong 
magnetic field B,, is investigated in two simplified situations: a) the beam is an infinite plane parallel to By, 
and b) the beam is an infinite cylindrical surface parallel to B,, the radius of the cylinder being the gyro-radius 
of the beam particles. This latter beam can be ideally constructed by injecting into a plasma a linear beam, 
modulated at a frequency which is an integral multiple of the gyrofrequency of the beam particles and incident 
in a direction which forms a very small angle with a plane perpendicular to B,. 


In both situations some resonances of the Poynting vector are found. The resonance, which occurs when the 
modulation frequency is equal to the “‘ion-resonance”’ frequency, is specifically investigated, taking into 
account the finite electric conductivity of the plasma. It is shown that, under appropriate conditions, the beam- 
plasma interaction at this resonance becomes very strong. 
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Theory of Cerenkov and cyclotron radiations in plasmas, T. Kruara. O, Aono, R. Sucmmara (Department of 
Physics, University of Tokyo, Tokyo, Japan) Nuclear Fusion 1 (1961) 181—188 


Radiation from a charge gq moving in a helix in magneto-plasma is investigated theoretically. When its speed 
v is much larger than the thermal velocity (m~! kT')'/? of the plasma electrons and the gyration frequency 
is much smaller than the plasma frequency ,, the radiation power from the charge is (q? )/2v) In (v?/m-! kT). 
Cyclotron radiation from an electron with non-relativistic speed decreases to zero as plasma density increases. 
For a positron in a dilute plasma, however, the radiation is strengthened. This strengthened radiation from 
a positron again decreases with increasing *,/my* and becomes zero for w*,/wy? > 2 (my = gyration frequency 
of the plasma electron). Damping of the Cerenkov radiation due to collisions of plasma electrons is also discussed. 


Cyclotron radiation by ions in plasma, V. I. Pisrunovicu, V.D.SHarranov (“J. V. Kurchatov” Institute 
of Atomic Energy, USSR Academy of Sciences, Moscow, USSR) Nuclear Fusion 1 (1961) 189—194 


In this work the intensity of radiation by fast ions in cold plasma is determined. This problem arose in con- 
nection with observations made on the OGRA device of resonance peaks of the electric field strength, i.e. at 
the ion cyclotron frequency and its overtones. Although the abovementioned observations were made outside 
the wave zone, it is natural to expect that the experimentally revealed dependence of the number of observed 
peaks on the density of cold plasma must also be reflected in the intensity of radiation. Calculations showed 
that, with an increase of the ion velocity and of the plasma density, the maximum of the intensity of radiation 
actually shifts towards the high frequencies. 


This shift is analogous to the shift of the maximum of the radiation intensity for synchrotron radiation by an 
electron which has a velocity approaching the velocity of light. In the case of ions in plasma the role of the 
velocity of light is played by the phase velocity of the electro-magnetic waves in the plasma. This velocity 


approaches the Alfvén velocity c,4 = By/\/4= mj nm in the region which is important in the calculations. 
Therefore, the high overtones already become effective at a comparatively small ion velocity, O~c,<c. 


Measurement of the plasma electron temperature in a strong shock wave, T. I. Fitippova, N. V. Frivrppoy, 
V. V. V.P. Vinocrapov (“J. V. Kurchatov” Institute of Atomic Energy, USSR Academy of 
Sciences, Moscow, USSR) Nuclear Fusion 1 (1961) 195—197 

The electron conductivity of deuterium plasma behind a strong shock wave has been measured by the method 
of displaced magnetic flux. At shock-wave velocities from 0.9 x 10’ to 1.25 x 107 em/sec electron temperatures 
of 50 and 90 eV have been found. A comparison between the observed electron temperature and the electron 


temperature computed from the shock-wave velocity has been made. Piezoelectric measurements and com- 
pression of plasma by a strong magnetic field yield exactly the same value of the gas-kinetic pressure of plasma. 


Loss of particles in a pinched discharge in an axial magnetie field, R. K. Jacat (Institute for Fluid Dyna- 
mics and Applied Mathematics University of Maryland, College Park, Maryland, U.S.A.) 

Nuclear Fusion 1 (1961) 198—200 

A caleulation of the loss of deuterons from a pinched current to the wall of the container has been given by 


G. P. Thomson [Phil. Mag. 32 (1958) 886]. This calculation is extended so that account is taken of an axial 
magnetic field. It is found that such a field can materially reduce the particle loss. 
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RESUMES EN FRANCAIS 


Comportement du plasma dans un champ magnétique tournant, A. Lecarowrcz (Institut de recherche nuclé- 
aire, Varsovie, Pologne) Fusion Nucléaire 1 (1961) 155—159 


L’auteur étudie le mouvement non relativiste d’une particule chargée (dans un plasma) dans un champ électro- 
magnétique tournant externe de la forme: Hx =H, cos (wy/c) cos wt, Hy = Hy cos (w2/c) sin wt, Hz =cons- 
tant, Ey=E,y=0, Ez=H, [sin (w2/c) cos wot +sin (wy/c) sin wt], avee wy/e<l. Si la con- 
dition — 1< (eH,z/mem) < |—1-+ (eH,/mem)?} est réalisée, la particule s’éloigne de laxe z. L’énergie 
de la particule dans °K est ~H, e*r/mc? k), oi * est la distance moyenne & partir de l’axe z. 
L’auteur étudie ensuite le mouvement du plasma en tenant compte de son champ électromagnétique propre. 
A cet régard, il utilise ’équation de transport: nm ¥ +vy - ¥=nq (E+v x H/c) — yy — Il admet 
les hypothéses suivantes: le plasma est composé de populations égales d’électrons et de deutérons d’une 
densité de ~10!5/em*, < 10!°/s, Hy~10® G, v < 0,1 ¢; pour t=:0, T= 10% °K et v=0 et les dérivées de 
la densité du plasma par rapport aux variables spatiales sont négligeables si on les compare aux autres termes 
de Péquation de transport. Il utilise un développement en série de puissances de v/c. I] calcule les approxi- 
mations d’ordre zéro, un et deux au moyen de la transformation de Laplace. Jusqu’é lapproximation du 
deuxiéme degré, le champ ne donne lieu ni & un changement durable dans la densité du plasma ni & une 
séparation de charges. 


Des oscillations de quatre fréquences différentes apparaissent dans le plasma. A une fréquence déterminée du 
champ tournant, on constate un phénoméne de résonance dans lequel l’amplitude de Voscillation augmente 
d'une facon linéaire avec le temps. A la résonance, |’énergie moyenne par ion (en °K) transférée directement 
& la partie ionique du plasma augmente avec le temps comme suit: (1/192 =*) e® m/ct M4 k) Hy* Hz*/n,?) &. 
Cela signifie, par exemple, que, dans un champ axial de 10 G et un champ tournant d’une amplitude de 
10° G, le temps nécessaire pour fournir une énergie correspondant & 10° °K (abstraction faite des pertes) est 
~ 0,3 secondes. 


Pertes dues 4 la conduction axiale et au rayonnement dans une striction linéaire stabilisée, A. H. pe Borpe 
et F.A. Haas (English Electric Co. Ltd., Nelson Research Laboratories, Stafford, Royaume-U ni) 
Fusion Nucléaire 1 (1961) 160—166 


Le mémoire expose une théorie concernant les pertes dues & la conducton et au rayonnement dans une 
décharge linéaire pincée, dans des conditions stables. Le modéle choisi est concu de telle sorte que le chauffage 
ohmique dans une mince couche de courant compense les pertes dues au rayonnement et & la conduction axiale, 
la décharge étant considérée & une pression uniforme qui est déterminée par le rapport Bennett modifié de 
maniére & tenir compte de l’effet d’°un champ magnétique axial éventuel entiérement capturé. Les auteurs 
présentent des formules concernant la température et d’autres quantités physiques et examinent les facteurs 
de limitation dans diverses circonstances. Les auteurs examinent enfin l’effet des phénoménes thermoélectri- 
ques et les conditions dans lesquelles le traitement est susceptible d’application. 


Oscillations longitudinales dans un faisceau d’électrons neutralisé et limité, M. Yosuikawa (Département de 
physique, Université de Tokyo, Tokyo, Japon) Fusion Nucléaire 1 (1961) 167—171 


L’auteur examine les oscillations longitudinales d'un faisceau d’électrons neutralisé dans un champ magnétique 
axial intense. Il suppose le faisceau cylindrique et contenu dans une surface cylindrique conductrice, ou de 
forme plane et compris entre deux plaques conductrices. I] réussit & déterminer pour la stabilité une con- 
dition suffisante et, dans certains cas, presque nécessaire. Pour que la stabilité soit réalisée, le courant du faisceau 
doit étre inférieur & une certaine valeur qui dépend de la vitesse des électrons, du rapport entre la masse des 
électrons et celle des ions, ainsi que des dimensions géométriques du faisceau et du conducteur. L’auteur étudie 
également la validité de l’approximation. 


Rayonnements émis par un faiseeau modulé de particules chargées pénétrant dans un plasma dans un champ 
magnétique uniforme, E. Canospsio (Maz-Planck-Institut fiir Physik und Astrophysik, République 
fédérale d’ Allemagne ) Fusion Nucléaire 1 (1961) 172—180 


Les rayonnements émis par un faisceau d’ions & densité modulée, qui pénétre dans un plasma perpendiculaire 
& un champ magnétique intense B,, sont étudiés dans deux situations simplifiées: a) le faisceau est un plan 
infini paralléle & B,; b) le faisceau est une surface cylindrique infinie, paralléle & B,, dont le rayon est le rayon 
de rotation des particules du faisceau. On peut construire idéalement ce dernier faisceau en injectant dans un 
plasma un faisceau linéaire, modulé & une fréquence qui est un multiple entier de la fréquence de giration des 


particules du faisceau, et tombant suivant une direction qui forme un trés petit angle avec un plan perpendicu- 
laire & By. 


Dans les deux situations, on observe quelques résonances du vecteur de Poynting. L’auteur étudie plus parti- 
culiérement la résonance qui se produit lorsque la fréquence de modulation est égale & la fréquence de la 
«résonance ionique», en tenant compte de la conductibilité électrique finie du plasma. I] démontre que, dans 
des conditions appropriées, l’interaction faisceau-plasma devient trés forte & cette résonance. 
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La théorie du rayonnement de Teherenkoy et du rayonnement cyclotronique dans un plasma, T. Kiara, 
O. Aono et R.SuGmara (Départment de physique, Université de Tokyo, Tokyo, Japon) 


Fusion Nucléaire 1 (1961) 181—188 


Les auteurs étudient, du point de vue théorique, les rayonnements émis par une charge q se déplacant en 
= hélice dans des magnétoplasmas. Lorsque la vitesse v est de beaucoup supérieure & la vitesse thermique 
Ao (m-! kT')'/? des électrons du plasma et que la fréquence de rotation est de beaucoup inférieure & la fréquence 
du plasma w,, énergie du rayonnement émanant de la charge est (q? ,?/2v) In (v?/m=! kT’). L’intensité du 
rayonnement cyclotronique émis par un électron ayant une vitesse non relativiste tend vers zéro lorsque la 
densité du plasma augmente. Cependant, dans le cas d’un positon se trouvant dans un plasma dilué, le rayonne- 
ment se trouve intensifié. Ce rayonnement plus intense émis par un positon diminue ensuite & mesure que 
ee w*,/oy” augmente et tombe a zéro pour w*,/oq?> 2 (wy = la fréquence de rotation de l’électron plasmatique). 
Ee Les auteurs examinent en outre l’amortissement du rayonnement de Tcherenkov par suite de collisions 
délectrons plasmatiques. 


Rayonnement cyclotronique des ions dans un plasma, V. 1. Pisrounovircn et V. D. SHarranov (/nstitut 
I. V. Kourtchatov de Vénergie nucléaire, Académie des sciences de l'Union soviétique, Moscou, U RSS) 


Y Fusion Nucléaire 1 (1961) 189—194 


Dans cette étude les auteurs déterminent l’intensité du rayonnement des ions rapides dans un plasma froid. 
Le probléme s’est posé & propos de observation, au moyen d’un dispositif OGRA a des régimes déterminés, 
des pies de résonance de l’intensité du champ électrique 4 partir de la fréquence cyclotronique des ions et de 
ses harmoniques supérieurs. Bien que les observations mentionnées n’aient pas été faites dans la zone de l’onde, 
il était naturel de prévoir que la relation de dépendance, révélée par l’expérience, entre le nombre de pics ob- 
servés et la densité du plasma froid se refléterait aussi dans l’intensité du rayonnement. Les calculs ont montré 


que lorsque la vitesse des ions et la densité du plasma augmentent, l’intensité maximum se déplace effective- 
ment vers les hautes fréquences. 


i Ce déplacement est analogue & celui du maximum d’intensité du rayonnement synchrotronique dia & un électron 
Seay ayant une vitesse proche de celle de la lumiére. Lorsqu’il s’agit d’ions dans un plasma, la réle de la vitesse de 
la lumiére est joué par la vitesse de phase des ondes électromagnétiques dans le plasma; dans la région impor- 
tante pour les calculs, cette vitesse s’approche de la vitesse d’Alfvén: c4 = Bo/(4 = mj no)!?. 


Il s’ensuit que les harmoniques supérieurs élevés se manifestent déja & une vitesse relativement faible des 
ions: vu~ca<c. 


Mesure de la température électronique du plasma dans une forte onde de choe, T. 1. Firrprova, N. V. Fiviepoy, 
V.V. Jourty et V. P. VinoGrapov (Institut I.V. Kourtchatov de Vénergie nucléaire, Académie des sciences 
de l Union soviétique, Moscou, U RSS) Fusion Nucléaire 1 (1961) 195—197 


La conductivité électronique du plasma au deutérium, & larriére d’une forte onde de choc, a été mesurée par 
la méthode du flux magnétique déplacé. On a constaté qu’aux vitesses de propagation de l’onde de choe com- 
prises entre 0,9 x 10° et 1,25 x 107 em/s, correspondaient des températures électroniques de 50 & 90 eV. La 
température électronique observée a été comparée & celle que l’on a caleculée & partir de la vitesse de l’onde 
de choc. Les mesures piézo-électriques et la compression du plasma par un champ magnétique intense ont 
donné une seule et méme valeur pour la pression cinétique des gaz du plasma. 


Perte de particules au cours d’une décharge pincée dans un champ magnétique axial, R. K. Jacer (Jnstitut 
de la dynamique des fluides et de mathématiques appliquées, Université du Maryland, College Park, Maryland, 
Etats-Unis) Fusion Nucléaire 1 (1961) 198—200 


Ge La perte de deutérons entre le courant pincé et la paroi de la chambre a été calculée par G. P. Thomson 
wie [Phil. Mag. 32 (1958) 886]. L’auteur du présent article a perfectionné le calcul de maniére & tenir compte de 


la présence d’un champ magnétique axial. On a constaté qu’un tel champ peut réduire sensiblement la perte 
de particules. 
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AHHOTALIMWA TO-PYCCKU 


[lope enve 1.1a3MbI BO MarHHTHOM nose, A. JIETATOBUY (Mucmumym uccaedoeanuii, 
Bapwaea, [Toavwa) cunte3 1 (1961) 155—159 


Ilpopeneno YaACTHUbI (B M1a3Me) BO BHELIHEM 
nose Bula: Hy= H,cos(wy/c)cos wt, Hy= Hy cos(wx/c)sin wt, H»z=const., 
Ex = Ey =0, Ez = Hoy [sin (x/c) cos wt + sin (my/c) sinwt), ox/c<1, wy/ce<1. ynaosnetsopsetca ycnosue 
—1<(e Hx/mceo)< [—1+ (e Hy/mco)*), To yactuua OT OCH-Z. YaCTHUbI B rpadycax wikKanbt Kenb- 
BuHa CocTaBsinet ~ H, e*r?/(m c? k), rae CpemHee paccTosHHe OT OCH-z. 


ypaBHeHHve mepeHoca: 

+p. ABIAIOTCA: 111a3Ma, COCTOAUIAA M3 PaBHOrO 
QICKTPOHOB DeiTpOHOB Cc NOTHOCTLIO ~ 10!5/cm~, w< 10" cex, Hy~ 108 rayce v< 0,1 c. t= 0 npeanonaraetca, 
T= 108 °K v=0 4TO NPOK3BOAHbIMH MAOTHOCTH M1a3MbI MpOCTPaHCTBeHHbIM MePpeMeCHHbIM MOXHO 
CpaBHeHulo C 4“JieHaMH ypaBHeHHA NepeHoca. Ucnonb3yeTca B NO cTeneHaM v/c. C NOMOLIbIO 
nO BTOpOrO NopsaKa He BbIZbIBACT HH H3MEHECHHA B MIOTHOCTH M1a3Mbl, HH 
3apana. 


B n1a3Me NOABAAIOTCA YETbIPeX PazIM4HbIX YacTOT. YaCTOTe NOABIIA- 
eTCA ABJICHHe Pe3OHAaHCa, MPH KOTOPOM BO3PacTaeT JIMHEMHO CO BpeMeHem. pesoHaHca 
CpeAHAA Ha HOH (B wikanbI Kenbpuna), HENOCpeACTBEHHO B HOHHY!O BO3Pa- 
cTaeT CO BpeMeHeM O6pa3om: (1/192 =*) M4 k) (H,* Hz*/n,”) t?. o3HayaeT, HANPUMeEp, 4TO B AKCHaJIb- 
HOM nosie 10* rayce BO None 10° raycc, BpeMa HEOGXOAMMOEe JIA NOpaKa 
10° B rpagycax wikanbt Kenbauna (6e3 notepb), ~0,3 cek. 


AKcHaJIbHaa NpOBOAHMOCTb NOTepH B A. X. JIE 
BOP, ®. A. XAAC Daexmpux K°, aabopamopuu umenu Heavcona, Cmaddopo, 
Aneausn) Anepubiii cunte3 1 (1961) 160—166 


B HacTOAUIeEM W310%KeHa TEOPHA NPOBOAMMOCTH NOTepb B HHEMHOM CKATOM pa3paue 
B OMHMYeCKOe HarpeBaHHe B TOHKOM TOKA K 
COOTHOWEeHHeM beHHeTTAa, H3MCHEHHBIM C BKIIOYEHHA KaK MOXKHO GosIee BMOPOXeHHOTO B aKCHalIb- 
HOrO MarHuTHOrO POpMysJIbI 110 TeMMepatypbi AIA Apyrux paccMaTpuBaloTcA 
YCNOBUA, MPH KOTOPbIX, NO-BHAMMOMY, MOXKHO MeTOL. 


omdea, Tokuiickuii ynueepcumem, Toxuo, Hnonua) Anepubiii cunte3 1 (1961) 167—171 


NpORONbHOe KONeOaHHe IIEKTPOHHOTO B MOULHOM aKCHaJIbHOM MarHHTHOM Moule. 
MexKLY JBYMA NPOBOAALIMMH Co3maeTcA MOCTaTOYHOe, a B HEKOTOPbIX NOYTH 
BeJIHYHHbI, KOTOPas 3aBHCHT OT CKOPOCTH 3IEKTPOHOB, COOTHOLICHHA WICEKTPOHHOM HOHHOM Macc, a TakKxKe OT reOMeTpH- 


Ny4Ka YACTHI, B 11a3My B MarHHTHOM 
nose, E. KAHOBBUO (®u3zuyeckuii u acmpoguzuyeckuii uncmumym Makca-Ilaauxa, Mionxen, 
cunte3 1 (1961) 172—180 


Uccneayetca OT Ny4Ka HOHOB MOLYJHPOBaHHOH NIOTHOCTH, KOTOPbIi NPOHHKAaeT Yepe3 
CHJIBHOMY MarHHTHOMY By B : a) NYYOK ABIAeTCA GecKOHeEYHOH Napasieb- 
Ho B,; 6) ny4oK GeckOHe4YHOK NOBEPXHOCTbH, NapaswienbHoK By, pannyc KOTOpOK aBAeTCA 
PaMycoM BpallleHHsA YaCcTHL NyyKa. DTOT NyYOK MOxeT CO32aH B 
ny4ka, NpH YacTOTe, KOTOPaA ABJIACTCA LICJIbIM KPaTHbIM YAaCTOTHI BPallleHHA YaCTHL My4Ka, 
B KOTOPOe CO32aeT Malibii C By. 


B o60nx Hai HECKONbKO pe3oHaHcos BeKTopa CneunasbHo pe3oHaHC, 
Mbiii B TOM Clly4ae, KOr a MOAYIAUMOHHAA YACTOTa PaBHa “HOH-Pe3OHAaHCHOM» YaCTOTe, C YYeTOM WIeKTpO- 
3TOM Pe3OHAHCe CTAHOBUTCA OYCHb CHJIbHbIM. 
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Teopua H3J1y4eHHA B T. KAXAPA, O. AOHO, P. CYTUXAPA ( 


yeckui omdea, TokuiicKuii yuueepcumem, Toxuo, Anonus) Anepubii cunte3 1 (1961) 181 —188 


Teoperuyecku 3apAda g, NO B Korga ero CKOpocTh v 
3HAYMTEIBHO CKOpocTb (m~!kT)!/* sneKTpOHOB M1a3Mbl, a YACTOTA 3HAYHTEJIBHO 
4aCTOTHI M1a3MbI MOULHOCTb paBHa (q?~,"/2v) In (v?/m! kT). 2eKTpoHa 
NO3KTPOHA B pa3pexeHHOH M1a3Me HHTCHCHBHOCTh ITA yBesIMYeHHaA MHTCHCHBHOCTh 


PaccmaTpuBaeTca TakxKe 3aTYXaHHe YEPEHKOBCKOTO BCIEACTBHE CTOJIKHOBEHHH 
QsIEKTPOHOB. 


Lluk10Tponnoe HOHOB B B. U. B. J]. LIA®@PAHOB (Hucmumym amomnuoii 
gnepeuu umenu H. B. Kypuamoea Akademuu Hayx CCCP, Mocxea, CCCP) 


Anepubii cunte3 1 (1961) 189-—194 


B paOote onpesesieHa HHTCHCHBHOCTh H3J1y4eHHA GbICTpbIX HOHOB B XOJIONHOM 3a1a4a BO3HMKJIA B CBA3H C 
Ha ycTaHoBKe Orpa B PC3OHAHCHbIX MMKOB HalpAKEHHOCTH 
Ha WMKJIOTPOHHOM HOHHOK YacTOoTe ee OGepTOHAaX. XOTA yKa3aHHble HAGMIONCHHA HE B BOJIHOBOH 30He, 
€CTeCTBEHHO YTO OOHAPy2xKeHHAaA B IKCMEPHMEHTE MMKOB OT MAOTHOCTH 
HOM M1a3Mbl HalTH CBOe OTPaxkeHHe B HHTCHCHBHOCTH PacueTbI MOKa3a/IM, YTO NPM 
CKOPOCTH HOHOB H MJIOTHOCTH MAKCHMYM HHTCHCHBHOCTH CMellaeTCA B CTOPOHy 
BbICOKHX 4acTOT. 

STO CMelleHHe CMeLICHHIO MaKCMMYMa MHTeCHCHBHOCTH MPH CHHXPOTPOHHOM WIEKTPOHA, 
CKOpOCThb, K CKOPOCTH cBeTa. B HOHOB B 11J1a3Me POsIb CKOPOCTH CBeTa HrpaeT CKOPOCTL 
BOJIH B KOTOPaA B pacyeToB OGsacTH K CKOpOCTH 


cA = By|(4 = mj Tlostomy sbicoxue oGepToHbI NpOABIAIOTCA yxKe NPM CpaBHUTebHO CKOPOCTH HOHOB 
U~ca<c. 


Vismepenne TeMNepaTypbi 1.1a3MbI B BoaHe, T. H. B. 
nos, B. B. KKYPHH, B. IT. BAHOrPAJLOB (Mucmumym amomnoii snepeuu umenu B. Kypuamoea Axademuu 
Hayx CCCP, Mocxea, CCCP) Anepubii cunte3 1 (1961) 195— 197 


208 


MeToOM BbITeCHEHHOrO MarHHTHOrO MOTOKa M3MepeHa IIEKTPOHHAA MPOBOAMMOCTh 3a CHIbHOI 
ckopoctax yoapHoi OT 0,9- 107 no 1,25-107 cm/cek HaiiaeHa Temmepatypa B 50x 
90 eV. Ilpospeneno cpaBHenue WeKTPOHHOK TeMMepaTypbI C MO CKOPOCTH 


H3IMCPCHHA H CHATHC CHJIBHbIM M@rHUTHbIM MOJICM OHO TO Ke 3HAYCHHE TAIOKHHCTHYCCKOTO 
AABNCHHA 


Tlorepa 4uactuy B B OCeBOM MarHHTHOM nove, P. K. (Mucmumym dunamuxu 
Kocmeti u npukAadnoii Mamemamuku Mepuasndckoeo yuueepcumema, Koaaed# Tlapx, Mepuagnd, CIA) 


Anepupiii cunte3 1 (1961) 198 —200 


I. TI. Tomconom [Phil. Mag. 32 (1958) 886] coenaHbI pacyeTbI M3 NHHY-TOKa K CTeHKe cocyia. 
STH pacieTbI pacnpocTpaHeHbl WIA yueTa BIIMAHHA OCeBOrO MarHHuTHOrO NosIA. HaiigeHo, 4TO TaKOe NOJIC MOXKET CyLecT- 
BeHHO YMCHbLUMTb MOTeplO YaCTHLL. 
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RESUMENES EN ESPANOL 


Comportamiento del plasma en un campo magnético giratorio, A. Lecatowicz (Instituto de Investigaciones 
Nucleares, Varsovia, Polonia) Fusion Nuclear 1 (1961) 155—159 


El autor estudia el movimiento no relativista de una particula cargada (en un plasma) en un campo electro- 
magnético giratorio externo, que responde a la forma: H, = Hy, cos (w y/c) cos wt, Hy = Hy cos (w x/c) sen wt, 

z=constante, H,=E,y=0, Ez=H, [sen (w2z/c) cosmt+sen (wy/c) sen oft), siendo wy/e<l. 
Si se cumple la condicién — 1 < (eHx/me w) <[—1+ (eHy/mew)*], la particula se aleja del eje z. La energia 
de la particula, expresada en °K, es~H, e? ¥?/(mc? k), donde 7 es la distancia media al eje z. 


Seguidamente, el autor estudia el movimiento del plasma teniendo en cuenta su campo electromagnético 
propio. Utiliza la siguiente ecuacién de transporte: nmv +vy-v=ng (E+v x H/c) —- yy —nmyy +p, 
formulando las siguientes hipdtesis basicas: el plasma consta de ntimero igual de electrones y deuterones con 
una densidad de ~10!5/em', w < 10!°/s, Hy ~108 G, v< 0,1 c. Para t=0, se supone que = 108 °K y v=90 
y que las derivadas de la densidad del plasma con respecto a las variables espaciales son despreciables en com- 
paracion con los demas términos de la ecuacién de transporte. Desarrolla v/c en serie de potencias y calcula 
las aproximaciones de orden cero, primero y segundo, usando la transformacién de Laplace. Hasta la aproxi- 
macion de segundo orden, el campo no causa cambios duraderos en la densidad del plasma ni provoca ninguna 
separaciOn de cargas. 


Se observan en el plasma oscilaciones de cuatro frecuencias diferentes. Cuando la frecuencia del campo giratorio 
alcanza un valor determinado, se manifiesta un fenémeno de resonancia en el que la amplitud de la oscilacion 
crece en forma directamente proporcional al tiempo. Durante la resonancia, la energia media por ion (expresada 
en °K) que se transfiere directamente a la parte iénica del plasma aumenta con el tiempo segtin la siguiente 
funcién: (1/192 =?) (e#m/c* M4 k) H*, H*z/n,?) t?. Esto significa que, por ejemplo, en un campo axial de 10* G 
y un campo giratorio de una amplitud 10° G el tiempo necesario para suministrar la energia correspondiente 
a 10° °K (despreciando las pérdidas) es del orden de los ~0.3 s. 


Pérdidas debidas a la conduecién axial y a la radiacién en una estriccién lineal estabilizada, A. H. pe Borpe, 
F. A. Haas (English Electric Co. Lid., Nelson Research Laboratories, Stafford, Reino Unido) 
Fusion Nuclear 1 (1961) 160—166 


Los autores exponen una teoria acerca de las pérdidas debidas a la conduccién y a la radiacién en una descarga 
con estriccién lineal en condiciones de estabilidad. Han elegido un modelo en el que se admite que el calenta- 
miento éhmico en una capa delgada de corriente compensa las pérdidas por radiacién y por conduccién axial, 
considerdndose que la descarga se verifica a una presién uniforme que esta determinada por la relacién de 
Bennett, modificada a fin de tener en cuenta el efecto de un posible campo magnético axial totalmente captu- 
rado. Presentan formulas referentes a la temperatura y otros pardmetros fisicos y consideran factores de 
limitacién en diversas circunstancias. 

Por ultimo, estudian el efecto de los fenémenos termoeléctricos y analizan las condiciones en que el tratamiento 
es susceptible de aplicarse. 


Oscilaciones longitudinales en un haz electrénico neutralizado y limitado, M. Yosurkawa (Departamento de 

Fisica, Universidad de Tokio, Tokio, Japén) Fusi6n Nuclear 1 (1961) 167—171 

El autor estudia las oscilaciones longitudinales de un haz de electrones neutralizado en un campo magnético 

axial de gran intensidad. Supone que el haz es cilindrico y se encuentra encerrado en un cilindro conductor o 

bien que es de forma plana y esta intercalado entre dos placas conductoras. Determina la condicién suficiente 

y, en algunos casos, cuasi necesaria, para que se mantenga la estabilidad. Para aleanzar esa estabilidad, es 

preciso que la corriente del haz sea inferior a cierto valor, que depende de la velocidad de los electrones y de 

la razén masa electrénica/masa idnica, asi como de las dimensiones geométricas del haz y del conductor. El 
autor estudia también la validez de la aproximacidén. 


Radiacién emitida por un haz modulado de particulas cargadas al penetrar en un plasma en un campo magné- 
tico uniforme, E. Canosppio (Maz-Planck-Institut fiir Physik und Astrophysik, Munich, Repiblica Federal 
de Alemania) Fusion Nuclear 1 (1961) 172—180 


La radiaci6n emitida por un haz idnico de densidad modulada al penetrar en un plasma perpendicular a un 
campo magnético intenso B, es estudiada en dos casos simplificados: a) cuando el haz es un plano infinito 
paralelo a B, y b) cuando el haz forma una superficie cilindrica infinita paralela a B,, cuyo radio es el radio 
de giro de las particulas del haz. Este ultimo puede obtenerse idealmente inyectando en un plasma un haz lineal 
modulado a una frequencia que sea un multiplo entero de la frecuencia de giro de las particulas del haz y que 
incida en una direccién que forme un angulo muy pequefio con un plano perpendicular a By. 


En ambos casos se observan algunas resonancias del vector de Poynting. El autor estudia especialmente la 
resonancia que se produce cuando la frecuencia de modulacién es igual a la frecuencia de “‘resonancia idnica”’, 
teniendo en cuenta la conductividad eléctrica finita del plasma. Demuestra que en condiciones favorables la 
interaccién haz-plasma se hace muy fuerte para esa resonancia. 
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Teoria de la radiacién de Cerenkov y de la radiacién ciclotrénica en los plasmas, T. Kraara, O. Aono, 
ae R. Sucmmara (Departamento de Fisica, Universidad de Tokio, Tokio, Japén) 


Fusion Nuclear (1961) 181—Iss8 


Los autores estudian desde el punto de vista tedrico, la radiacién emitida por una carga q que describe una 
trayectoria helicoidal en un magnetoplasma. Cuando su velocidad v es mucho mayor que la velocidad térmica 
(m-! k7’)'/? de los electrones del plasma y la frecuencia de rotaci6n mucho menor que la frecuencia w, del 
plasma, la energia de la radiacién procedente de la carga viene dada por la expresién (q? w,?/2v) In (v?/m kT). 
La radiacién ciclotrénica emitida por un electrén de velocidad no relativista decrece hasta cero a medida que 
la densidad del plasma aumenta. Sin embargo, esa radiacion se intensifica para un positrén en un plasma diluido. 
Esta radiacion intensificada del positrén disminuye seguidamente al aumentar el cociente ,?/my? y se anula 
para 2 (my = frecuencia de giro del electrén plasmatico). Los autores discuten también la amorti- 
Gamat guaciOn de la radiacién de Cerenkov debida a las colisiones de electrones plasmaticos. 


Radiacién ciclotrénica de los iones en un plasma, V. 1. Pisrunovicn, V.D.SHarranov (Jnstituto de 
Energia Atomica “I. V. Kurchatov”, Academia de Ciencias de la U.R.S.S., Moseti, U.R.S.S.) 


Fusion Nuclear 1 (1961) 189—194 


Los autores determinaron la intensidad de radiacién de los iones rdpidos en un plasma frio. El problema se 
planted con motivo de observaciones, efectuadas en el aparato OGRA, de los maximos de resonancia de la 
intensidad de un campo eléctrico en regimenes determinados, a saber, a la frecuencia ciclotr6énica de los iones 
y sus armonicos. Aunque las citadas observaciones no se realizaron en la zona ondulatoria, es l6gico suponer 
que la relacién de dependencia entre el ntimero de maximos registrados y la densidad del plasma frio, hallada 
a experimentalmente, se reflejaria también en la intensidad de radiacién. Los calculos indicaron que, al aumentar 
Bets la velocidad de los iones y la densidad del plasma, el maximo de intensidad se desplaza, efectivamente, en el 
sentido de las frecuencias elevadas. 


Este desplazamiento es andlogo al desplazamiento del maximo de intensidad de la radiacién sincrotrénica 
de un electr6n animado de una velocidad préxima a la de la luz. En el caso de los iones plasmaticos, el papel 
de la velocidad de la luz lo desempenia la velocidad fasica de las ondas electromagnéticas en el plasma; en la 
iY. regién que reviste importancia para los cdlculos, esta velocidad se aproxima a la velocidad de Alfvén: 
uate cx = By/V/4e mj no. Por consiguiente, los arménicos de orden superior se manifiestan ya en correspondencia 
ae con una velocidad idnica relativamente baja, a saber, v~ca<c. 


e Medicién de la temperatura electrénica de un plasma en una onda de choque de gran intensidad, T. 1. Frtrprova, 
eit N. V. Frurepov, V. V. Yurty, V. P. VinoGrapov (Instituto de Energia Atémica “I. V. Kurchatov’ Academia 
ae de Ciencias de la U.R.S.S., Moscti, U.RS.S) Fusion Nuclear 1 (1961) 195—197 


F Aplicando el método de desplazamiento del flujo magnético, los autores han medido la conductividad electr6nica 
Pts de un plasma de deuterio detras de una fuerte onda de choque. Comprobaron que a velocidades de propagacién 
a de 0,9 x 107 y 1,25 x 10? em/s correspondian temperaturas electrénicas de 50 y 90 eV, respectivamente. Se 
eh establecié una comparacion entre la temperatura electronica registrada y la calculada a partir de la velocidad 
i de la onda de choque. Las mediciones piezoeléctricas y la constriccién del plasma mediante un campo magnético 

. intenso condujeron a un mismo valor para la presién cinética del gas plasmiatico. 


Pérdida de particulas en una desearga constrefida en un campo magnético axial, R. K. Jacer (Instituto 
de Dindmica de los Fliidos y Matemdticas Aplicadas, Universidad de Maryland, College Park, Mary- 
land, Estados Unidos) Fusion Nuclear 1 (1961) 198—200 


La pérdida de deuterones entre una corriente constrenida y la pared del recipiente fue caleulada por G. P. Thom- 
son [ Phil. Mag. 82 (1958) 886]. El autor del presente trabajo extiende los limites de dicho calculo a fin de tener 
en cuenta el efecto de un campo magnético axial. Llega a la conclusién de que tal campo es capaz de reducir 
sensiblemente la pérdida de particulas. 


: 
7 
Ve 
‘ 
P 
hy 
210 | 
| 
Py 


NO, 
188 


una 
nica 

del 
$y 
que 
ido, 
orti- 


de 
194 


A se 
e la 
nes 
yner 
ada 
ntar 
n el 


nica 
Apel 
n la 
én: 
ncia 


VA, 
mia 
197 


nica 
. Se 
dad 
tico 


tuto 
ury- 
200 


om- 
ucir 


Vol. 1, No. 2 ERRATA ET ADDENDA Tom 1, Bbinyck 2 


Page 79 


Page 134 


Figure 2: delete the factor 10? in the ordinates label. Similarly, delete the factor 10 in Figure 3. 


The approximate expression in Eq. (21) and Eq. (22) are incorrect even to first order, inasmuch as 
the quantity y, assumed constant, contains n,. A useful expression for the upper critical current, 
where background gas density N, can be neglected, comes from Eq. (23): 
(Ty) >= Pox v 
It has been pointed out by I. N. Golovin (P. R. Bell, private communication) that a sufficiently 
high initial pumping speed, 9), can thus reduce the upper critical current significantly. 
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TRANSLATIONS 


A limited number of copies of the translated texts (without figures and equations) of articles which have appeared 
in “Nuclear Fusion” stet available for subscribers to the journal or for those official agencies, exchange centers, govern- 
mental libraries, ete. who receive free copies of the journal. 


If you are subscriber, or an official recipient, and wish to receive a copy of one or more translated texts, address your 


request to “The Editor, Nuclear Fusion, International Atomic Energy Agency, Kaerntnerring 11, Vienne I, Austria”, 


Include your address and the ‘TT’? number (see list below) of the translated texts which you wish. 


TRADUCTIONS 


Un nombre limité d’exemplaires des traductions (sans figures ni équations) des articles parus dans «Fusion nucléaire» 
seront prochainement & la disposition des abonnés et des institutions officielles, centres d’échange, bibliothéques 
nationales, ete. qui regoivent des exemplaires gratuits de la revue. 

Si vous étes abonné & la revue ou si elle vous parvient a titre officiel, et si vous désirez un exemplaire d’une ou de 
plusieurs traductions, adressez votre demande au Rédacteur de «Fusion nucléaire», Agence internationale de I’énergie 
atomique, Kaerntnerring 11, Vienne I, Autriche. Indiquez votre adresse et le numéro «TT» (voir liste ci-dessous) 
des traductions qui vous seraient utiles. 


TIEPEBOJIbI 


B cCKOpOM BpeMeHH OFpaHHYeHHOe YHCIIO IK3EMMIAPOB NMepeBeAeHHbIX TeKCTOB CTaTei (6e3 4YepTexeli (hOpMy.1), NOABHBLUMXCA 
B *KypHase CHHTe3", NOANMCYMKAaM WIM TEM LEHTpam n0 


HbIX TeKCTOB, OOpaulaiitecb no anpecy: «The Editor, Nuclear Fusion, International Atomic Energy Agency, Kaerntner- 
ring 11, Vienna I, Austria». Ykaxute Baul HOMep «TT» mepeBeneHHbIX TekcTOB (cM. HH%Ke) KOTOPbIC BbI 


TRADUCCIONES 


En breve se dispondraé de un ntimero limitado de ejemplares de los textos traducidos (sin cifras ni ecuaciones) de 
los articulos aparecidos en «Fusién Nuclear» para las personas suscritas a la revista o para los organismos oficiales, 
centros de intercambio, bibliotecas ptblicas, etc., que reciben ejemplares gratuitos. 

Si esta usted suscrito a la revista, 0 la recibe a titulo oficial, y desea recibir un ejemplar de uno 0 mas textos traducidos, 


pidalos a la direccién siguiente: «Redactor de Fusién Nuclear, Organismo Internacional de Energia Atémica, 


Kaerntnerring 11, Viena I, Austria». Indique su direccién y el ntiimero «TT» de los textos traducidos (véase la 
lista a continuacién) que desee. 


* Translated into 
(1960) 3—41 R 
(1960) 42—46 R 
(1960) 47—5: E 
(1960) 47—53 R 
(1960) 54—61 R 
(1960) 62—63 R 
(1961) 82—100 E 
(1961) 121—124 E 
(1961) 189—194 E 
(1961) 195—197 E 


TT-1 Nuclear Fusion 
TT-2 Nuclear Fusion 
TT-3 Nuclear Fusion 
TT-4 Nuclear Fusion 
TT-5 Nuclear Fusion 
TT-6 Nuclear Fusion 
TT-7 Nuclear Fusion 
TT-8 Nuclear Fusion 
TT-9 Nuclear Fusion 
TT-10 Nuclear Fusion 


* EK — English, R — Russian 
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